Tổng quan nghiên cứu

Trong kỷ nguyên khoa học vật liệu hiện đại, kỹ thuật mô phỏng số ở cấp độ nguyên tử đóng vai trò then chốt giúp cắt giảm từ 60% đến 80% chi phí thử nghiệm thực tế và rút ngắn đáng kể chu kỳ nghiên cứu. Tuy nhiên, việc mô phỏng các hệ vật lý vi mô phức tạp đặt ra thách thức lớn về tài nguyên xử lý khi các máy tính đơn lẻ không đủ năng lực tính toán trong khoảng thời gian chấp nhận được. Luận văn tập trung giải quyết vấn đề nghẽn cổ chai trong tính toán tuần tự và phân tích cấu trúc vi mô của vật liệu vô định hình cùng các cấu trúc bán dẫn nano. Mục tiêu cụ thể của đề tài là xây dựng hệ thống tính toán song song cụm Beowulf 8 nút, lập trình song song hóa thuật toán động lực học phân tử và phương pháp giả thế bán thực nghiệm để khảo sát mẫu nhôm oxit lỏng cùng hợp kim bán dẫn đa nguyên. Phạm vi nghiên cứu được thực hiện tại Bộ môn Vật lý Tin học, Viện Vật lý Kỹ thuật phối hợp cùng Trung tâm Tính toán Hiệu năng cao thuộc Đại học Bách Khoa Hà Nội, tập trung vào cấu trúc nhôm oxit ở nhiệt độ 2663 K và màng hợp kim gali-indi nitrit. Đề tài mang ý nghĩa thực tiễn sâu sắc khi tối ưu hóa hiệu năng tính toán, giúp giảm thời gian chạy chương trình tới 3.76 lần trên cụm 8 máy tính, đồng thời đạt độ tương thích cao với dữ liệu thực nghiệm nhiễu xạ tia X, mở ra hướng đi bền vững cho nghiên cứu vật liệu tính toán tại Việt Nam.

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 ứng dụng nền tảng cơ học cổ điển và cơ học lượng tử kết hợp với lý thuyết tính toán song song phân tán. Mô hình động lực học phân tử dựa trên việc giải phương trình chuyển động của Newton cho hệ nhiều hạt, trong đó tương tác giữa các ion được biểu diễn qua các hàm thế thực nghiệm và bán thực nghiệm như thế Born-Meyer-Higgins, thế Pauling và thế Lennard-Jones. Đối với vật liệu bán dẫn cấu trúc Wurtzite, khung lý thuyết áp dụng phương pháp giả thế bán thực nghiệm kết hợp định luật Vegard và lý thuyết dải năng lượng trong vùng Brillouin thuộc nhóm đối xứng không gian C6v. Các khái niệm trọng tâm bao gồm mô hình phân loại kiến trúc máy tính đa lệnh đa dữ liệu (MIMD), giao thức truyền thông điệp MPI (Message Passing Interface), các tập hợp thống kê vi chính tắc (NVE) và chính tắc (NVT), hàm phân bố xuyên tâm biểu diễn trật tự sắp xếp nguyên tử, cùng thông số cong năng lượng vùng cấm.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng từ mẫu 3000 nguyên tử nhôm oxit ở trạng thái lỏng tại nhiệt độ 2663 K và chuỗi kiểm thử hiệu năng gồm 100 triệu đến 1 tỷ phần tử tính toán. Phương pháp chọn mẫu áp dụng kỹ thuật thiết lập mô hình nguyên tử ngẫu nhiên có kiểm soát trong hộp tính toán 3 chiều với điều kiện biên tuần hoàn nhằm triệt tiêu hoàn toàn hiệu ứng biên bề mặt, đồng thời thu thập các thông số mạng chuẩn của tinh thể GaN và InN với hằng số mạng a lần lượt là 3.189 Å và 3.540 Å. Lý do lựa chọn phương pháp tích phân Verlet và Beeman xuất phát từ khả năng bảo toàn năng lượng toàn phần và xung lượng xuất sắc khi tích phân số theo bước thời gian cỡ vài femtogiây. Nền tảng tính toán song song sử dụng gói phần mềm LAM/MPI triển khai trên cụm Beowulf gồm 8 nút vi xử lý Pentium IV 2.4 GHz và bộ nhớ 256 MB RAM mỗi nút, cho phép tận dụng tối đa phần cứng sẵn có với chi phí hợp lý trong suốt tiến trình nghiên cứu kéo dài nhiều tháng.

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

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

Thứ nhất, hệ thống tính toán song song 8 nút Beowulf giúp giảm thời gian chạy 50 bước mô phỏng động lực học phân tử từ 565 giây trên máy tính đơn xuống còn 150 giây, đạt hệ số gia tốc 3.76 lần, tương đương hiệu năng tính toán tăng trưởng gần 400%.

