Tổng quan nghiên cứu

Trên quy mô toàn cầu, hơn 2 triệu ca phẫu thuật ghép xương được thực hiện mỗi năm, đồng thời ngành năng lượng nguyên tử đang phải quản lý lượng chất thải phóng xạ mức cao chiếm tới 95% tổng hoạt độ phóng xạ phát tán. Thực tiễn này đặt ra nhu cầu cấp thiết về việc phát triển các hệ vật liệu thủy tinh sinh học và thủy tinh lưu trữ chất thải có độ bền hóa học cao, khả năng tương thích sinh học tối ưu và tính ổn định cấu trúc vượt trội. Hệ vật liệu ba thành phần CaO – P2O5 – SiO2 là một trong những nền tảng quan trọng nhất đáp ứng cả hai lĩnh vực y sinh và kỹ thuật môi trường.

Vấn đề nghiên cứu cốt lõi tập trung vào việc làm rõ cơ chế biến đổi vi cấu trúc của mạng lưới thủy tinh khi thay đổi tỷ lệ thành phần giữa chất tạo mạng P2O5 và chất biến tính mạng CaO. Mục tiêu cụ thể của luận văn là mô phỏng động lực học phân tử hệ vật liệu CaO – P2O5 – SiO2 ở nhiệt độ 3000 K với dải nồng độ P2O5 biến thiên từ 5 mol% đến 40 mol%, tương ứng tỷ lệ CaO giảm từ 35 mol% xuống 0 mol% trên nền 60 mol% SiO2 cố định.

Phạm vi nghiên cứu được thực hiện tại Bộ môn Vật lý Tính toán, Viện Vật lý Kỹ thuật, Đại học Bách khoa Hà Nội, hoàn thành vào tháng 10 năm 2023. Luận văn đã tính toán chi tiết 7 mô hình mô phỏng với quy mô từ 5270 đến 5520 nguyên tử trong hộp lập phương kích thước 53 Å. Kết quả nghiên cứu mang ý nghĩa học thuật và thực tiễn sâu sắc, cung cấp bộ thông số định lượng về độ dài liên kết, góc liên kết, tỷ lệ oxy cầu nối và sự phân bố cụm tứ diện, tạo cơ sở dữ liệu quan trọng giúp kiểm soát tốc độ hòa tan sinh học và nâng cao khả năng cố định ion phóng xạ trong vật liệu mới.

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 nền tảng lý thuyết mạng lưới vô định hình của Zachariasen và lý thuyết trường phối trí vô cơ, phân loại các thành phần thủy tinh thành ba nhóm: chất tạo mạng (SiO2, P2O5), chất biến tính mạng (CaO) và các đơn vị cấu trúc trung gian. Cấu trúc của hệ thủy tinh silicat - phosphat được mô tả thông qua hai mức độ trật tự không gian: trật tự gần (SRO) thể hiện qua các khối đa diện phối trí và trật tự trung bình (IRO) thể hiện qua liên kết giữa các đơn vị cấu trúc cơ bản.

Năm khái niệm then chốt được áp dụng xuyên suốt mô hình nghiên cứu bao gồm:

  1. Hàm phân bố xuyên tâm (RDF - g(r)): Đại lượng thống kê xác định xác suất tìm thấy hạt ở khoảng cách r so với hạt mốc, từ đó suy ra khoảng cách liên kết trung bình.
  2. Số phối trí (CN): Số lượng nguyên tử oxy bao quanh một cation trung tâm trong bán kính cắt xác định.
  3. Oxy cầu nối (BO): Nguyên tử oxy liên kết đồng thời với ít nhất hai cation tạo mạng Si hoặc P để duy trì tính liên tục của khung mạng.
  4. Oxy không liên kết cầu (NBO): Nguyên tử oxy chỉ liên kết với một cation tạo mạng và mang điện tích âm cục bộ được bù trừ bởi cation Ca2+.
  5. Mức độ trùng hợp mạng lưới: Tỷ lệ tương đối giữa BO và NBO phản ánh mật độ liên kết không gian ba chiều của hệ vật liệu.

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ử (MDS) cổ điển kết hợp điều kiện biên tuần hoàn ba chiều nhằm loại bỏ hoàn toàn hiệu ứng bề mặt trên mẫu kích thước nano. Cỡ mẫu nghiên cứu gồm 7 hệ vi mô với số lượng từ 5270 đến 5520 nguyên tử, được tạo cấu hình ban đầu bằng cách phân bố ngẫu nhiên các hạt trong hộp mô phỏng lập phương kích thước 53 Å. Lý do lựa chọn phương pháp này là vì MDS cho phép quan sát chi tiết sự dịch chuyển vị trí của từng nguyên tử ở cấp độ femtogiây, điều mà các phương pháp thực nghiệm nhiễu xạ tia X hay phổ hồng ngoại rất khó phân tách rõ ràng ở pha vô định hình nhiệt độ cao.

