计算化学公社

标题: Gas和PCM优化问题 [打印本页]

作者
Author: CruiseBend 时间: 2020-10-21 09:24
标题: Gas和PCM优化问题
本帖最后由 CruiseBend 于 2020-10-21 10:17 编辑

老师们好，想问下我现在有一个带电体系（4,1），在gas下b3lyp-d3/6-311g**下优化无虚频，结构合理，而在pcm下从第3步开始结构就不合理，且最终存在虚频。问题如下：
1. 这种情况是什么问题？感觉不像是设置收敛限就能解决的
2. 是否能使用gas下的优化构型，pcm计算sp，以计算相互作用能。
PS-此体系是以静电作用为主导的。

作者
Author: biogon 时间: 2020-10-21 10:07
1 换个溶剂模型试试，也有可能结构在溶剂模型下不稳定了
2 溶剂模型影响这么大这么做就不一定能行

作者
Author: sobereva 时间: 2020-10-21 11:12
光从你的描述上没法准确猜测实际情况
整体带那么多净电荷，优化不带溶剂模型是绝对不能接受的，就算你觉得结果没问题，内行审稿人也不会答应

作者
Author: CruiseBend 时间: 2020-10-22 15:32

sobereva 发表于 2020-10-21 11:12
光从你的描述上没法准确猜测实际情况
整体带那么多净电荷，优化不带溶剂模型是绝对不能接受的，就算你觉得 ...

老师好，这是我的input文件，使用了opt=recalc=4消除虚频还是不行，您看要怎么处理

作者
Author: wzkchem5 时间: 2020-10-22 16:57

CruiseBend 发表于 2020-10-22 15:32
老师好，这是我的input文件，使用了opt=recalc=4消除虚频还是不行，您看要怎么处理

虚频有多大？
优化出来的结构合理吗？

作者
Author: CruiseBend 时间: 2020-10-22 16:59

wzkchem5 发表于 2020-10-22 16:57
虚频有多大？
优化出来的结构合理吗？

-30
还算合理

作者
Author: CruiseBend 时间: 2020-10-22 19:51
各位老师，是不是因为我在放置复合物时使用的是气态稳定的单体而不是PCM下稳定的单体造成的呢

作者
Author: wzkchem5 时间: 2020-10-23 18:12

CruiseBend 发表于 2020-10-22 16:59
-30
还算合理

-30完全可能是格点误差或溶剂模型的数值误差
可以试试提高格点精度

作者
Author: CruiseBend 时间: 2020-10-24 09:56

wzkchem5 发表于 2020-10-23 18:12
-30完全可能是格点误差或溶剂模型的数值误差
可以试试提高格点精度

老师好，我使用了int=superfine opt=tight 但还是有虚频

作者
Author: sobereva 时间: 2020-10-24 10:02

CruiseBend 发表于 2020-10-22 19:51
各位老师，是不是因为我在放置复合物时使用的是气态稳定的单体而不是PCM下稳定的单体造成的呢

跟着没关系。优化完了都一样
看虚频对应的是什么模式，微调试试。实在没辙就用calcall

作者
Author: CruiseBend 时间: 2020-10-27 14:56
老师们好，我使用了微调+opt-recalc-superfine还是有-30的虚频，请问要怎么处理呢？（每次都是频率计算时候位移不收敛）

作者
Author: CruiseBend 时间: 2020-10-28 14:53
有老师能解答下吗？不胜感激

作者
Author: sobereva 时间: 2020-11-1 13:34

CruiseBend 发表于 2020-10-27 14:56
老师们好，我使用了微调+opt-recalc-superfine还是有-30的虚频，请问要怎么处理呢？（每次都是频率计算时候 ...

加tight啊
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

另外，如果用的是G09 D.01，不要用BJ阻尼的D3校正

作者
Author: CruiseBend 时间: 2020-11-1 16:51

sobereva 发表于 2020-11-1 13:34
加tight啊
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（ht ...

加了tight
但是用的em=gd3bj，g16下的

作者
Author: sobereva 时间: 2020-11-1 17:32

CruiseBend 发表于 2020-11-1 16:51
加了tight
但是用的em=gd3bj，g16下的

解决办法我博文里都写了，没有遗漏。反复阅读

作者
Author: CruiseBend 时间: 2020-11-2 11:12

sobereva 发表于 2020-11-1 17:32
解决办法我博文里都写了，没有遗漏。反复阅读

方法1和方法2与微调我都尝试了

作者
Author: sobereva 时间: 2020-11-2 13:07

CruiseBend 发表于 2020-11-2 11:12
方法1和方法2与微调我都尝试了

怎么就不考虑calcall

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