ÔN THI CUỐI KỲ TÍNH TOÁN VÀ MÔ PHỎNG VẬT LIỆU
Câu 1 Mô phỏng động lực học phân tử dùng để làm gì ? .......................................................................................................................... 1
Câu 2 Cơ sở lý thuyết của mô phỏng động lực học phân tử ? .................................................................................................................... 2
Câu 3 Nguyên lý nào trong vật lý thống kê là cơ sở của ý tưởng phát triển phương pháp mô phỏng động lực học phân tử ? Phát biểu
nguyên lý và trình bày những hiểu biết của em về mối tương quan giữa hệ vi mô và hệ vi mô thông qua nguyên lý đó . ....................... 3
Câu 4 Trong mô phỏng MD trường lực có quan trọng không , tại sao? Trường lực cung cấp nhữung thông tin gì , Hãy nêu các thành
phần cơ bản của một trường lực mô tả tương tác giữa các nguyên tử trong một hệ sinh học. ................................................................... 3
Câu 5 Liệt kê các điều kiện biên trong mô phỏng mà em biết và nêu rõ điều kiện biên nào sẽ sử dụng cho mục đích mô phỏng nào ..... 4
Câu 6 Nêu các chế độ mô phỏng và các thống số, làm sao để duy trì toàn bộ áp suất , nhiệt độ trong quá trình ? ................................... 5
Câu 7 Phương pháp nào được sử dụng để đưa hệ về trước khi mô phỏng ................................................................................................. 6
Câu 8 Phương pháp nào được sử dụng để cấu trúc thay đổi, chức năng thay đổi. ..................................................................................... 7
Câu Tính .................................................................................................................................................................................................... 8
Câu 10. Hãy nêu quá trình thực hiện mô phỏng MD cho hệ nhiều hạt, cho mô phỏng protein, Graphene ................................................ 9
Câu 11. Cách phân tích tính chất vật liệu xem đã hoàn toàn chuyển pha chưa, nêu ý nghĩa của các hàm tính. ....................................... 11
Đề cuối kỳ 07/01/2025
Câu 1 (4đ) Cơ sở lý thuyết của mô phỏng động lực học phân tử ? nguyên lý nào trong vật lý thống kê là cơ sở
của ý tưởng phát triển phương pháp mô phỏng động lực học phân tử ? Phát biểu nguyên lý và trình bày những
hiểu biết của em về mối tương quan giữa hệ vi mô và hệ vi mô thông qua nguyên lý đó .
Câu 2 (2đ) Liệt kê các điều kiện biên trong mô phỏng mà em biết và nêu rõ điều kiện biên nào sẽ sử dụng cho
mục đích mô phỏng nào
Câu 3(4đ) Giả sử em cần mô phỏng quá trình chuyến pha từ rắn sang lỏng của màng mỏng graphene đơn lớp,
em cần thực hiện qua những bước nào . Sau mô phỏng các tính chất nào của vật liệu cần được phân tích để
khẳng định hệ đã chuyển pha hoàn toàn.
Câu 1 Mô phỏng động lực học phân tử dùng để làm gì ?
Tiên đoán và giải mã cơ chế của các hệ phức tạp (Bài toán đa hạt)
Giải quyết "bài toán đa hạt": Dự đoán tính chất của cả hệ thống lớn dựa trên sự tương tác giữa hàng tỷ nguyên
tử, điều mà việc cộng gộp đơn giản không làm được
Cầu nối Vi mô - Vĩ mô: Giải thích nguồn gốc các tính chất vĩ mô (nhiệt độ, áp suất, độ nhớt...) thông qua
chuyển động vi mô của nguyên tử.
Hiểu rõ bản chất: Không chỉ cho biết kết quả mà còn trực quan hóa quá trình, giúp trả lời câu hỏi "tại sao" vật
liệu có tính chất đó.
Kiểm chứng các lý thuyết và mô hình (Thí nghiệm trên máy tính)
Đóng vai trò cầu nối, kiểm tra tính đúng đắn của các lý thuyết thuần túy trong môi trường ảo được kiểm soát
chặt chẽ từ đó xác định xem các giả thuyết đơn giản hóa của mô hình lý thuyết có phù hợp với thực tế hay
không. → Hiệu chỉnh mô hình: Đối chiếu kết quả mô phỏng với thực nghiệm để tinh chỉnh các tham số tương
tác cho chính xác hơn.
Hỗ trợ và thay thế cho các thí nghiệm khó khăn
MD là công cụ đắc lực để vượt qua các giới hạn của thực nghiệm truyền thống:
Điều kiện bất khả thi (Impossible): Mô phỏng các môi trường con người không thể tiếp cận (lõi Trái Đất, bề mặt
Mặt Trời, lò phản ứng hạt nhân...).
Môi trường nguy hiểm (Dangerous): Nghiên cứu an toàn các vật liệu phóng xạ, chất độc hóa học, virus, hoặc
các quá trình cháy nổ.
Vấn đề chi phí và thời gian (Expensive): Tiết kiệm chi phí: Giảm thiểu việc mua hóa chất/thiết bị đắt tiền nhờ
quá trình "sàng lọc" vật liệu trên máy tính trước khi làm thực nghiệm.
Điều kiện khắc nghiệt: Tránh tốn kém năng lượng để duy trì áp suất siêu cao hoặc nhiệt độ siêu thấp.
Giới hạn quan sát: Thang thời gian: Quan sát các sự kiện cực nhanh (femtogiây) như quá trình đứt gãy liên kết
mà máy móc không bắt kịp. Thang không gian: Nhìn thấy chi tiết cấu trúc hạt nhân, quark hoặc tương tác
nguyên tử mà kính hiển vi không soi tới được.
Câu 2 Cơ sở lý thuyết của mô phỏng động lực học phân tử ?
Cơ sở lý thuyết của Mô phỏng Động lực học Phân tử (MD) được xây dựng dựa trên sự kết hợp giữa Cơ học Cổ
điển, Toán học tính toán và Cơ học Thống kê
Nguyên lý cốt lõi cơ học Cổ điển (Newton) :Mô phỏng MD dựa trên giả thuyết rằng các nguyên tử và phân tử
tuân theo các định luật của cơ học cổ điển, bỏ qua các hiệu ứng lượng tử đối với chuyển động của hạt nhân.
Trạng thái của hệ tại một thời điểm được xác định bởi tọa độ và vận tốc của từng hạt, với sự chuyển động được
chi phối bởi định luật II Newton
𝐹𝑖(𝑡)= 𝑚𝑖𝑎𝑖(𝑡)= 𝑚𝑖𝑑2𝑟𝑖(𝑡)
𝑑𝑡2 = 𝐹𝑖= −∇𝑖𝐸 (r )Trong đó Fi là lực, mi là khối lượng, và ri là tọa độ của hạt i.
Quan hệ Lực - Thế năng: Lực không phải là một tham số ngẫu nhiên mà được dẫn xuất trực tiếp từ bề mặt thế
năng của hệ điều này có nghĩa là các nguyên tử sẽ chịu lực đẩy chúng về phía có mức năng lượng thấp hơn để
đạt trạng thái bền vững 𝐹𝑖= −∇𝑖𝐸 (r )
Trong đó, lực Fi được sinh ra từ trường thế năng E (Trường lực). Vì lực xác định gia tốc, và gia tốc xác định sự
thay đổi vị trí, nên nếu biết trạng thái ban đầu, ta có thể xác định được quỹ đạo của hệ tại mọi thời điểm.
Mô hình hóa tương tác: Trường lực (Force Field) Để giải phương trình chuyển động, thành phần quan trọng
nhất là hàm thế năng E(r). Trong MD, hàm này được gọi là Trường lực (Force Field), bao gồm tổng hợp của các
tương tác liên kết và không liên kết 𝐸(𝑟)= 𝐸𝑠𝑡𝑟𝑒𝑡𝑐ℎ +𝐸𝑏𝑒𝑛𝑑 +𝐸𝑑𝑖ℎ𝑒𝑑𝑟𝑎𝑙 +𝐸𝑒𝑙𝑒𝑐𝑡𝑟𝑜𝑠𝑡𝑎𝑡𝑖𝑐 +𝐸𝑣𝑑𝑊
Trong đó, các tương tác liên kết bao gồm năng lượng co giãn liên kết Estretch, uốn góc Ebend và xoắn Edihedral giúp
định hình cấu trúc phân tử. Các tương tác không liên kết bao gồm lực tĩnh điện Eelectrostatic và lực Van der Waals
EvdW, đóng vai trò quan trọng trong việc mô tả sự sắp xếp của các phân tử trong không gian. Các tham số cho các
hàm này được cung cấp bởi các bộ trường lực phổ biến như AMBER, CHARMM, GROMOS hoặc OPLS.
Cơ học thống kê: Cầu nối vi mô - vĩ mô Để kết nối các thông tin vi mô (tọa độ, vận tốc) thu được từ mô phỏng
với các tính chất vĩ mô quan sát được (như nhiệt độ, áp suất) MD dựa trên nền tảng của Cơ học thống kê. Hệ
thống được mô tả trong một không gian pha, và theo giả thuyết Ergodic, giá trị trung bình của một đại lượng
vật lý theo thời gian mô phỏng dài sẽ tương đương với giá trị trung bình trên tập hợp thống kê ⟨𝐴⟩𝐸𝑛𝑠 =⟨𝐴⟩𝑇𝑖𝑚𝑒
Thay vì phải tính tích phân phức tạp trên toàn bộ không gian pha, MD thực hiện việc "lấy mẫu" (sampling) thông
qua việc tạo ra một quỹ đạo các cấu hình theo thời gian để ước lượng các tính chất nhiệt động
𝐴=1
𝑁𝑠𝑡𝑒𝑝 ∑ 𝐴({𝑟𝑖(𝑡)})
𝑁𝑠𝑡𝑒𝑝
𝑡=1
Để mô phỏng sát với thực tế, hệ thống thường được đặt trong các tập hợp thống kê với các điều kiện biên cụ thể
như Vi chính tắc (NVE - năng lượng không đổi), Chính tắc (NVT - nhiệt độ không đổi) hoặc Đẳng nhiệt - Đẳng
áp (NpT).
Giải thuật tích phân số Phương trình chuyển động của hệ nhiều hạt là một hệ phương trình vi phân phức tạp
không có lời giải giải tích chính xác. MD giải quyết vấn đề này bằng phương pháp sai phân hữu hạn, chia nhỏ
thời gian thành các bước tin t. các thuật toán phổ biến được sử dụng giải thuật Verlet, Giải thuật Nhảy cóc (Leap-
frog).
Yêu cầu tiên quyết của các thuật toán tích phân trong MD là phải đảm bảo tính thuận nghịch theo thời gian
(time-reversibility). Đây là tính đối xứng cơ bản của các phương trình Hamilton/Newton, đảm bảo định luật bảo
toàn năng lượng được tuân thủ trong suốt quá trình mô phỏng dài.
Câu 3 Nguyên lý nào trong vật lý thống kê là cơ sở của ý tưởng phát triển phương pháp mô phỏng động
lực học phân tử ? Phát biểu nguyên lý và trình bày những hiểu biết của em về mối tương quan giữa hệ vi
mô và hệ vi mô thông qua nguyên lý đó .
Nguyên lý đó chính là nguyên lý Ergodic: Nguyên lý Ergodic phát biểu rằng: "Giá trị trung bình theo thời gian
của một đại lượng vật lý (khi quan sát hệ thống trong một thời gian đủ dài) sẽ tương đương với giá trị trung bình
theo tập hợp thống kê của đại lượng đó: ⟨𝐴⟩𝐸𝑛𝑠 =⟨𝐴⟩𝑇𝑖𝑚𝑒 Ens là trung bình theo tập hợp, Time trung bình theo
thời gian
Trong thực tế, chúng ta đo đạc các tính chất vĩ mô (Nhiệt độ T, Áp suất P, Thể tích V, Năng lượng nội tại U...).
Tuy nhiên, mô phỏng máy tính lại làm việc với các chi tiết vi mô (tọa độ r, vận tốc v của từng nguyên tử).
Hệ vi mô: Được mô tả bởi các trạng thái vi mô, tức là tập hợp tọa độ (r) và xung lượng (p) của tất cả các nguyên
tử tại một thời điểm. Tập hợp tất cả các trạng thái vi mô khả dĩ tạo thành một Không gian pha. Theo lý thuyết
thống kê thuần túy, để tính tính chất vĩ mô, ta phải tính tích phân trên toàn bộ không gian pha này (điều này bất
khả thi trong thực tế vì không gian pha quá lớn.
Mô phỏng MD tạo ra một quỹ đạo pha (phase trajectory) bằng cách giải phương trình chuyển động theo thời
gian. Nhờ Nguyên lý Ergodic, thay vì phải lấy mẫu trên toàn bộ không gian pha tĩnh (Ensemble), ta chỉ cần lấy
trung bình các giá trị của hệ trên quỹ đạo mô phỏng theo thời gian (Time average).
Từ đó, những con số vô tri về vị trí và vận tốc của từng nguyên tử (vi mô) được chuyển hóa thành các đại lượng
đo đạc được (vĩ mô) với độ chính xác phụ thuộc vào thời gian mô phỏng.
Tại sao kết quả được chấp nhận kết quả thực tế
Trong thực tế, các thí nghiệm không diễn ra trong môi trường ngẫu nhiên mà luôn chịu sự ràng buộc của các điều
kiện môi trường (nhiệt độ phòng, áp suất khí quyển, bình kín...). Để mô phỏng máy tính phản ánh đúng hiện tượng
thực tế đó, ta phải áp đặt các điều kiện biên tương ứng thông qua các tập hợp thống kê (Ensembles).
Câu 4 Trong mô phỏng MD trường lực có quan trọng không , tại sao? Trường lực cung cấp nhữung
thông tin gì , Hãy nêu các thành phần cơ bản của một trường lực mô tả tương tác giữa các nguyên tử
trong một hệ sinh học.
Trường lực (Force Field) đóng vai trò là thành phần cốt lõi, quyết định sự thành bại của một phép mô phỏng
Động lực học Phân tử (MD). Tầm quan trọng này được thể hiện qua ba khía cạnh chính sau đây:
+ Thứ nhất, là cơ sở vật lý sinh ra chuyển động.
Trường lực cung cấp hàm thế năng (U) để tính toán lực tác dụng lên từng nguyên tử thông qua đạo hàm. Đây là
dữ liệu đầu vào bắt buộc để giải phương trình Newton (F=ma). Nếu không có trường lực, hệ thống không thể
xác định gia tốc và quỹ đạo chuyển động của các hạt.
+ Thứ hai, quyết định độ tin cậy của kết quả
Trường lực chứa các tham số thực nghiệm mô tả bản chất hóa học của vật liệu. Nếu bộ tham số này không chính
xác, kết quả mô phỏng sẽ sai lệch hoàn toàn so với thực tế, bất kể thuật toán tính toán có tốt đến đâu. Do đó,
việc chọn đúng trường lực là yếu tố tiên quyết để kết quả mô phỏng có giá trị khoa học.
+ Thứ ba, giải quyết bài toán hiệu năng tính toán.
Trường lực là sự xấp xỉ hiệu quả thay thế cho các phương trình cơ học lượng tử phức tạp. Bằng cách đơn giản
hóa tương tác electron thành các mô hình cơ học cổ điển, trường lực cho phép mô phỏng các hệ thống lớn (hàng
triệu nguyên tử) trong thời gian dài mà vẫn đảm bảo độ chính xác chấp nhận được.
Trường lực cung cấp thông tin :
+ Hệ thống các hàm toán học : Đây là các phương trình giải tích định nghĩa cách tính thế năng tương tác
dựa trên tọa độ của các nguyên tử. Nó quy định cụ thể mô hình vật lý nào được sử dụng cho từng loại
tương tác, Các hàm này chính là "luật" để máy tính chuyển đổi vị trí không gian thành giá trị năng
lượng.
+ Cung cấp bộ cơ sở dữ liệu các tham số thực nghiệm : Đây là phần thông tin quan trọng nhất, chứa các
hằng số vật lý cụ thể cho từng loại nguyên tử. Các tham số này bao gồm: khối lượng nguyên tử, điện
tích riêng phần (q), độ dài liên kết cân bằng (r0), góc cân bằng (theta0), các hằng số lực đàn hồi (k), và
các tham số đặc trưng cho kích thước và độ sâu hố thế Van der Waals. Các giá trị này thường được xác
định thông qua thực nghiệm hoặc tính toán cơ học lượng tử độ chính xác cao.
+ Cung cấp quy tắc định danh và phân loại nguyên tử: Trường lực phân biệt các nguyên tử không chỉ dựa
trên nguyên tố hóa học mà còn dựa trên môi trường lai hóa của chúng. Ví dụ, trường lực sẽ cung cấp
thông tin để phân biệt nguyên tử Carbon trong nhóm methyl (sp3) khác với nguyên tử Carbon trong vòng
thơm (sp2) hay trong nhóm carbonyl. Điều này giúp máy tính áp dụng đúng các tham số vật lý cho từng
ngữ cảnh hóa học cụ thể, đảm bảo tính chính xác của mô hình.
Thành phần trường lực :
Trường lực (Force Field), bao gồm tổng hợp của các tương tác liên kết và không liên kết 𝐸(𝑟)= 𝐸𝑠𝑡𝑟𝑒𝑡𝑐ℎ +
𝐸𝑏𝑒𝑛𝑑 +𝐸𝑑𝑖ℎ𝑒𝑑𝑟𝑎𝑙 +𝐸𝑒𝑙𝑒𝑐𝑡𝑟𝑜𝑠𝑡𝑎𝑡𝑖𝑐 +𝐸𝑣𝑑𝑊
Trong đó, các tương tác liên kết bao gồm năng lượng co giãn liên kết Estretch, uốn góc Ebend và xoắn Edihedral giúp
định hình cấu trúc phân tử.
+ Thế năng co giãn liên kết: Mô tả dao động đàn hồi của độ dài liên kết giữa hai nguyên tử nối trực tiếp.
Thành phần này thường được mô hình hóa bằng định luật Hooke (lò xo), giữ cho độ dài liên kết dao động
quanh giá trị cân bằng.
+ Thế năng uốn góc: Mô tả năng lượng liên quan đến sự thay đổi góc hóa trị tạo bởi ba nguyên tử liên tiếp.
Thành phần này giúp duy trì cấu trúc lai hóa không gian của nguyên tử.
+ Thế năng góc xoắn: Mô tả chuyển động quay quanh trục liên kết. Đây là thành phần quan trọng nhất
trong hệ sinh học, quyết định cấu hình xoắn alpha hay phiến gấp beta của protein.
Các tương tác không liên kết bao gồm lực tĩnh điện Eelectrostatic và lực Van der Waals EvdW, đóng vai trò quan trọng
trong việc mô tả sự sắp xếp của các phân tử trong không gian.
+ Tương tác Van der Waals: bao gồm lực đẩy Pauli (giữ cho các nguyên tử không chồng lấn lên nhau) và
lực hút phân tán London (giúp ổn định phần lõi kỵ nước của protein).
+ Tương tác Tĩnh điện : Được mô tả bằng định luật Coulomb giữa các điện tích điểm riêng phần. Đây là
tương tác tầm xa, chịu trách nhiệm cho việc hình thành các liên kết hydro, cầu muối (salt bridges) và định
hướng phân tử trong môi trường nước.
Các tham số cho các hàm này được cung cấp bởi các bộ trường lực phổ biến như AMBER, CHARMM, GROMOS
hoặc OPLS.
Câu 5 Liệt kê các điều kiện biên trong mô phỏng mà em biết và nêu rõ điều kiện biên nào sẽ sử dụng cho
mục đích mô phỏng nào
1. Biên tuần hoàn : Đây là thủ thuật sao chép khối vật chất thành vô hạn các bản sao bao quanh mô hình gốc.
Nguyên tắc là nếu một nguyên tử đi ra khỏi biên bên này thì nó (hoặc ảnh của nó) sẽ đi vào từ biên đối diện.
• Mục đích sử dụng: Sử dụng để mô phỏng vật liệu khối Giúp loại bỏ các hiệu ứng bề mặt (surface
effects), mô phỏng hệ thống như một mạng tuần hoàn vô tận, giúp tính chất của hệ tiến gần với thực tế
hơn khi số lượng nguyên tử đủ lớn.
2. Biên tuần hoàn theo hai hướng và không tuần hoàn theo một hướng
Đây là dạng điều kiện biên hỗn hợp: tuần hoàn theo trục x, y nhưng không tuần hoàn theo trục z. Có 4 cách áp
dụng tùy mục đích:
- Biên tự do cho cả hai đầu trục z:
o Mục đích: Dùng để mô phỏng vật liệu dạng màng mỏng tự do. Giúp khảo sát cấu trúc, tính chất
của màng mỏng hoặc ảnh hưởng của bề mặt tự do lên vật liệu.
- Biên giả tự do theo hai đầu trục z:
o Mục đích: Tương tự biên tự do nhưng thêm một khoảng chân không đủ lớn để tránh lỗi tràn bộ
nhớ, thường dùng cho hệ có tương tác đơn giản để tiết kiệm chi phí tính toán.
- Biên tự do/giả tự do ở một đầu trục z + Biên cố định ở đầu kia:
o Mục đích: Dùng để mô phỏng màng mỏng tựa lên một đế cố định (ví dụ: màng mỏng nằm trên
đế tinh thể). Một đầu là bề mặt tự do, đầu kia bị giữ cố định bởi lớp đế.
- Biên cứng (Hard wall) / Phản xạ đàn hồi:
o Mục đích: Mô phỏng màng mỏng bị giới hạn giữa hai bức tường cứng.
3. Biên không tuần hoàn cho cả ba hướng (Biên tự do 3D)
Thường dùng cho các hệ trong không gian bị hạn chế hoặc vật liệu nano.
- Biên tự do:
- Mục đích: Dùng để khảo sát hạt nano hoặc các cụm nguyên tử cô lập
- Giúp khảo sát ảnh hưởng của bề mặt tự do lên cấu trúc hạt nano, hoặc mô phỏng các quá trình
chuyển pha như sự bay hơi (vì nguyên tử có thể thoát ra khỏi biên).
- Biên cứng :
- Mục đích: Mô phỏng vật liệu bị lưu giữ trong không gian hạn chế (ví dụ: khí trong bình kín, vật
liệu trong lỗ trống). Hạt sẽ va chạm và bật ngược lại khi gặp biên.
- Biên giả tự do : Mô hình được đặt trong không gian tính toán có biên cố định có ưu điểm giúp chống
tràn bộ nhớ khi tính toán.
- Biên có tương tác : Điều kiện biên có tương tác giữa các nguyên tử khác loại, đây là biên tiếp xúc giữa
2 vật liệu khác nhau
- Mục đích : Dùng mô phỏng màng mỏng đa lớp , tương tác giữa hai loại vật liệu khác nhau, khảo sát bề
mặt tiếp xúc giữa 2 loại vật liệu. Vd giữa H20 và SiO2, H2O và TiO2
Câu 6 Nêu các chế độ mô phỏng và các thống số, làm sao để duy trì toàn bộ áp suất , nhiệt độ trong quá
trình ?
1. Chế độ vi chính tắc NVE: Đây là chế độ cơ bản và tự nhiên nhất của phương pháp MD khi giải phương
trình Newton mà không có sự can thiệp từ bên ngoài. Tương tự quá trình đoạn nhiệt Hệ cô lập và không có trao
đổi nhiệt với bên ngoài.
Cần xác định vận tốc ban đầu và bước thời gian để thiết lập hệ
2.Chế độ NVT : Đây là chế độ mô phỏng hệ thống tiếp xúc với một bể nhiệt, thường dùng để nghiên cứu cấu
trúc vật liệu tại một nhiệt độ xác định. Hệ không cô lập, nhiệt được nhận vào hay lấy đi từ biên của hệ bằng
nhiều cách tùy thuộc cơ chế hoạt động của bể nhiệt (thermostat) .
Điều chỉnh nhiệt độ: bằng cách điều chỉnh vận tốc hoặc dùng các cơ chế bể nhiệt như Nosé-Hoover
thermostat, Berendsen thermostat, Andersen thermostat, Langevin dynamics.