Quy trình mô phỏng sử dụng thế tương tác cặp Born-Mayer-Huggins (BMH) bao gồm ba thành phần: tương tác Coulomb tầm xa, lực đẩy Born tầm ngắn và tương tác hút Van der Waals. Phương trình chuyển động của các hạt được tích phân số theo thuật toán Verlet với bước thời gian 1 femtogiây trong tập hợp đẳng nhiệt - đẳng áp (NPT).

Timeline quy trình nhiệt động học được chuẩn hóa qua các bước: mẫu được gia nhiệt lên 5000 K trong 100.000 bước thời gian để phá vỡ hoàn toàn cấu hình ban đầu, sau đó làm nguội dần qua các mức nhiệt độ 4500 K, 4000 K, 3500 K với thời gian hồi phục 100.000 bước tại mỗi mức, và cuối cùng đưa về trạng thái cân bằng nhiệt động ở 3000 K trong 1.000.000 bước thời gian để trích xuất dữ liệu vi cấu trúc phục vụ phân tích thống kê.

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

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

Phân tích số liệu từ 7 mô hình mô phỏng cho thấy các quy luật biến đổi cấu trúc rõ rệt:

Thứ nhất, tính ổn định của cấu trúc trật tự gần (SRO): Khoảng cách liên kết trung bình giữa các cặp nguyên tử hầu như không đổi khi tăng nồng độ P2O5 từ 5 mol% lên 40 mol%. Cụ thể, khoảng cách Si-O duy trì ổn định ở mức 1.60 đến 1.62 Å, P-O duy trì ở mức 1.52 đến 1.56 Å và Ca-O dao động quanh 2.28 đến 2.32 Å. Phối trí của Si và P chủ yếu tồn tại ở dạng tứ diện SiO4 (chiếm 74.47% đến 80.00%) và PO4 (chiếm 66.00% đến 100.00%), với góc liên kết O-Si-O đạt 108 đến 112 độ và góc O-P-O đạt 110 đến 112 độ.

Thứ hai, sự gia tăng mạnh mẽ của mức độ trùng hợp mạng lưới: Khi nồng độ P2O5 tăng từ 5 mol% lên 40 mol% (đồng thời CaO giảm từ 35 mol% về 0 mol%), tỷ lệ oxy cầu nối (BO) tăng vọt từ 55.75% lên 94.75% (tăng 39.00 điểm phần trăm). Ngược lại, tỷ lệ oxy phi cầu nối (NBO) sụt giảm mạnh từ 44.25% xuống chỉ còn 5.25%, trong khi tỷ lệ oxy tự do (FO) gần như triệt tiêu hoàn toàn (xấp xỉ 0%).

Thứ ba, sự chuyển dịch các dạng liên kết oxy: Liên kết Si-O-Si giảm từ mức cực đại 87.50% ở mẫu 5 mol% P2O5 xuống khoảng 50.00% ở mẫu 40 mol% P2O5. Liên kết lai P-O-Si tăng nhanh từ 12.46% đạt đỉnh 51.31% tại nồng độ 20 mol% P2O5, sau đó giảm về 32.00%. Đồng thời, liên kết đồng loại P-O-P tăng đều đặn từ 0% lên 18.20%.

Thứ tư, sự phân mảnh và tái cấu trúc cụm tứ diện: Tại 5 mol% P2O5, toàn bộ mạng lưới SiO4 liên kết thành một cụm siêu lớn duy nhất chứa 3947 nguyên tử. Khi nồng độ P2O5 tăng lên 40 mol%, cụm SiO4 lớn nhất bị phân mảnh chỉ còn 123 nguyên tử, số lượng các cụm nhỏ 5 nguyên tử tăng từ 9 cụm lên 124 cụm. Ngược lại, các cụm PO4 phát triển từ các đơn vị cô lập 5 nguyên tử (186 cụm tại 5 mol%) thành các đại phân tử chứa hơn 20 nguyên tử ở dải nồng độ cao.

