牧生 发表于 2022-3-26 09:16
你这个结构式，确保是全原子模型么？
我认为也有可能是初始结构不合理导致，两个原子如果距离比较近，或被 ...

谢谢您的回复，我这个体系因为组装结构的原因导致有的烷烃链上的氢会挨得比较近，我试试修改一下模型。

Rosefinch 发表于 2022-3-14 22:56
好的，谢谢sob老师，我那无法识别的二面角都集中在下图这6个红圈内，应该是可以视为刚性的。

你这个结构式，确保是全原子模型么？
我认为也有可能是初始结构不合理导致，两个原子如果距离比较近，或被认为成键。画结构时，把两个原子距离拉远一点再优化，不然可能导致程序认错。优化后的位置和距离就合理了。
我有过类似的错误，看第五楼，和你现在的情况有点相似。
http://bbs.keinsci.com/thread-28435-1-1.html   

Rosefinch 发表于 2022-3-25 22:04
谢谢sob老师，那请问这个力常数一般设置为多少能够刚好保持刚性又不会使环受力太大导致体系崩溃呢？
模 ...

没具体文件没法说
检查拓扑文件里相关参数
用1fs步长，不设约束

Rosefinch 发表于 2022-3-25 12:28
还想请问一下sob老师：
1. 对于这种大环分子该如何去获取它的原子电荷？总共500+的原子，只有C、S、O、H ...

1 “选择单点能”语义不明

取出一个大环分子，把烃链短，构造较小的模型体系，然后算RESP电荷。用懒人脚本里的级别就可以
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）

然后再构造个模型体系（不考虑所有的环），带着完整的烃链，再计算RESP电荷，级别同上。

以上两次得到的原子电荷恰当拼接构成单个完整环状+烃链的原子电荷。算RESP电荷期间注意恰当设置片段电荷约束，以便于完美拼接，做法参考下文。
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

如果追求傻瓜化，用Sobtop自动调用Openbabel算MMFF94电荷，但质量比RESP有明显差距。

2 跟原子电荷没关系。参数问题

Rosefinch 发表于 2022-3-14 22:56
好的，谢谢sob老师，我那无法识别的二面角都集中在下图这6个红圈内，应该是可以视为刚性的。

这肯定都是刚性的
在Sobtop的bonded_param.dat的$HARM_DIH里手动加入缺失的力场项，力常数随便设一个足够大的能维持刚性的值即可。

Rosefinch 发表于 2022-3-13 13:15
谢谢您的回复，我利用sobtop生成了完全基于GAFF力场的文件，但是在二面角那里没有生成参数，显示missing  ...

看具体缺哪些，没结构和拓扑文件别人没法说
缺的二面角可以自己弄个模型体系做扫描并拟合得到，论坛上ztop可以拟合（以后sobtop也可能直接内置拟合功能），对于不重要的二面角可以凑合一个。如果二面角是刚性的就更简单了，用优化后的结构当平衡角度，力常数随便猜一个够大的即可。

sobtop没二面角的话就用acpype获得拓扑，别让它优化，自己算电荷

Rosefinch 发表于 2022-3-13 13:16
还想请问一下，如果我的有机分子采用了GAFF力场，那么两层石墨烯基底是否也需要用sobtop来生成GAFF力场的参 ...

1、如果你能看得懂这个帖子，就按照这个来   http://sobereva.com/635   但这个也有一个问题，就是你用的分子中，原子数较多，量化计算过程耗时间比较长。
2、如果想要简单易行，就参照sobtop的第二个例子 ，如果有些参数missing，还得要回到上一个的操作
3、原则上，都建议使用相同的力场，amber和gaff是兼容的。

还想请问一下，如果我的有机分子采用了GAFF力场，那么两层石墨烯基底是否也需要用sobtop来生成GAFF力场的参数？

牧生 发表于 2022-3-12 17:46
考虑一下 http://sobereva.com/soft/Sobtop

谢谢您的回复，我利用sobtop生成了完全基于GAFF力场的文件，但是在二面角那里没有生成参数，显示missing and treated as free rotor，请问我该如何填补参数，或者是否需要填补

考虑一下 http://sobereva.com/soft/Sobtop

What is Sobtop?

Very very very easily generating forcefield parameters and GROMACS topology file. Applicable to almost any system!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!