GIÁO TRÌNH: MÔ PHỎNG TRONG VẬT LÝ

Thông tin xuất bản và tác giả:

  • Chủ biên: GS. Võ Văn Hoàng (Trường Đại học Bách Khoa – ĐHQG TP.HCM)
  • Đồng tác giả: PGS. Huỳnh Kim Lâm (Trường Đại học Quốc tế – ĐHQG TP.HCM), TS. Nguyễn Trung Hải (Viện Công nghệ Illinois), TS. Nguyễn Hà Hùng Chương (Trường Đại học Tôn Đức Thắng)
  • Nhà xuất bản: Nhà xuất bản Đại học Quốc gia TP. Hồ Chí Minh (Năm xuất bản: 2016)
  • Thông tin đề tài: Sản phẩm của đề tài ĐHQG-HCM loại B (Mã số: B2014-20-01)

Tổng quan về giáo trình

Giáo trình Mô phỏng trong Vật lý được biên soạn làm tài liệu giảng dạy chính khóa cho sinh viên đại học và học viên cao học tại Trường Đại học Bách Khoa – ĐHQG TP.HCM, Trường Đại học Cần Thơ cùng nhiều cơ sở đào tạo chuyên ngành khoa học tự nhiên và kỹ thuật trên cả nước. Trong cơ cấu chương trình đào tạo, tài liệu đóng vai trò cầu nối giữa cơ sở lý thuyết vật lý, toán ứng dụng và khoa học máy tính, hình thành nên phương pháp nghiên cứu then chốt của lĩnh vực Khoa học Tính toán (Computational Science).

Mục tiêu học tập của giáo trình tập trung vào việc trang bị cho người học hệ thống phương pháp luận mô hình hóa và mô phỏng số. Sau khi hoàn thành nội dung tài liệu, người học có khả năng:

  • Hiểu rõ bản chất của quá trình mô phỏng trên máy tính như một phương thức thực nghiệm số độc lập.
  • Phân loại, lựa chọn và thiết lập các điều kiện biên, trường lực tương tác (thế tương tác cặp, thế nhiều hạt, thế tương tác Coulomb) cho các hệ vật liệu cụ thể.
  • Nắm vững nguyên lý và thuật toán của các phương pháp mô phỏng cổ điển (Monte-Carlo, Động lực học phân tử - MD) và các phương pháp lượng tử xuất phát từ các nguyên lý ban đầu (ab initio, DFT).
  • Ứng dụng thành thạo các phần mềm mô phỏng chuyên dụng như LAMMPS, GROMACS, Quantum ESPRESSO, VMD để khảo sát cấu trúc, động học và nhiệt động học của vật liệu khối, vật liệu nano, vật liệu 2 chiều và protein.

Giáo trình được cấu trúc thành 07 chương trọng tâm và 02 phần phụ lục. Cách tiếp cận của tài liệu đi từ các nguyên lý toán - lý cơ bản, xây dựng thuật toán, thiết kế chương trình mô phỏng đến việc phân tích các đối tượng vật lý cụ thể. Điểm đặc sắc của giáo trình là sự kết hợp chặt chẽ giữa lý thuyết vi mô và kỹ thuật thực thi tính toán, đồng thời cung cấp các dẫn chứng dữ liệu cụ thể từ các nghiên cứu thực nghiệm và mô phỏng đã được công bố.


Nội dung kiến thức cốt lõi

                                  KHOA HỌC TÍNH TOÁN
                     (Toán học + Khoa học máy tính + Chuyên ngành)
PHƯƠNG PHÁP CỔ ĐIỂN (Classical)                                  PHƯƠNG PHÁP LƯỢNG TỬ (Ab Initio)
                             ĐỐI TƯỢNG VÀ ỨNG DỤNG