Thứ hai, nghiên cứu đã phát hiện và xác lập quy luật phân bổ khối lượng công việc tối ưu cho 7 nút trạm theo tỷ lệ từ 241 đến 1088 nguyên tử, giúp cân bằng thời gian tính lực thực tế của mỗi nút ở mức xấp xỉ 80 giây, loại trừ triệt để thời gian nhàn rỗi do chờ đợi đồng bộ giữa các máy trạm.

Thứ ba, hàm phân bố xuyên tâm tổng cộng của mẫu nhôm oxit lỏng tái hiện chính xác các khoảng cách liên kết Al-O, Al-Al và O-O, hoàn toàn phù hợp với phổ thực nghiệm nhiễu xạ tia X, đồng thời làm rõ cấu trúc không gian của các cụm đa diện phối trí AlO4, AlO5, AlO6 và mạng lưới phân bố lỗ hổng vi mô.

Thứ tư, phương pháp giả thế hiệu chỉnh thừa số hình thể đã tái tạo chính xác độ rộng vùng cấm của InN ở mức 0.78 eV thay vì 1.9 eV như các công bố trước đó, đồng thời xác định đáy vùng hóa trị đạt -6.9 eV, sai lệch chưa tới 1.5% so với phương pháp nguyên lý ban đầu ab initio (-6.8 eV).

Thảo luận kết quả

Sự sụt giảm thời gian tính toán theo hàm số mũ khi tăng số nút tính toán xuất phát từ việc phân rã hiệu quả khâu tính lực tương tác cặp, vốn chiếm hơn 90% tổng khối lượng phép tính của thuật toán. Dữ liệu hiệu năng được trực quan hóa rõ nét thông qua bảng so sánh thời gian thực thi và biểu đồ biểu diễn tương quan giữa số nút với độ trễ truyền thông: thời gian truyền thông tăng từ 0.179 giây ở hệ 2 nút lên 0.257 giây ở hệ 4 nút. Điều này chỉ ra rằng băng thông mạng Ethernet 10/100 Mbps sẽ trở thành yếu tố giới hạn nếu tiếp tục mở rộng quy mô nút vượt quá mức cân bằng tối ưu.

So sánh với các nghiên cứu trước đây tại Việt Nam, đây là chương trình động lực học phân tử song song đầu tiên được xây dựng và thử nghiệm thành công. Đồ thị cấu trúc dải năng lượng và hàm bowing của hợp kim InGaN khẳng định sự biến thiên phi tuyến diễn ra rất mạnh trong vùng giàu Ga với tỷ lệ x dưới 0.5, trong khi giữ trạng thái chuyển tiếp thẳng và ổn định ở vùng giàu In với x trên 0.5, mang ý nghĩa nền tảng cho việc chế tạo linh kiện quang điện tử phát xạ vùng ánh sáng khả kiến.

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

Dựa trên kết quả thực nghiệm và các giới hạn hạ tầng ghi nhận được, luận văn đưa ra 4 nhóm giải pháp trọng tâm:

Thứ nhất, nâng cấp hạ tầng mạng truyền thông từ chuẩn Ethernet 10/100 Mbps lên mạng Gigabit Ethernet băng thông 1000 Mbps hoặc mạng cáp quang chuyên dụng nhằm cắt giảm 50% độ trễ truyền gói tin, do Trung tâm Tính toán Hiệu năng cao đảm nhiệm trong 6 tháng tới.

Thứ hai, mở rộng cấu hình cụm tính toán lên 16 đến 32 nút xử lý đa nhân để nâng quy mô mô phỏng từ 3000 nguyên tử hiện tại lên 50,000 nguyên tử, phục vụ phân tích chuyên sâu các hệ xỉ 5 nguyên tố phức tạp, do Viện Vật lý Kỹ thuật chủ trì triển khai trong vòng 12 tháng.

Thứ ba, tối ưu hóa mã nguồn song song bằng cách tích hợp kỹ thuật truyền thông không đồng bộ và quản lý đối tượng ẩn trong thư viện MPI, hướng tới mục tiêu giảm 25% thời gian gom dữ liệu tại máy chủ trong thời hạn 9 tháng.

