计算化学公社

标题: 使用reaxff力场模拟C=N双键的形成该如何设置控制参数？ [打印本页]

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

作者
Author: wildon 时间: 2024-7-25 11:25
得确保力场文件用的合适

作者
Author: mtpeng 时间: 2024-7-25 16:07

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

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

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