Các chương và chủ đề chính

  1. Chương 1: Cơ sở của quá trình mô phỏng
    Trình bày bản chất của khoa học tính toán, phân loại hệ thống (gián đoạn và liên tục), so sánh giải pháp giải tích và mô phỏng. Trình bày chi tiết các kỹ thuật mô phỏng: trạng thái ban đầu, các loại điều kiện biên (tuần hoàn 3 hướng cho vật liệu khối; tuần hoàn 2 hướng cho màng mỏng; biên tự do, giả tự do, biên cứng cho hạt nano), các loại thế tương tác (Lennard-Jones, Morse, Pak-Doyama, Bohr-Mayer, EAM) và đơn vị ước lượng (đơn vị LJ, fs/ps). Giới thiệu đặc trưng cấu trúc chất lỏng, vô định hình và tinh thể qua hàm phân bố xuyên tâm $g(r)$, hệ số cấu trúc $S(K)$, số phối vị $Z$ và đa diện Voronoi. Thiết lập biểu đồ tương quan kích thước - thời gian - phương pháp mô phỏng.

  2. Chương 2: Một số phương pháp mô phỏng cổ điển dùng trong vật lý
    Tập trung vào phương pháp Metropolis Monte-Carlo (MC) cổ điển: quá trình sinh số giả ngẫu nhiên (RANDOM_SEED, RANDOM_NUMBER), bài toán bước ngẫu nhiên một chiều và đa chiều, điều kiện cân bằng chi tiết và chuỗi Markov. Trình bày phương pháp Động lực học phân tử cổ điển (MD) với các tập hợp thống kê vi chính tắc (NVE), đẳng nhiệt - đẳng áp (NPT), đẳng nhiệt - đẳng tích (NVT). Giới thiệu tổng quan cấu trúc phần mềm mã nguồn mở LAMMPS.

  3. Chương 3: Tổng quan về các phương pháp tính từ các nguyên lý ban đầu
    Xây dựng cơ sở cơ học lượng tử thông qua phương trình Schrödinger nhiều hạt, phương pháp gần đúng Hartree-Fock. Trình bày chi tiết Lý thuyết phiếm hàm mật độ (DFT) với hai định lý nền tảng Hohenberg-Kohn và hệ phương trình Kohn-Sham. Khảo sát phương pháp động lực học phân tử lượng tử: Born-Oppenheimer MD (BOMD) và Car-Parrinello MD (CPMD).

  4. Chương 4: Mô phỏng các hệ đơn giản đơn nguyên tử
    Khảo sát các hệ tương tác mẫu: Lennard-Jones, quả cầu cứng/mềm, đĩa cứng/mềm, thế Dzugutov, Lennard-Jones-Gauss, thế Jagla, thế Honeycomb. Phân tích các hiện tượng nhiệt động học và cấu trúc: quá trình làm lạnh nhanh tạo trạng thái vô định hình, hiện tượng tinh thể hóa và nóng chảy (đồng nhất ở vật liệu khối, không đồng nhất ở hệ có bề mặt tự do), cùng cơ chế chuyển pha lỏng - vô định hình.

  5. Chương 5: Mô phỏng vật liệu khối, vật liệu có kích thước nano và vật liệu 2 chiều
    Ứng dụng mô phỏng trên các đối tượng vật liệu cụ thể: vật liệu bán dẫn, aluminosilicate, kim loại và hợp kim khối. Khảo sát hiệu ứng kích thước và bề mặt trên màng mỏng, hạt nano. Phân tích cấu trúc và tính chất của vật liệu 2 chiều đơn lớp nguyên tử gồm Graphene, Silicene và các hệ 2D có thế tương tác đơn giản.

  6. Chương 6: Mô phỏng protein
    Trình bày cơ sở lý thuyết mô phỏng đại phân tử sinh học: các mô hình trường lực cơ học phân tử (MM), trường lực lai lượng tử/cơ học phân tử (QM/MM), trường lực Coarse-grained và mô hình dung môi không tường minh (Implicit Solvent). Giới thiệu các thuật toán nâng cao: cực tiểu hóa năng lượng, phương pháp trao đổi bản sao (Replica Exchange), các phương pháp tính năng lượng tự do (Thermodynamic Integration - TI, Free Energy Perturbation - FEP, Umbrella Sampling, Metadynamics).

  7. Chương 7: Giới thiệu về các phần mềm mô phỏng protein bằng phương pháp động lực học phân tử cổ điển và các phần mềm xử lý dữ liệu
    Phân loại các phần mềm MD phổ biến (GROMACS, NAMD, Materials Studio) và phần mềm xử lý trực quan dữ liệu (VMD, ORIGIN). Hướng dẫn quy trình thực hành chuẩn mô phỏng hệ protein trong nước bằng GROMACS qua các tệp cấu hình tham số .mdp (đưa ion vào hệ, cực tiểu hóa năng lượng, chạy cân bằng NVT, cân bằng NPT và chạy mô phỏng MD chính thức).

  • Phụ lục 1 & 2: Giới thiệu các phần mềm tính toán nguyên lý ban đầu (VASP, SIESTA, Quantum ESPRESSO) và bảng tra cứu thuật ngữ tiếng Anh chuyên ngành.