Thảo luận kết quả

Các phát hiện trên có thể được biểu diễn trực quan thông qua đồ thị đường biểu diễn hàm phân bố xuyên tâm RDF với các đỉnh sắc nét của Si-O và P-O, kết hợp cùng biểu đồ cột chồng phân bố tỷ lệ BO/NBO và biểu đồ phân bố kích thước cụm tứ diện.

Về mặt bản chất vật lý, việc bổ sung P2O5 đóng vai trò kép: vừa là chất tạo mạng bổ sung mạng lưới liên kết phosphat, vừa làm giảm tỷ lệ cation biến tính Ca2+. Khi nồng độ CaO giảm, số lượng liên kết Si-O-Ca bị phá vỡ, các liên kết Si-O tự do tái kết hợp với P để tạo thành mạng liên kết hỗn hợp Si-O-P bền vững. Sự tồn tại đỉnh liên kết lai Si-O-P tại điểm 20 mol% P2O5 chứng minh đây là ngưỡng nồng độ chuyển pha cấu trúc tối ưu giữa hai phân mạng silicat và phosphat.

So sánh với các nghiên cứu trước đây của Antonio Tilocca về thủy tinh sinh học 45S5 ở 3000 K và nghiên cứu của Jiang Diao trên hệ CaO-P2O5-SiO2 ở 1673 K, kết quả của luận văn hoàn toàn tương thích về khoảng cách liên kết Si-O (1.62 Å) và P-O (1.56 Å). Tuy nhiên, luận văn đã làm sáng tỏ thêm cơ chế phân bố dị thể cục bộ: hệ vật liệu xuất hiện sự phân tách vi mô thành các vùng giàu phosphat, vùng giàu silicat và các vùng tập trung cation Ca2+ tại ranh giới giữa hai phân mạng. Cation Ca2+ có số phối trí trung bình từ 5.27 đến 5.36 (chủ yếu là phối trí 5 đỉnh CaO5 chiếm hơn 35%), ưu tiên định vị gần các liên kết NBO và liên kết cầu nối P-O-Si, đóng vai trò then chốt trong việc kiểm soát độ linh động của các ion mang điện.

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

Dựa trên các kết luận định lượng về vi cấu trúc, bốn giải pháp kỹ thuật cụ thể được đề xuất nhằm ứng dụng hệ vật liệu CaO – P2O5 – SiO2 vào sản xuất thực tế:

  1. Tối ưu hóa tỷ lệ thành phần cho vật liệu cấy ghép y sinh: Thiết lập dải nồng độ P2O5 ở mức 15 mol% đến 20 mol% kết hợp 20 mol% đến 25 mol% CaO trên nền 60 mol% SiO2. Mục tiêu đạt tỷ lệ liên kết lai Si-O-P cực đại trên 50% và tỷ lệ NBO duy trì trong khoảng 25% đến 35%. Chỉ số này đảm bảo cân bằng tối ưu giữa độ bền cơ học và tốc độ giải phóng ion Ca2+, PO4 3- để kích thích hình thành khoáng xương Hydroxyapatite trong thời gian 12 tháng. Chủ thể thực hiện: Kỹ sư nghiên cứu phát triển vật liệu y sinh tại các viện nghiên cứu và doanh nghiệp sản xuất thiết bị y tế.

  2. Cố định chất thải phóng xạ mức cao trong công nghệ hạt nhân: Áp dụng tỷ lệ P2O5 từ 30 mol% đến 40 mol% nhằm nâng tỷ lệ oxy cầu nối (BO) lên trên 90% và hạn chế tối đa các liên kết NBO yếu dưới 10%. Cấu trúc mạng lưới liên tục mật độ cao này cho phép khóa chặt các cation kim loại phóng xạ và actinide bên trong khung -O-P-O-Si-O-, đạt mục tiêu nâng độ bền bức xạ và chống rửa trôi lên trên 95% với chu kỳ lưu trữ an toàn từ 20 đến 50 năm. Chủ thể thực hiện: Các trung tâm xử lý chất thải hạt nhân và viện công nghệ môi trường.

  3. Mở rộng quy mô mô phỏng và ứng dụng trí tuệ nhân tạo: Xây dựng các mô hình động lực học phân tử quy mô lớn từ 50.000 đến 100.000 nguyên tử kết hợp thế tương tác máy học (Machine Learning Potentials) trong giai đoạn 2024–2026. Hướng đi này nhằm mục tiêu cắt giảm 40% chi phí tính toán đồng thời mô phỏng chính xác hành vi cơ lý ở dải nhiệt độ phòng 300 K. Chủ thể thực hiện: Nhóm nghiên cứu vật lý tính toán và khoa học vật liệu tại các trường đại học.

  4. Triển khai thử nghiệm thực nghiệm kiểm chứng nhiệt động: Tiến hành nấu luyện mẫu thủy tinh thực tế theo 7 dải tỷ lệ nồng độ tương ứng và kiểm tra bằng phương pháp phổ cộng hưởng từ hạt nhân (NMR) cùng phổ Raman trong vòng 18 tháng. Mục tiêu kiểm chứng độ sai lệch vi cấu trúc thực tế so với mô hình mô phỏng dưới 5%. Chủ thể thực hiện: Các phòng thí nghiệm trọng điểm quốc gia về vật liệu tiên tiến.

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

