AIMD và tỷ lệ vận tốc trong VASP

Question

Thông thường người ta thực hiện chia tỷ lệ vận tốc để "làm nóng" mẫu của họ với MD.

1. Làm thế nào để chúng tôi thực hiện điều này trong VASP? Có phải bằng cách đặt các thẻ SMASS=-1và MDALGO=0không?

2. Tại sao chúng ta không thể "sưởi ấm" bằng cách sử dụng máy điều nhiệt Nose-Hoover? SMASS=0, MDALGO=2?

* khi tôi thử cả hai phương pháp "gia nhiệt" (100 đến 600 K), tôi thấy tổng E tăng trong phương pháp 1 và giảm trong phương pháp 2. Tại sao lại như vậy? Tôi mong đợi năng lượng sẽ tăng lên kể từ khi tôi bắt đầu ở trạng thái cơ bản của mình. Nhân tiện, tôi đang sử dụng một nhóm NVT ( ISIF=2).

Chỉnh sửa: Tôi nghĩ rằng bài đăng của tôi đã quá lâu trước đó nên tôi đã đơn giản hóa nó

Answer 1

Tôi nghĩ rằng tôi đã tìm ra một phần câu hỏi của mình. Nhưng thực sự sẽ được hưởng lợi nếu có chuyên gia VASP / AIMD tham gia.

1. Đây dường như là cách chính xác để thực hiện quy mô vận tốc trong VASP. Và bạn có thể chia tỷ lệ vận tốc mỗi NBLOCKbước (nhưng tôi không chắc liệu / khi nào thì thích hợp hơn để chia tỷ lệ mỗi 1 bước, 10 bước, v.v.)

2. Về mặt kỹ thuật, bạn có thể tôi nghĩ. Tôi đặt các thông số được liệt kê trong câu hỏi của mình. Và cũng như TEBEG = 100và TEEND = 600. Nhiệt độ và năng lượng tăng lên. Tôi chỉ không biết liệu điều này có chính xác để "làm nóng" theo cách này hay không vì tôi thường thấy mọi người chỉ điều chỉnh nhiệt độ không đổi với bộ điều chỉnh nhiệt Nose-Hoover và bộ ghép NVT

* tổng E thực tế giảm trong phương pháp 2 vì thế năng của bồn tắm giảm. Vì vậy, về mặt kỹ thuật, tôi chỉ nên nhìn vào E0 + EK (thế năng + động năng), sẽ tăng khi "nhiệt độ tăng".