Kiến thức nền tảng được xây dựng

  • Lý thuyết vật lý: Cơ học cổ điển Newton, cơ học thống kê Boltzmann-Gibbs, cơ học lượng tử, nhiệt động lực học các quá trình chuyển pha, cấu trúc chất rắn và vật liệu vô định hình.
  • Mô hình hóa tương tác: Nguyên lý xây dựng trường lực cổ điển (thế cặp, thế ion Coulomb, thế nhúng nhiều hạt EAM) và tương tác phiếm hàm mật độ điện tử.
  • Toán học và thuật toán: Phương pháp số giải phương trình vi phân chuyển động, lý thuyết xác suất thống kê, chuỗi Markov, phép biến đổi Fourier chuyển đổi giữa $g(r)$ và $S(K)$, hình học không gian phân tích đa diện Voronoi tuân theo quy tắc Euler ($V - E + F = 2$).

Kỹ năng phát triển

  • Kỹ năng kỹ thuật: Viết và tối ưu hóa mã nguồn mô phỏng; thiết lập các điều kiện biên 1D, 2D, 3D; cấu hình tham số mô phỏng nhiệt độ, áp suất, bước thời gian tích phân ($\Delta t$).
  • Kỹ năng phân tích: Trích xuất và xử lý dữ liệu vi mô (tọa độ, vận tốc) thành các đại lượng vĩ mô; tính toán phân bố số phối vị, bán kính mặt cầu phối vị thứ nhất ($R_I$), khoảng cách liên kết trung bình ($r_1$).
  • Kỹ năng ứng dụng công cụ: Vận hành các phần mềm chuẩn quốc tế (LAMMPS, GROMACS, VMD, ORIGIN) và khai thác hệ thống tính toán hiệu năng cao chạy song song.

Phương pháp giảng dạy và học tập

Giáo trình triển khai phương pháp sư phạm tiếp cận theo tiến trình: Lý thuyết mô hình $\rightarrow$ Xây dựng thuật toán $\rightarrow$ Viết/chạy chương trình $\rightarrow$ Kiểm chứng tới hạn $\rightarrow$ Khảo sát ứng dụng thực tế.

  • Phương pháp kiểm chứng mô hình: Giáo trình hướng dẫn kỹ thuật kiểm tra tính đúng đắn của mã nguồn thông qua các trường hợp tới hạn. Ví dụ: khi viết chương trình tính số phối vị, người học phải thử nghiệm trên mô hình tinh thể lập phương tâm diện (fcc) để kiểm chứng số phối vị chính xác $Z = 12$. Nếu sai lệch, thuật toán cần được hiệu chỉnh trước khi áp dụng cho pha lỏng hoặc vô định hình.
  • Bài tập và thực hành quy trình: Giáo trình cung cấp quy trình từng bước mô phỏng protein trong nước bằng GROMACS thông qua hệ thống tệp .mdp:
    1. Đưa ion trung hòa điện tích vào hệ.
    2. Chạy cực tiểu hóa năng lượng (Energy Minimization).
    3. Cân bằng nhiệt độ hệ ở tập hợp NVT.
    4. Cân bằng áp suất hệ ở tập hợp NPT.
    5. Chạy mô phỏng động lực học chính (Production MD).
  • Đánh giá kết quả: Đánh giá độ tin cậy của mô hình bằng cách so sánh dữ liệu đầu ra với thực nghiệm (ví dụ: so sánh nhiệt độ nóng chảy của mô hình Ni dạng màng mỏng, phân bố hàm $g(r)$ và hệ số $S(K)$ của Fe vô định hình ở 300 K).
  • Hướng dẫn tự học: Giáo trình khuyến khích người học nắm vững cấu trúc chương trình con, cách đặt tên biến tường minh, phân tách mã nguồn và kết hợp công cụ trực quan hóa 3D (VMD) với công cụ phân tích đồ thị (ORIGIN).

