使用reaxff力场模拟C=N双键的形成该如何设置控制参数？

mtpeng

大家好，本人想模拟氨基和醛之间的席夫碱反应，Lammps输入文件是通过packmol在3×3×3盒子中放入两种反应分子各50个，随后利用VMD转的data文件，先遇到如下问题：
1.species结果表明没有形成C=N双键，甚至输入的初始结构分子在species都没有出现过。再者，输入的反应物单体仅含有一个氨基（N1）,当species结果偏向出现N3分子而没有N1;
2.不知建模方式是否有问题才导致无法正常形成初始结构？
3.该如何调整键级判断标准以促使形成C=N双键？

wildon

得确保力场文件用的合适

mtpeng

wildon 发表于 2024-7-25 11:25
得确保力场文件用的合适

你好，再请问一下我是在lammps自带的力场里找了个包含C, H, N , O的力场，这难道不是通用的吗？