Gaussian几何优化达到四个YES仍继续优化问题

Qing · 发表于 Post on 2022-5-6 14:04:04

sob老师，我在使用Gaussian进行几何优化时，出现了达到四个YES后没有显示找到不动点，仍然在继续优化的问题。体系是一个过渡金属阴离子（FeCl4-）和两个有机磷酸酯分子（C12H27O4P）结合，初始构型见附图，采用的关键词为#p opt(maxstep=5,notrust,cartesian) freq b3lyp/gen scrf=(smd,solvent=tributylphosphate) em=gd3bj，基组为ma-SVP（从您的博文《给def2以ma-方式加弥散函数的Gaussian格式的基组定义文件（含所有def2支持的元素）》中下载得到）。具体信息可以在log文件的11283看到（已上传至附图），出现了四个YES，但是没有显示找到不动点，优化仍在继续。Gaussian输入文件，输出文件和基组文件都已经上传到附件，麻烦老师看一下是怎么回事，谢谢sob老师！

ionexchangeC · 发表于 Post on 2022-5-6 15:00:11

优化到这一步时Hessian矩阵仍有本征值为负。
目测是SMD溶剂模型的锅（有可能导致假虚频），换成默认的PCM试一试

Qing · 发表于 Post on 2022-5-6 21:37:29

ionexchangeC 发表于 2022-5-6 15:00
优化到这一步时Hessian矩阵仍有本征值为负。
目测是SMD溶剂模型的锅（有可能导致假虚频），换成默认的PCM ...

请问老师如果采用PCM收敛了之后，还需要以该几何结构和波函数为初猜在SMD下计算到收敛吗？

wzkchem5 · 发表于 Post on 2022-5-6 22:06:41

Qing 发表于 2022-5-6 14:37
请问老师如果采用PCM收敛了之后，还需要以该几何结构和波函数为初猜在SMD下计算到收敛吗？

不需要，PCM对于结构优化来说精度一般够了

Qing · 发表于 Post on 2022-5-6 22:20:23

wzkchem5 发表于 2022-5-6 22:06
不需要，PCM对于结构优化来说精度一般够了

好的老师，还有一个问题，因为我想要计算的是A+B=C这个过程的自由能变化，A和B都已经在SMD模型下收敛了，而且得到了自由能校正值，C就是现在有问题的这个构型，如果对它在PCM下优化，那产生的自由能校正项，和前面A和B在SMD计算得到的自由能校正项，能同时带入到各个单体的总自由能计算，然后进行做差吗？

wzkchem5 · 发表于 Post on 2022-5-6 23:01:25

Qing 发表于 2022-5-6 15:20
好的老师，还有一个问题，因为我想要计算的是A+B=C这个过程的自由能变化，A和B都已经在SMD模型下收敛了， ...

不能，必须统一用同一个溶剂模型。要么仍然用SMD继续试试消虚频，要么把A和B都换成用PCM优化

Qing · 发表于 Post on 2022-5-6 23:45:13

wzkchem5 发表于 2022-5-6 23:01
不能，必须统一用同一个溶剂模型。要么仍然用SMD继续试试消虚频，要么把A和B都换成用PCM优化

好的谢谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian几何优化达到四个YES仍继续优化问题

评分 Rate