Điểm nổi bật và cập nhật

  • Mức độ đổi mới nội dung: Phiên bản năm 2016 được biên soạn mới hơn 90% so với phiên bản năm 2004, mở rộng phạm vi từ các bài toán vật liệu cổ điển sang vật lý nano tính toán, vật lý sinh học tính toán và tính toán cấu trúc điện tử lượng tử.
  • Mở rộng các mô hình tương tác hiện đại: Bổ sung các thế tương tác tiên tiến như thế BKS (cho hệ $\text{SiO}_2$), thế nhiều hạt EAM/Finnis-Sinclair/Sutton-Chen (cho kim loại và hợp kim), và trường lực phản ứng (ReaxFF).
  • Cập nhật các hướng nghiên cứu mới:
    • Khảo sát hiện tượng chuyển pha và cấu trúc của vật liệu 2 chiều đơn lớp nguyên tử (Graphene, Silicene).
    • Ứng dụng các thuật toán tính năng lượng tự do tiên tiến trong sinh học phân tử: Metadynamics, Replica Exchange, Umbrella Sampling.
    • Phân tích cấu trúc vô định hình qua hiện tượng tách đỉnh thứ hai của hàm phân bố xuyên tâm $g(r)$ gắn liền với đơn vị cấu trúc đa giác 20 mặt (icosahedron).
  • Chuẩn hóa công cụ mô phỏng: Tích hợp hướng dẫn sử dụng các phần mềm mã nguồn mở chuẩn quốc tế được sử dụng rộng rãi trong nghiên cứu học thuật như LAMMPS, GROMACS, Quantum ESPRESSO.

Đối tượng sử dụng giáo trình

Nhóm đối tượng Chuyên ngành / Lĩnh vực Điều kiện tiên quyết (Prerequisites) Mục đích sử dụng
Sinh viên đại học (Năm 3, 4) Vật lý kỹ thuật, Khoa học vật liệu, Công nghệ Nano, Hóa học tính toán Vật lý đại cương, Cơ học lý thuyết, Nhiệt động lực học, Lập trình cơ bản (C/Fortran/Python) Học môn chính khóa Mô phỏng trong Vật lý, làm khóa luận tốt nghiệp
Học viên cao học & Nghiên cứu sinh Vật lý chất rắn, Vật lý lý thuyết và tính toán, Sinh học tính toán Cơ học lượng tử, Cơ học thống kê, Phương pháp tính toán số Thực hiện các đề tài nghiên cứu chuyên sâu, mô phỏng cấu trúc và tính chất vật liệu
Giảng viên & Nhà nghiên cứu Khoa học tính toán, Vật lý vật liệu, Khoa học vật liệu tính toán Kiến thức chuyên sâu về mô hình hóa và cấu trúc điện tử Thiết kế bài giảng, làm tài liệu tham khảo cho các đề tài nghiên cứu khoa học

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

1. Giáo trình này phù hợp với đối tượng nào?

