计算化学公社

标题: 求助，CP2K做碳控压降温时温度降下不去 [打印本页]

作者
Author: qaz707000195 时间: 2024-6-5 15:38
标题: 求助，CP2K做碳控压降温时温度降下不去
本帖最后由 qaz707000195 于 2024-6-5 15:52 编辑

目前使用CP2K软件做AIMD计算。将石墨烯片加热，然后降温形成无定形结构。升温采用GFN1x-TB算的，降温采用PBE算的，均采用CSVR热浴，降温时采用各项同性控压。在第二步冷却过程中，温度从4000K降到700K，跑了2000步温度变化如下图，感觉变化不大。 (, 下载次数 Times of downloads: 21)
目前有以下几个问题：（1）GFN1x-TB和PBE是否能够合理描述碳原子的断键和成键过程（只要没有重大bug就好）。（2）因为在做第二步的cool.inp文件时忘记了把&VELOCITYS部分从第一步得到的heat.reatart文件中复制过来了（修正后的cool.inp还在跑），不知道温度降不下来是没有&VELOCITYS部分的原因还是说跑的时间不够的原因。下面附带第一次跑的cool.inp（没有&VELOCITYS部分，后续已经添加了）。

作者
Author: sobereva 时间: 2024-6-5 16:06
GFN-xTB是否适合描述，看轨迹判断。记得开smearing

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)