Tổng quan nghiên cứu

Vật liệu nhôm ôxít (Al2O3) là một trong những hệ gốm kỹ thuật quan trọng hàng đầu với nhiệt độ nóng chảy lên tới 2327 K (xấp xỉ 2054 độ C), sở hữu độ bền cơ học, độ bền nhiệt vượt trội cùng khả năng ứng dụng rộng rãi trong quang học, vi điện tử, y sinh và vật liệu xúc tác. Trong chu trình chuyển pha cấu trúc, Al2O3 thủy tinh và lỏng đóng vai trò là tiền chất then chốt để hình thành pha không bền gamma-Al2O3 trước khi chuyển hóa hoàn toàn sang pha bền nhiệt động alpha-Al2O3 (corundum). Tuy nhiên, việc khảo sát thực nghiệm đối với hệ Al2O3 ở trạng thái lỏng và vô định hình luôn đối mặt với rào cản kỹ thuật lớn do nhiệt độ cực cao gây ra phản ứng hóa học giữa mẫu và vật chứa, dẫn đến nguy cơ lẫn tạp chất và tạo ra sự sai lệch đáng kể giữa các kết quả công bố. Các nghiên cứu thực nghiệm trước đây ghi nhận khoảng cách liên kết Al-O dao động từ 1,80 Å đến 2,02 Å, trong khi số phối trí trung bình của nhôm biến thiên từ 4,1 đến 4,8 mà chưa có mô hình giải thích thấu đáo.

Để giải quyết mâu thuẫn này, luận văn thạc sĩ chuyên ngành Vật lý chất rắn của tác giả Nguyễn Hải Yến, dưới sự hướng dẫn khoa học của PGS.TS Phạm Hữu Kiên tại Trường Đại học Sư phạm thuộc Đại học Thái Nguyên (bảo vệ năm 2022, đạt độ tương đồng 23% qua kiểm định Turnitin), đã tập trung nghiên cứu định lượng mối tương quan giữa tính chất cấu trúc và các đơn vị phối trí trong hệ ôxít Al2O3 thủy tinh bằng phương pháp mô phỏng tính toán vi mô. Phạm vi nghiên cứu được triển khai toàn diện trên các mô hình nguyên tử có quy mô từ 3000 đến 5500 hạt, quét qua dải nhiệt độ từ 300 K đến 6000 K và dải áp suất nén từ 0 đến 30 GPa (tương ứng với mật độ khối lượng từ 2,85 g/cm3 đến 3,82 g/cm3). Nghiên cứu mang ý nghĩa học thuật và thực tiễn to lớn khi thiết lập thành công các phương trình giải tích tường minh liên kết phân bố góc liên kết với tỷ phần đa diện phối trí, đồng thời xác lập mô hình hai pha giúp nâng cao độ chính xác dự báo cấu trúc vật liệu lên mức sai số dưới 2% so với thực nghiệm nhiễu xạ.

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 nền tảng cơ học thống kê cổ điển kết hợp với định luật chuyển động Newton để mô tả quỹ đạo động học của hệ đa hạt trong không gian pha. Tương tác vi mô giữa các ion Al3+ và O2- được mô hình hóa thông qua thế thế năng cặp hiệu dụng Born-Mayer kết hợp tương tác Coulomb tầm xa. Trong đó, tương tác lực đẩy bán thực nghiệm có các tham số đặc trưng: năng lượng tương tác đẩy B(Al-O) = 1479,86 eV với bán kính suy giảm R(Al-O) = 3,4483 Å^-1; năng lượng B(O-O) = 1500,00 eV với R(O-O) = 3,4483 Å^-1; điện tích hiệu dụng quy ước của nhôm là +3e và ôxy là -2e.

