求助：关于用cp2k跑有机自由基分解的MD过程中XC泛函的选择

秋心

各位大佬，本人刚刚入坑cp2k做AIMD没多久，做了一个在35x35x35Å盒子中有机自由基分解（HOOCH2SC•=O，该自由基本身存在两种分解路径，一种是CO，能垒只有2.9kcal/mol,另一种是OCS,能垒很高）的真空相MD，由于刚开始并不是很了解xc,便直接选用了文献中很常用的BLYP泛函，得到的结果是自由基发生了分解（1ps左右），生成了OCS和H2C=O以及OH•。后面随着深入的了解，发现涉及过氧化物在界面上的光解，像H2O2和CH3OOH都是使用QM/MM做的，QM用的都是B3LYP来描述过氧化物，所以后面就补充了将cp2k中BLYP泛函改为B3LYP泛函的真空相MD，预料之中的计算量和耗时成倍加大，但是截至目前为止跑了6.9 ps有机自由基（HOOCH2SC•=O）都没有发生分解。这就让我对BLYP泛函得到的结果产生了怀疑，但是一个有机自由基并不稳定，况且分解成CO这条路径，能垒只有2.9kcal/mol,按照常理来说，这个分解路径应该是可以跑出来的。现在的想法是模拟时长不够？B3LYP和BLYP处理有机过氧化物的准确性？有没有大佬有类似的经验，希望不吝赐教一下，万分感谢！！！

sobereva

小分子的AIMD用CP2K效率极低，最好选择是用ORCA，看
使用ORCA做从头算动力学(AIMD)的简单例子
http://sobereva.com/576（http://bbs.keinsci.com/thread-20800-1-1.html）

秋心

sobereva 发表于 2022-11-2 10:06
小分子的AIMD用CP2K效率极低，最好选择是用ORCA，看
使用ORCA做从头算动力学(AIMD)的简单例子
http://sob ...

好的，谢谢sob老师，我去了解一下。

madhatter

sobereva 发表于 2022-11-2 10:06
小分子的AIMD用CP2K效率极低，最好选择是用ORCA，看
使用ORCA做从头算动力学(AIMD)的简单例子
http://sob ...

sob老师，我想问一下，例如我想模拟电池中锂电极表面硝酸根离子的分解并与锂反应生成新物质的这类反应动力学的研究，我使用pbe跑了10ps，并没有出现变化。目前思路一个是否是因为pbe的原因，需要选择其他泛函吗？还是需要把电解液的其他组分一起加入之后再去模拟？还是这类反应动力学需要设置其他参数？

wzkchem5

madhatter 发表于 2022-11-10 16:06
sob老师，我想问一下，例如我想模拟电池中锂电极表面硝酸根离子的分解并与锂反应生成新物质的这类反应动 ...

这个反应实验已知在10ps内能反应吗？做过瞬态光谱吗？
可能这个反应本来就不是ps级的时间尺度里能发生的，此时AIMD本来就该跑不出来反应，跑出来有反应反倒是错的