Luận văn là tài liệu tham khảo chuyên sâu, cung cấp giá trị học thuật và ứng dụng thiết thực cho bốn nhóm đối tượng chính:

Thứ nhất, Nhà nghiên cứu và phát triển vật liệu y sinh: Nắm bắt chính xác cơ chế liên kết của các đơn vị PO4 và SiO4 để thiết kế các thế hệ thủy tinh hoạt tính sinh học mới, tối ưu hóa khả năng tương thích mô xương và nâng cao hiệu quả điều trị trong phẫu thuật chỉnh hình nha khoa.

Thứ hai, Kỹ sư công nghệ hạt nhân và quản lý môi trường: Sử dụng các dữ liệu định lượng về độ trùng hợp mạng lưới và cơ chế tích hợp cation kim loại kiềm thổ Ca2+ vào khung phosphat để tính toán công thức thủy tinh hóa tối ưu cho việc lưu trữ chất thải nguy hại và chất thải phóng xạ mức cao an toàn dài hạn.

Thứ ba, Chuyên gia kỹ thuật luyện kim và sản xuất xỉ thép: Khai thác quy luật biến đổi pha giữa 2CaO.SiO2 và 3CaO.P2O5 để tối ưu hóa quá trình khử photpho trong công nghệ luyện thép chất lượng cao, giúp giảm chi phí phụ gia và kiểm soát độ nhớt của xỉ lỏng ở nhiệt độ cao.