Nhằm triệt tiêu sự phân kỳ và tối ưu hóa thời gian tính toán cho lực Coulomb tầm xa vốn chiếm hơn 80% tài nguyên xử lý, nghiên cứu áp dụng kỹ thuật gần đúng Ewald-Hansen. Kỹ thuật này phân tách thế tĩnh điện thành tổng thành phần hội tụ nhanh trong không gian thực và không gian đảo thông qua ô mạng tuần hoàn 3x3. Khung lý thuyết của luận văn tập trung vào 4 khái niệm cấu trúc cốt lõi:

  • Hàm phân bố xuyên tâm g(r): Xác định xác suất tìm thấy nguyên tử tại khoảng cách r, phản ánh trật tự tầm gần và tầm trung.
  • Thừa số cấu trúc S(q): Biến đổi Fourier của hàm tương quan cặp, đối sánh trực tiếp với phổ tán xạ tia X và neutron.
  • Đơn vị phối trí đa diện AlOx (với x = 4, 5, 6, 7): Các khối hình học cơ sở cấu thành mạng vô định hình.
  • Phân bố góc liên kết O-Al-O và Al-O-Al: Đặc trưng cho mức độ biến dạng mạng và kiểu liên kết chung đỉnh, chung cạnh giữa các khối đa diện.

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

Nghiên cứu sử dụng phương pháp mô phỏng động lực học phân tử (ĐLHPT) cổ điển làm công cụ chủ đạo. Cỡ mẫu nghiên cứu được chuẩn hóa chính xác ở quy mô 3000 nguyên tử, bao gồm 1200 nguyên tử Al và 1800 nguyên tử O, bảo đảm tuyệt đối tỷ lệ hợp thức 2:3 của phân tử Al2O3. Phương pháp chọn mẫu ban đầu là phân bố ngẫu nhiên các hạt vào ô mô phỏng hình lập phương có kích thước xác định sao cho mật độ khối lượng khởi điểm đạt 2,85 g/cm3. Điều kiện biên tuần hoàn 3 chiều (3D periodic boundary conditions) được thiết lập xuyên suốt quá trình tính toán nhằm loại bỏ triệt để các hiệu ứng bề mặt và mô phỏng chính xác trạng thái của vật liệu khối.

Lý do lựa chọn phương pháp ĐLHPT kết hợp thế Born-Mayer xuất phát từ ưu thế vượt trội về mặt quy mô tính toán so với mô phỏng nguyên lý ban đầu (ab initio) vốn bị giới hạn ở hệ dưới 300 nguyên tử, đồng thời vượt trội hơn mô phỏng Monte Carlo ở khả năng theo dõi liên tục tiến trình động học theo thời gian thực. Phương trình vi phân chuyển động được giải số bằng thuật toán tích phân Verlet với bước thời gian siêu nhỏ dt = 0,45 femtô giây (fs), đảm bảo độ hội tụ năng lượng và tính ổn định pha tối ưu.

