计算化学公社

标题: 求助：优化后的小分子能量最小化提示受力过大的问题 [打印本页]

作者
Author: MGRACE 时间: 2024-11-28 17:11
标题: 求助：优化后的小分子能量最小化提示受力过大的问题
各位老师好，目的是要在水中模拟离子液体。首先使用Gaussian16优化PF6-分子，基组选择的是社长在《RESP2原子电荷的思想以及在Multiwfn中的计算》中推荐的基组，电荷使用resp2.0，拓扑参考sobtop例1生成，但在gromacs中提示F原子受力达到10^14，尝试进行npt时程序崩溃。不知道怎么调整才合理

作者
Author: sobereva 时间: 2024-11-28 17:23
要么RESP2、要么RESP2(0.5)，没有RESP 2.0

http://sobereva.com/soft/Sobtop里专门加入了以下文字：

注意：如果基于上面的DRIH方法产生拓扑文件进行模拟遇到报错，可以尝试后面提到的mSeminario或m2Seminario方法产生力常数，或者干脆让Sobtop自动猜谐振势形式的力场参数而不基于Hessian计算（即让你选择如何生成力场参数时，选择4 Same as option 3, but missing ones are arbitrarily guessed。反正对于这种刚性体系，力常数不重要，只要大小足够维持住体系的刚性就够了）。

作者
Author: MGRACE 时间: 2024-11-28 21:13

sobereva 发表于 2024-11-28 17:23
要么RESP2、要么RESP2(0.5)，没有RESP 2.0

http://sobereva.com/soft/Sobtop里专门加入了以下文字：

好的，谢谢社长，尝试m2Seminario生成力常数之后可以正常模拟

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