Thứ tư, Giảng viên, học viên cao học và sinh viên ngành Vật lý Kỹ thuật, Khoa học Vật liệu Tính toán: Sử dụng luận văn như một tài liệu quy chuẩn về phương pháp mô phỏng động lực học phân tử (MDS), kỹ thuật thiết lập thế tương tác BMH, thuật toán tích phân Verlet và phương pháp phân tích thống kê cấu trúc vô định hình.

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

  1. Tại sao nghiên cứu lại lựa chọn mô phỏng hệ vật liệu ở nhiệt độ cao 3000 K? Nhiệt độ 3000 K đưa hệ vật liệu về trạng thái nóng chảy hoàn toàn, giúp các nguyên tử đạt độ linh động cao nhằm loại bỏ cấu hình ban đầu ngẫu nhiên trong 1.000.000 bước mô phỏng. Trạng thái lỏng cân bằng ở 3000 K phản ánh chính xác các đặc trưng nhiệt động học và là cơ sở nền tảng để nghiên cứu quá trình đông đặc tạo thủy tinh vô định hình.

  2. Tỷ lệ P2O5 ảnh hưởng trực tiếp đến mức độ trùng hợp mạng lưới như thế nào? Khi tăng nồng độ P2O5 từ 5 mol% lên 40 mol%, tỷ lệ oxy cầu nối (BO) tăng tuyến tính từ 55.75% lên 94.75%, trong khi oxy không liên kết cầu (NBO) giảm từ 44.25% xuống 5.25%. Sự gia tăng này chứng minh P2O5 đóng vai trò chất tạo mạng mạnh mẽ, thúc đẩy sự hình thành khung liên kết ba chiều đồng nhất và bền vững.

  3. Vị trí và vai trò vi mô của ion Ca2+ trong cấu trúc mạng là gì? Ion Ca2+ đóng vai trò chất biến tính mạng với số phối trí trung bình từ 5.27 đến 5.36 (chủ yếu là đơn vị phối trí 5 oxy CaO5 chiếm hơn 35%). Ion Ca2+ ưu tiên định vị tại ranh giới chuyển tiếp giữa phân mạng POx và SiOx, liên kết trực tiếp với các nguyên tử NBO mang điện tích âm để ổn định hóa cấu trúc mạng lưới.

  4. Sự biến đổi kích thước cụm tứ diện SiO4 mang ý nghĩa cấu trúc gì? Tại 5 mol% P2O5, toàn bộ tứ diện SiO4 liên kết thành một đại cụm duy nhất gồm 3947 nguyên tử. Khi nồng độ P2O5 đạt 40 mol%, đại cụm này bị phân mảnh thành nhiều cụm nhỏ với kích thước cực đại chỉ còn 123 nguyên tử. Hiện tượng này minh chứng cho sự thâm nhập của mạng phosphat chia cắt mạng silicat thành các cấu trúc vi dị thể.

  5. Kết quả mô phỏng cấu trúc có thể áp dụng thực tế vào sản xuất vật liệu y sinh ra sao? Dữ liệu tỷ lệ liên kết Si-O-P đạt đỉnh 51.31% tại 20 mol% P2O5 cho phép xác định chính xác tỷ lệ phối trộn tối ưu để tạo ra thủy tinh sinh học có độ bền cơ học cao nhưng vẫn giải phóng ion Ca2+ và phosphat đều đặn, kích thích quá trình tái tạo xương nhanh hơn so với công thức thủy tinh sinh học 45S5 truyền thống.

Kết luận

Luận văn Thạc sĩ của tác giả Nguyễn Mai Anh đã giải quyết trọn vẹn bài toán mô phỏng và phân tích vi cấu trúc hệ vật liệu ba thành phần CaO – P2O5 – SiO2. Năm đóng góp khoa học cốt lõi bao gồm:

  • Xác lập bộ thông số cấu trúc trật tự gần (SRO) ổn định với khoảng cách liên kết Si-O (1.62 Å), P-O (1.56 Å) và Ca-O (2.30 Å) không phụ thuộc vào nồng độ thành phần.
  • Định lượng hóa sự gia tăng của mức độ trùng hợp mạng lưới, trong đó tỷ lệ oxy cầu nối (BO) tăng từ 55.75% lên 94.75% khi tăng P2O5 từ 5 mol% lên 40 mol%.
  • Phát hiện quy luật phân mảnh cụm tứ diện SiO4 từ 3947 nguyên tử xuống 123 nguyên tử, đồng thời làm rõ sự phát triển của cụm PO4 ở nồng độ phosphat cao.
  • Làm sáng tỏ cơ chế phân bố dị thể vi mô và vai trò biến tính mạng của cation Ca2+ với số phối trí trung bình xấp xỉ 5.30 định vị tại ranh giới hai phân mạng.
  • Cung cấp cơ sở dữ liệu tính toán chuẩn xác từ 7 mô hình với hơn 5200 nguyên tử, tạo tiền đề cho việc thiết kế vật liệu thủy tinh sinh học và gốm thủy tinh lưu trữ phóng xạ thế hệ mới.

Kế hoạch phát triển tiếp theo trong giai đoạn 2024–2026 sẽ tập trung vào việc mở rộng mô phỏng quá trình làm nguội từ 3000 K về nhiệt độ phòng 300 K và kiểm chứng thực nghiệm bằng các phương pháp phổ học tiên tiến. Hãy liên hệ với nhóm nghiên cứu tại Đại học Bách khoa Hà Nội hoặc tham khảo toàn văn luận văn để hợp tác phát triển các ứng dụng vật liệu chức năng cao cấp!