Quy trình nhiệt học và nén ép mẫu diễn ra nghiêm ngặt qua nhiều giai đoạn: Hệ ban đầu đạt cân bằng hồi phục tại 0 K sau 8000 bước, nung nóng cưỡng bức lên 5000 K trong 50.000 bước để phá vỡ hoàn toàn cấu trúc tinh thể, sau đó làm nguội dần qua các mốc 2000 K, 1000 K và đạt 300 K với tốc độ làm nguội ổn định 0,5 K/ps. Sau 100.000 bước hồi phục ở áp suất thường, mẫu tiếp tục được nén đẳng nhiệt tại 300 K qua 750.000 bước để tạo ra 6 mô hình mật độ từ 2,85 g/cm3 đến 3,82 g/cm3. Đối với nghiên cứu áp suất cao, 12 mô hình ở dải áp suất 5, 10, 15, 20, 25, 30 GPa tại hai mốc nhiệt độ 300 K và 3000 K được hồi phục cân bằng trong tập hợp vi chính tắc (NVE) qua 50.000 bước tính toán. Bán kính ngắt cấu trúc Al-O được cố định ở giá trị 2,40 Å tương ứng với điểm cực tiểu đầu tiên của hàm phân bố xuyên tâm.

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 ĐLHPT đã cung cấp bức tranh chi tiết ở cấp độ nguyên tử về sự biến đổi mạng Al2O3 dưới tác động của mật độ và áp suất nén, thể hiện qua 4 phát hiện quan trọng:

  1. Chuyển dịch số phối trí trung bình: Khi mật độ khối lượng tăng từ 2,85 g/cm3 lên 3,82 g/cm3, số phối trí trung bình Z(Al-O) tăng mạnh từ 4,20 lên 5,08 (mức tăng ròng 20,95%), trong khi số phối trí Z(O-Al) tăng từ 2,79 lên 3,35. Kết quả này phản ánh xu hướng các nguyên tử xếp chặt hơn khi không gian tự do bị thu hẹp.

  2. Sự tái cấu trúc các đơn vị đa diện phối trí: Tỉ phần của cấu trúc tứ diện AlO4 suy giảm mạnh từ 79,0% ở mật độ 2,85 g/cm3 xuống chỉ còn 19,0% ở mật độ 3,82 g/cm3 (mức sụt giảm lên tới 60,0 điểm phần trăm). Ngược lại, đơn vị phối trí 5 đỉnh AlO5 tăng mạnh từ 19,9% lên chiếm ưu thế tuyệt đối ở mức 54,8%, và đơn vị bát diện AlO6 tăng vọt từ 0,9% lên 26,5%. Đồng thời, tỷ lệ liên kết cầu ôxy OAl2 giảm từ 22,0% xuống 1,5%, trong khi liên kết OAl4 tăng từ 2,5% lên 38,0%.

  3. Thu hẹp phân bố góc liên kết: Góc liên kết trung bình O-Al-O giảm liên tục từ 106,5 độ xuống 88,5 độ (thu hẹp 18,0 độ), tiệm cận dần về góc 90 độ đặc trưng của khối bát diện đều. Trong khi đó, góc liên kết cầu Al-O-Al duy trì tương đối ổn định trong khoảng 120,5 độ đến 121,3 độ trước khi giảm nhẹ xuống 115,2 độ ở mật độ cao nhất. Luận văn đã thiết lập thành công phương trình giải tích biểu diễn hàm phân bố góc liên kết gAl(theta) và gO(theta) thông qua tổ hợp tuyến tính tỷ phần các đơn vị AlOx và OAly với độ khớp gần như tuyệt đối.

  4. Xác lập trạng thái hai pha và áp suất chuyển pha: Cấu trúc Al2O3 tồn tại đồng thời hai pha riêng biệt: pha mật độ thấp (MĐT) hình thành từ các cụm tứ diện AlO4 kết nối qua liên kết OAl2 có mật độ riêng 2,54 g/cm3; và pha mật độ cao (MĐC) cấu thành từ các cụm AlO5, AlO6, AlO7 liên kết qua OAl2 và OAl3 có mật độ riêng 4,444 g/cm3. Điểm cắt chuyển pha của đơn vị AlO5 diễn ra tại áp suất 7,5 GPa và 16,0 GPa đối với pha lỏng (3000 K), nhưng dịch chuyển lên ngưỡng áp suất cao hơn nhiều là 12,0 GPa và 22,5 GPa đối với pha thủy tinh (300 K).

Thảo luận kết quả

Khi khảo sát sự co ngắn liên kết, dữ liệu tính toán chỉ ra rằng khoảng cách cặp O-O giảm từ 2,77 Å xuống 2,61 Å (giảm 5,78%) và khoảng cách Al-Al giảm từ 3,15 Å xuống 3,03 Å (giảm 3,81%). Trái lại, khoảng cách liên kết Al-O lân cận gần nhất hầu như không đổi, dao động trong biên độ hẹp từ 1,77 Å đến 1,80 Å (sai lệch dưới 1,7%). Điều này chứng minh rằng trật tự tầm trung trong mạng thủy tinh Al2O3 chịu tác động mạnh mẽ của áp suất nén hơn hẳn so với trật tự tầm gần, do các liên kết cầu góc Al-O-Al dễ dàng bị uốn cong và nén ép hơn độ dài liên kết hóa học trực tiếp.

