计算化学公社

标题: GAFF力场对200℃的体系是否适用 [打印本页]

作者
Author: Jiang_Yi 时间: 2025-11-12 14:47
标题: GAFF力场对200℃的体系是否适用
如题，我想要研究在200℃下有机分子在铁表面吸附行为，有机分子用GAFF力场描述，铁的参数来自文献（npj Comput Mater 7, 17 (2021). https://doi.org/10.1038/s41524-020-00478-1），但是对此我有一些困惑：
1、GAFF力场可以用于200℃的体系的研究吗？包括吸附能、径向分布函数、质量浓度分布和均方位移。我担心的是会不会描述不好。
2、如果不可以，还有什么方法可以研究吗？
希望各位老师不吝赐教，谢谢

作者
Author: 牧生 时间: 2025-11-12 15:59
本帖最后由牧生于 2025-11-12 16:01 编辑

第一：那篇文献中，Parameters are given for fcc γ-Fe in stainless steel alloys and can be used for bcc α-Fe. In the simulation of bcc α-Fe, atomic positions should be fixed to prevent a transition into fcc γ-Fe.
我一直怀疑这里是写错了，应该是can not be，要不然这句话也读不通啊。所以，如果你要做200℃的铁，要找 bcc α-Fe的铁胞，并fix这些铁原子（我建议用冻结）

第二：GAFF应该不太适用于200度的高温，但我觉得也不会有定性错误。你可以尝试一下就直接就用GAFF

作者
Author: Jiang_Yi 时间: 2025-11-12 20:30

牧生发表于 2025-11-12 15:59
第一：那篇文献中，Parameters are given for fcc γ-Fe in stainless steel alloys and can be used for b ...

明白了，谢谢老师

作者
Author: sobereva 时间: 2025-11-13 01:25
没有主流有机分子力场专门考虑高温的模拟。先跑起来，如果觉得有问题再说。也可以以CP2K跑的AIMD的结果作为参照进行验证

作者
Author: Jiang_Yi 时间: 2025-11-13 09:12

sobereva 发表于 2025-11-13 01:25
没有主流有机分子力场专门考虑高温的模拟。先跑起来，如果觉得有问题再说。也可以以CP2K跑的AIMD的结果作为 ...

好的，谢谢sob老师。我先跑起看看

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