求助：碳化钙与氢气反应的分子动力学模拟用什么软件合适呢？

含光君

Gromacs不能模拟反应力场，可考虑在Lammps中用reaxff势模拟。

sobereva

得用CP2K等跑AIMD，或者按2L说的用反应力场。gromacs支持的力场里没有能描述反应的，对这种体系顶多算个物理吸附之类。

含光君

keyanxiaobai 发表于 2022-9-22 11:34
请问有没有Lammps计算反应的相关资料吗，我对lammps不是很了解

分享几个我觉得不错的学习资源：
LAMMPS Manual (https://docs.lammps.org/Manual.html)
LAMMPS Tutorials (https://lammpstutorials.github.io/)
LAMMPS和Gromacs聚合物建模方法 (http://bbs.keinsci.com/thread-18520-1-1.html)
Introduction to LAMMPS (https://www.youtube.com/watch?v=QdygibMW1AI)
peachRL的小站 (https://peachrl.github.io/archives/)

cyx98

用lammps可以，ReaxFF反应力场，但是不一定有你这个Ca/C/H体系的力场。
力场可以在这个里面找 https://www.scm.com/doc/ReaxFF/Included_Forcefields.html
而且就算有，你也需要先验证力场是否符合你这个体系

wzkchem5

此外注意一下这个反应的速率。AIMD一般只能跑得动ps~ns级的动力学（到ns级可能得用超算大力出奇迹才行），用ReaxFF（或者AIMD用机器学习加速）可能可以跑到μs级，然后基本就到头了。但是很多反应都是秒级甚至小时级的时间尺度。如果你的反应在比如说1微秒之内反应不完，那么要么在高于实验温度的条件下模拟（但这样就只能定性看个大概，不好计算反应速率常数了，也不好估计不同反应机理的比例），要么做metadynamics（但这时怎么定义合适的CV就是个问题），要么用过渡态理论。

含光君

keyanxiaobai 发表于 2022-9-23 17:00
我去网站发现CaSiAlO.ff: (C/H/O/Fe/Cl/Si/Al/Ca)，请问这个力场可以用吗？

在网站上会标注力场的原文，文章中会写明研究的体系。网站下面的介绍也会简单说明是适用于溶液还是燃烧。

cyx98

keyanxiaobai 发表于 2022-9-23 17:00
我去网站发现CaSiAlO.ff: (C/H/O/Fe/Cl/Si/Al/Ca)，请问这个力场可以用吗？

你可以先验证下。

含光君

可以是肯定可以，但依然要看你选择的反应力场是否合适。