求助，CP2K做碳控压降温时温度降下不去

qaz707000195 · 发表于 Post on 2024-6-5 15:38:31

本帖最后由 qaz707000195 于 2024-6-5 15:52 编辑

目前使用CP2K软件做AIMD计算。将石墨烯片加热，然后降温形成无定形结构。升温采用GFN1x-TB算的，降温采用PBE算的，均采用CSVR热浴，降温时采用各项同性控压。在第二步冷却过程中，温度从4000K降到700K，跑了2000步温度变化如下图，感觉变化不大。

目前有以下几个问题：（1）GFN1x-TB和PBE是否能够合理描述碳原子的断键和成键过程（只要没有重大bug就好）。（2）因为在做第二步的cool.inp文件时忘记了把&VELOCITYS部分从第一步得到的heat.reatart文件中复制过来了（修正后的cool.inp还在跑），不知道温度降不下来是没有&VELOCITYS部分的原因还是说跑的时间不够的原因。下面附带第一次跑的cool.inp（没有&VELOCITYS部分，后续已经添加了）。

sobereva · 发表于 Post on 2024-6-5 16:06:28

GFN-xTB是否适合描述，看轨迹判断。记得开smearing

跟初速度没关系。热浴的事

牵扯到成键断键的问题，TIMESTEP 1.0太大了，改成不超过0.5
减小热浴时间常数，以让热浴效果更显著。TIMECON 800这也太大了。改成50

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，CP2K做碳控压降温时温度降下不去