Những biến đổi tinh tế này có thể được trình bày và đối sánh trực quan rất hiệu quả thông qua đồ thị hàm phân bố xuyên tâm g(r) và phổ thừa số cấu trúc tán xạ tia X S(q). Trên đồ thị S(q), khi áp suất nén tăng từ 0 đến 30 GPa, đỉnh cực đại thứ nhất tại vùng góc tán xạ thấp dịch chuyển rõ rệt từ vị trí q1 = 2,1 Å^-1 sang 2,4 Å^-1, phản ánh sự co cụm của mạng nguyên tử. Đáng chú ý, đỉnh phụ thứ hai xuất hiện ở mức áp suất 5 GPa đối với Al2O3 lỏng nhưng chỉ xuất hiện khi áp suất đạt tới 20 GPa đối với Al2O3 thủy tinh. Hiện tượng này minh chứng rằng mạng chất lỏng có độ linh động cấu trúc cao hơn và nhạy cảm với áp suất chuyển pha sớm hơn trạng thái vô định hình rắn.

Đặc biệt, nghiên cứu động học khuếch tán trong Al2O3 lỏng làm sáng tỏ hiện tượng khuếch tán dị thường: Hệ số khuếch tán không suy giảm đơn điệu theo áp suất mà xuất hiện điểm cực trị cục bộ tại 15 GPa đối với nguyên tử nhôm và 21 GPa đối với nguyên tử ôxy. Nguyên nhân là do tại ngưỡng áp suất này, các liên kết Al-O trong khối đa diện AlO5 không bền bị kéo căng và đứt gãy hàng loạt để tái cấu trúc sang mạng bát diện AlO6, làm suy yếu tương tác giam giữ tạm thời. So sánh với các công bố thực nghiệm nhiễu xạ tia X và neutron của Lamparter (1996), Ansell (1997) và Skinner (2008), các đường cong cấu trúc tính toán từ mô hình hoàn toàn tương thích về vị trí các đỉnh cực đại (1,80 Å cho Al-O và 2,80 Å cho O-O), khẳng định tính đúng đắn và độ tin cậy vượt bậc của luận văn.

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

Dựa trên các phát hiện cấu trúc và quy luật chuyển pha hai pha đã xác lập, luận văn đưa ra 4 khuyến nghị mang tính ứng dụng cụ thể:

  1. Chuẩn hóa quy trình thiêu kết vật liệu gốm kỹ thuật cao: Các viện nghiên cứu vật liệu và doanh nghiệp luyện kim cần áp dụng quy trình ép nóng đẳng tĩnh ở dải áp suất tối ưu từ 5 đến 10 GPa trong chu kỳ thử nghiệm 12 đến 18 tháng tới. Mục tiêu là kiểm soát tỷ phần đơn vị tứ diện AlO4 duy trì trên 60%, giúp triệt tiêu các lỗ xốp vi mô và nâng cao độ bền uốn của gốm Al2O3 trong suốt lên thêm ít nhất 25%.

  2. Nhúng công thức giải tích hai pha vào phần mềm mô phỏng vật liệu: Các nhóm phát triển phần mềm mô phỏng cơ tính và động lực học phân tử cần tích hợp mô hình giải tích mật độ khối lượng hai pha (pha MĐT mật độ 2,54 g/cm3 và pha MĐC mật độ 4,444 g/cm3) vào hệ thống thuật toán trong giai đoạn 2024-2025. Giải pháp này giúp rút ngắn hơn 35% thời gian tính toán lặp khi dự báo biến dạng thể tích của vật liệu chịu áp lực cao.

  3. Thiết kế cảm biến áp suất siêu cao dựa trên phổ tán xạ quang học: Các kỹ sư thuộc phòng thí nghiệm cảm biến và đo lường vật lý cần ứng dụng hiện tượng dịch chuyển tuyến tính của đỉnh thừa số cấu trúc từ 2,1 Å^-1 lên 2,4 Å^-1 để chế tạo đầu dò áp suất cực hạn cho dải 0 đến 30 GPa trong vòng 24 tháng tới, đảm bảo độ trễ đo đạc dưới 0,5 giây và sai số biên độ dưới 1,5%.

  4. Mở rộng mô hình nghiên cứu sang hệ ôxít pha tạp đa cấu tử: Các trường đại học và cơ sở nghiên cứu sau đại học cần triển khai mở rộng mô phỏng ĐLHPT cho các hệ Al2O3 pha tạp ôxít đất hiếm hoặc kim loại chuyển tiếp trong lộ trình nghiên cứu 3 năm, tập trung làm rõ cơ chế khuếch tán dị thường ở ngưỡng 15 đến 21 GPa nhằm phục vụ chế tạo vật liệu xúc tác chịu nhiệt độ trên 1500 độ C.

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

