分子动力学aimd_CP2K教程系列之第一性原理分子动力学(Pymatflow篇)

本系列CP2K教程是《CP2K菜根谭》的升级版，在旧版基础上添加了如何结合Pymatflow工具简化计算流程的内容。话不多说，本文将为您带来CP2K系列教程中的第一性原理分子动力学(AIMD)部分。

简介

CP2K支持第一性原理分子动力学(AIMD)以及混合量子经典分子动力学(QM/MM)。除了波恩-奥本海默分子动力学以及艾伦菲斯特分子动力学以外，CP2K还实现了一种类似于传统的Car-Parrinello分子动力学的第二代Car-Parrinello分子动力学，引起了许多学者的关注。

下面的内容，我们将对一个基本的AIMD计算进行示例分析，以让大家熟悉如何使用CP2K进行AIMD模拟。

使用CP2K进行分子动力学模拟的输入文件与进行几何结构优化的输入文件类似，同样都有三个大的一级输入模块GLOBAL、FORCE_EVAL与MOTION来组成。当然关于AIMD的参数设置主要位于MOTION模块中，而GLOBAL中通常只需要将GLOBAL/RUN_TYPE设置为MD即可。所以我们就MOTION中的设置进行分析。

下面是一个用于分子动力学模拟的MOTION模块的样例:

&MOTION
    &MD
        ENSEMBLE NVE
        STEPS 1000
        TEMPERATURE 300
        TEMP_TOL 0.0
        TIMESTEP 0.5
    &END MD
    &PRINT
    &END PRINT
&END MOTION