关于cp2k从头算体系平衡判断、密度计算等一些基本问题的请教

Charlie-Ma

大家好，最近刚开始学习cp2k，想用cp2k做从头算分子动力学模拟，想请教大家一些基本的问题。

1. 如何判断从头算模拟体系是否达到了平衡，有哪些方法或者输出参数可以参考判断？比如说电子步稳定在5步以内收敛是否可以说体系达到了平衡？

2. 由于很多时候不知道某体系在某一温度下的实际密度？如果建模时的密度高于或低于体系的实际密度，对平衡结构的准确性影响大吗？

3. 从头算的输出参数里哪里可以看到体系的密度，压强，模拟盒子边长等参数呢？lammps中可以通过thermo_style控制输出到log文件中，但是不知道cp2k咋输出？

4. 关于&MGRID中CUTOFF的选取，由于体系中含有O，是不是需要选取大于等于800Ry的值呢？选取这么大的值对计算速度会不会影响很大？

5. cp2k做从头算一般计算多少原子的体系？64核的服务器，考虑到精确性和效率，大概能跑多大体系合适呢？

6. 请问论坛里有没有用从头算做玻璃体系的同学，能否加个好友请教学习一下？

谢谢大家！

sobereva

1 具体拿什么说事看具体体系
诸如盒子体积或密度、势能、温度等都可以拿来作为判据。不清楚的话就多看点计算文章

2 控压，让盒子尺寸自发调节就完了

3 输出文件里就有盒子信息

4 没法一概而论，跟你用的基组关系很大，而且多大也看做什么计算

5 看具体计算级别，跑多长时间，体系是什么构成。没法简单说。多看些近几年CP2K跑AIMD的文章就大概有感觉了