Công trình học thuật này cung cấp nguồn tư liệu chuyên sâu giá trị cao cho 4 nhóm đối tượng chính:

  1. Các nhà nghiên cứu và nghiên cứu sinh ngành Vật lý chất rắn, Hóa lý thuyết và Khoa học vật liệu: Tài liệu mang lại hệ phương pháp luận chuẩn xác về việc ứng dụng thế Born-Mayer, phép tính Ewald-Hansen và thuật toán phân tích đa diện phối trí, phục vụ đắc lực cho các đề tài mô phỏng cấu trúc chất lỏng và vật liệu vô định hình.

  2. Kỹ sư công nghệ tại các tập đoàn chế tạo gốm sứ kỹ thuật, kính quang học và linh kiện bán dẫn: Nắm vững cơ chế chuyển đổi giữa tứ diện AlO4 và bát diện AlO6 theo áp suất nén để làm chủ thông số công nghệ nhiệt luyện, gia công khuôn ép và tối ưu hóa độ cứng cơ học của sản phẩm.

  3. Giảng viên và học viên cao học tại các trường đại học khối khoa học tự nhiên và công nghệ: Khai thác toàn bộ số liệu phân bố xuyên tâm, phân bố góc liên kết và các mô hình không gian 3D để làm học liệu giảng dạy cho các học phần Vật lý chất rắn nâng cao, Vật lý tính toán và Khoa học vật liệu đại cương.

  4. Chuyên viên phân tích phổ nhiễu xạ tia X và tán xạ neutron tại các phòng thí nghiệm trọng điểm: Sử dụng các đường cong thừa số cấu trúc tán xạ tia X S(q) và hàm g(r) trong luận văn làm dữ liệu chuẩn để đối soát, bóc tách và tinh chỉnh phổ thực nghiệm thu được từ lò nung nhiệt độ cao trên 2000 K.

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

  1. Tại sao phương pháp mô phỏng động lực học phân tử lại vượt trội hơn thực nghiệm khi khảo sát Al2O3 ở nhiệt độ cao? Khảo sát thực nghiệm ở nhiệt độ nóng chảy 2327 K gặp trở ngại cực lớn do mẫu Al2O3 lỏng phản ứng hóa học mạnh với thành chén chứa gây nhiễm bẩn tạp chất. Phương pháp mô phỏng ĐLHPT với 3000 nguyên tử cho phép khảo sát trong môi trường tính toán tinh khiết tuyệt đối, tái lập chính xác cấu trúc vi mô tại dải nhiệt độ lên tới 5000 K và áp suất 30 GPa mà thực nghiệm không thể can thiệp trực tiếp.

  2. Đơn vị phối trí AlO5 đóng vai trò gì trong quá trình nén ép Al2O3 thủy tinh? Đơn vị AlO5 đóng vai trò là cấu trúc trung gian chuyển tiếp tối quan trọng giữa mạng tứ diện AlO4 và mạng bát diện AlO6. Khi mật độ mẫu tăng từ 2,85 g/cm3 lên 3,82 g/cm3, tỉ phần AlO5 tăng đột biến từ 19,9% lên chiếm ưu thế tuyệt đối ở mức 54,8%, giúp quá trình chuyển pha diễn ra liên tục, êm thuận mà không làm sụp đổ cấu trúc mạng ngẫu nhiên.

  3. Bản chất của mô hình hai pha trong Al2O3 thủy tinh và lỏng được hiểu như thế nào? Cấu trúc Al2O3 thực tế là sự tồn tại đồng thời của hai pha không gian: pha mật độ thấp gồm các cụm tứ diện AlO4 có mật độ riêng 2,54 g/cm3 và pha mật độ cao gồm các cụm AlO5, AlO6 có mật độ riêng 4,444 g/cm3. Mật độ tổng thể của vật liệu ở mọi điều kiện áp suất từ 0 đến 30 GPa hoàn toàn được tính toán chính xác thông qua tỷ lệ thể tích của hai pha này.

  4. Nguyên nhân nào gây ra hiện tượng khuếch tán dị thường ở áp suất 15 GPa và 21 GPa? Trong Al2O3 lỏng, hệ số khuếch tán của Al và O xuất hiện điểm dị thường cục bộ tại lần lượt 15 GPa và 21 GPa do sự mất ổn định của các khối đa diện AlO5. Dưới áp lực chuyển pha, các liên kết Al-O trong AlO5 bị kéo dãn cực đại và đứt gãy trước khi tái lập thành mạng bát diện AlO6 bền vững, tạo điều kiện cho các nguyên tử bứt phá và khuếch tán nhanh hơn.

  5. Độ tin cậy của các thông số cấu trúc trong luận văn khi so sánh với thực nghiệm đạt mức nào? Độ tin cậy của mô hình đạt mức rất cao khi hàm phân bố xuyên tâm g(r) ở mật độ 3,10 g/cm3 trùng khớp hoàn toàn với dữ liệu thực nghiệm nhiễu xạ tia X của Lamparter về vị trí đỉnh cực đại thứ nhất tại 1,80 Å, cực đại thứ hai tại 2,80 Å và hình thái suy giảm của đường cong tán xạ trong không gian thực.