Giáo trình được thiết kế cho sinh viên năm cuối đại học, học viên cao học, nghiên cứu sinh và các nhà nghiên cứu thuộc các chuyên ngành Vật lý, Khoa học Vật liệu, Công nghệ Nano, Hóa tính toán và Sinh học tính toán.

2. Cần chuẩn bị những kiến thức nền tảng nào trước khi học?

Người học cần nắm vững kiến thức về Cơ học cổ điển, Cơ học lượng tử, Cơ học thống kê, Nhiệt động lực học, Toán ứng dụng (giải tích, phương pháp số) và có kỹ năng lập trình căn bản (sử dụng ngôn ngữ C, Fortran hoặc Python).

3. Điểm khác biệt của giáo trình này so với các tài liệu khác là gì?

Tài liệu bao quát liền mạch từ các phương pháp cổ điển (Metropolis Monte-Carlo, Molecular Dynamics) đến các phương pháp lượng tử nguyên lý ban đầu (ab initio, DFT), kết hợp phân tích đối tượng vật liệu đa dạng từ hệ đơn nguyên tử, vật liệu khối, màng mỏng, hạt nano, vật liệu 2D (graphene, silicene) cho đến các đại phân tử sinh học (protein). Đồng thời cung cấp quy trình cấu hình chi tiết cho các phần mềm chuẩn như GROMACS và LAMMPS.

4. Làm sao để tự học và thực hành giáo trình hiệu quả?

Người học nên kết hợp đọc lý thuyết với việc tự viết các đoạn mã nguồn thuật toán đơn giản (như bước ngẫu nhiên, thuật toán Metropolis), kiểm tra mã nguồn trên các hệ mẫu tới hạn (như mạng fcc với $Z = 12$), sau đó chuyển sang thực hành chạy mô phỏng các hệ thực tế trên các phần mềm chuyên dụng (LAMMPS, GROMACS) và hiển thị kết quả bằng VMD, ORIGIN.

5. Giáo trình có tài liệu tham khảo và phụ lục bổ trợ nào?

Mỗi chương trong giáo trình đều đi kèm danh mục tài liệu tham khảo quốc tế uy tín (như sách của Allen & Tildesley, Frenkel & Smit, D. C. Rapaport, C. Kittel, bài giảng MIT). Phần phụ lục cung cấp danh mục phần mềm mô phỏng ab initio (VASP, SIESTA, Quantum ESPRESSO) và bảng thuật ngữ tiếng Anh chuyên ngành.


Kết luận

Giáo trình Mô phỏng trong Vật lý (GS. Võ Văn Hoàng chủ biên, NXB ĐHQG-HCM, 2016) là tài liệu học thuật toàn diện về lĩnh vực Khoa học Tính toán ứng dụng trong vật lý chất rắn, vật lý vật liệu và vật lý sinh học.

Lộ trình tiếp cận kiến thức được tối ưu hóa theo các giai đoạn:

  1. Giai đoạn 1 (Nền tảng): Nắm vững kỹ thuật mô hình hóa, điều kiện biên, thế tương tác và lý thuyết cấu trúc ($g(r)$, $S(K)$, Voronoi) tại Chương 1.
  2. Giai đoạn 2 (Phương pháp luận): Nghiên cứu các thuật toán cốt lõi trong mô phỏng Monte-Carlo, Động lực học phân tử và tính toán cấu trúc điện tử DFT tại Chương 2 và Chương 3.
  3. Giai đoạn 3 (Ứng dụng chuyên sâu): Triển khai khảo sát trên các hệ vật liệu cụ thể (khối, màng mỏng, nano, 2D) tại Chương 4, 5 hoặc hệ sinh học (protein) và thực hành phần mềm GROMACS, LAMMPS tại Chương 6, 7.

Tài liệu cung cấp cơ sở lý thuyết vững chắc và công cụ thực hành chuẩn xác cho quá trình đào tạo và nghiên cứu khoa học chuyên ngành Vật lý tính toán.