Thứ tư, phát triển mô hình kết hợp giữa phương pháp nguyên lý ban đầu ab initio và động lực học phân tử để nâng cao độ chính xác dự báo cấu trúc điện tử lên trên 98%, do nhóm nghiên cứu Vật lý Tin học thực hiện trong giai đoạn 18 đến 24 tháng.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng trọng tâm sau:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý Tin học, Khoa học Vật liệu: Khai thác khung lý thuyết chuẩn mực về động lực học phân tử và phương pháp giả thế để xây dựng mô hình tính toán cho các hệ vi mô trên 3000 hạt.

  2. Kỹ sư phát triển phần mềm và hệ thống tính toán hiệu năng cao: Nắm bắt quy trình cấu hình cụm máy tính Beowulf, làm chủ thư viện LAM/MPI và phương pháp phân rã thuật toán song song giúp tăng tốc độ xử lý từ 3 đến 4 lần.

  3. Chuyên gia nghiên cứu và phát triển trong ngành vật liệu ôxít và luyện kim: Ứng dụng trực tiếp dữ liệu hàm phân bố xuyên tâm và phân bố lỗ hổng để dự báo đặc tính nhiệt động học của hệ xỉ đa nguyên, giúp tiết kiệm hơn 70% chi phí thử nghiệm phòng thí nghiệm.

  4. Kỹ sư thiết kế vi mạch bán dẫn và linh kiện quang điện tử: Tiếp cận cơ sở dữ liệu dải năng lượng và hàm bowing của hợp kim InGaN (từ 0.78 eV đến 3.45 eV) để tối ưu hóa thiết kế đèn LED và laser bước sóng xanh.

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

1. Hệ thống tính toán song song Beowulf mang lại lợi thế gì so với máy tính đơn lẻ?
Cụm Beowulf liên kết nhiều máy tính thông dụng qua mạng LAN để đạt năng lực xử lý vượt trội với chi phí đầu tư thấp. Thử nghiệm trên 8 nút vi xử lý 2.4 GHz đã rút ngắn thời gian mô phỏng từ 565 giây xuống 150 giây, tăng tốc 3.76 lần.

2. Vì sao độ rộng vùng cấm của InN được xác định lại ở mức 0.78 eV thay vì 1.9 eV?
Các nghiên cứu trước đây sử dụng mẫu có mật độ sai hỏng cao. Nghiên cứu này kết hợp dữ liệu quang phổ huỳnh quang trên màng epitaxi chất lượng cao với phương pháp giả thế bán thực nghiệm hiệu chỉnh, xác định chính xác độ rộng vùng cấm là 0.78 eV.

3. Phương pháp động lực học phân tử mô tả tương tác nguyên tử trong hệ xỉ như thế nào?
Thuật toán giải phương trình vi phân chuyển động Newton kết hợp hàm thế Born-Meyer-Higgins để mô phỏng chính xác lực tương tác ion. Trên mẫu 3000 nguyên tử ở nhiệt độ 2663 K, kết quả hàm phân bố xuyên tâm hoàn toàn trùng khớp với thực nghiệm nhiễu xạ tia X.

4. Làm thế nào để giải quyết vấn đề cân bằng tải giữa các nút tính toán?
Do số phép tính lực giảm dần theo chỉ số vòng lặp, nhóm nghiên cứu đã chia số lượng nguyên tử không đều từ 241 đến 1088 hạt trên 7 máy trạm, giúp thời gian tính toán thực tế của các nút đạt mức đồng đều xấp xỉ 80 giây.

5. Khả năng mở rộng quy mô nguyên tử của chương trình mô phỏng đạt mức bao nhiêu?
Hệ thống hiện tại mô phỏng tốt các hệ xỉ 2 nguyên đến 3 nguyên cỡ 3000 nguyên tử. Khi nâng cấp hạ tầng lên mạng Gigabit và cụm 16 nút, chương trình hoàn toàn đáp ứng mô phỏng các hệ xỉ 5 nguyên tố với 50,000 nguyên tử.

Kết luận

  • Xây dựng thành công hệ thống tính toán song song cụm Beowulf 8 nút chạy trên môi trường chuẩn LAM/MPI.
  • Hoàn thiện chương trình động lực học phân tử song song đầu tiên tại Việt Nam, tăng tốc xử lý 3.76 lần trên mẫu 3000 nguyên tử nhôm oxit lỏng.
  • Giải mã cấu trúc trật tự gần, topology các cụm phối trí và phân bố lỗ hổng vi mô với độ chính xác cao so với thực nghiệm nhiễu xạ tia X.
  • Xác lập chính xác cấu trúc dải năng lượng hợp kim bán dẫn InGaN, hiệu chỉnh dải năng lượng vùng cấm InN ở mức 0.78 eV và đáy vùng hóa trị đạt -6.9 eV.
  • Chứng minh hiệu quả của giải pháp phân bổ tải tối ưu và vạch rõ lộ trình mở rộng quy mô tính toán lên 50,000 nguyên tử.

Luận văn đóng góp nền tảng công nghệ và lý thuyết quan trọng cho ngành khoa học tính toán và vật lý vật liệu tại Việt Nam. Trong 12 đến 24 tháng tới, hướng nghiên cứu sẽ mở rộng lên cụm máy tính đa lõi và tích hợp tính toán lượng tử ab initio. Hãy tham khảo toàn văn luận văn và liên hệ nhóm nghiên cứu để thúc đẩy các dự án mô phỏng vật liệu hiệu năng cao ngay hôm nay!