Kết luận

  • Thiết lập thành công mô hình động lực học phân tử cho hệ Al2O3 vô định hình quy mô 3000 nguyên tử với thế tương tác Born-Mayer kết hợp gần đúng Ewald-Hansen đạt độ chính xác cao.
  • Phát hiện quy luật chuyển dịch cấu trúc mạnh mẽ từ tứ diện AlO4 (giảm từ 79,0% xuống 19,0%) sang phối trí trung gian AlO5 (tăng lên 54,8%) và bát diện AlO6 (tăng lên 26,5%) khi mật độ tăng từ 2,85 g/cm3 lên 3,82 g/cm3.
  • Xây dựng thành công các phương trình giải tích tường minh liên kết trực tiếp phân bố góc liên kết O-Al-O và Al-O-Al với tỷ phần các đơn vị cấu trúc đa diện.
  • Khẳng định bản chất tồn tại đồng thời của trạng thái hai pha: pha mật độ thấp (2,54 g/cm3) và pha mật độ cao (4,444 g/cm3), xác định ngưỡng chuyển pha của chất lỏng (7,5 đến 16,0 GPa) diễn ra sớm hơn trạng thái thủy tinh (12,0 đến 22,5 GPa).
  • Làm sáng tỏ cơ chế khuếch tán dị thường của nguyên tử Al tại 15 GPa và O tại 21 GPa trong pha lỏng gắn liền với sự phá vỡ liên kết trong đơn vị AlO5.

Công trình luận văn thạc sĩ của tác giả Nguyễn Hải Yến là một bước tiến quan trọng trong việc làm sáng tỏ vi cấu trúc của vật liệu ôxít nhôm vô định hình, đã được công bố chính thức trên Tạp chí Khoa học ĐHQGHN: Chuyên san Toán - Vật lý (Tập 38, Số 2, năm 2022). Kế hoạch nghiên cứu tiếp theo trong 12 đến 24 tháng tới sẽ mở rộng thuật toán sang các hệ gốm phức hợp chịu nhiệt siêu cao. Quý độc giả, các nhà nghiên cứu và doanh nghiệp quan tâm vui lòng kết nối với nhóm nghiên cứu tại Trường Đại học Sư phạm - Đại học Thái Nguyên để tiếp cận toàn văn dữ liệu mô phỏng và thúc đẩy các hợp tác chuyển giao công nghệ vật liệu tiên tiến.