Inv3 failed in PCMMkU (L502/L508 报错) 的成因与解决方法

Aridea · 发表于 Post on 2019-10-9 12:45:09

前辈，用surface=sas 果然能够opt完成，但是这个结构能直接用吗——我想接下来算freq得到其Gibbs自由能校正项，随后计算高精度气相单点能，会不会不准。还请前辈赐教

liyuanhe211 · 发表于 Post on 2019-10-9 13:46:11

Aridea 发表于 2019-10-9 12:45
前辈，用surface=sas 果然能够opt完成，但是这个结构能直接用吗——我想接下来算freq得到其Gibbs自由能校正 ...

你拿几个两种都能运行完的比较一下就知道了

Aridea · 发表于 Post on 2019-10-9 19:44:18

liyuanhe211 发表于 2019-10-9 13:46
你拿几个两种都能运行完的比较一下就知道了

嗯额，谢谢前辈

sunjk · 发表于 Post on 2021-5-11 15:40:53

请问老师，我在优化结构时，加surface=sas的和不加的最终优化构型，金属键长差了超过0.1nm。也就是说，加了surface=sas后，对结构优化的结果是有影响。这样加了surface=sas后的优化结果，其构型，能量等是不是就不能和其它的没加的相比较了？

量化小菜鸡 · 发表于 Post on 2021-5-11 16:22:45

sunjk 发表于 2021-5-11 15:40
请问老师，我在优化结构时，加surface=sas的和不加的最终优化构型，金属键长差了超过0.1nm。也就是说，加了 ...

你比较能量和自由能矫正了吗，我也遇到这种问题。我最后优化清一色全改IEFPCM去优化了，最后SMD做单点。

liyuanhe211 · 发表于 Post on 2021-5-11 17:28:34

sunjk 发表于 2021-5-11 15:40
请问老师，我在优化结构时，加surface=sas的和不加的最终优化构型，金属键长差了超过0.1nm。也就是说，加了 ...

键长一共多长，差0.1nm还有键吗？

CCCC · 发表于 Post on 2021-11-22 09:17:38

g09 D.01也报了同样的错误，来试试前辈们的解决办法

pzl · 发表于 Post on 2022-2-20 02:12:42

用bp86(d)/def2svp 算了一个乙烷的，发现在sas孔洞模型的时候，能量高一点，可能是体积大了，相互极化的小一点。然后ses和vdw差不多，因为体积几乎没差

ionexchangeC · 发表于 Post on 2022-3-19 22:04:29

我在计算中发现溶剂空腔问题会导致原本高对称性的结构出现不合理的虚频。（B3LYP/def2-TZVP级别）在气相下优化[Ca(H2O)6]2+水合钙离子的结构，得到的结构为Th点群，振动分析后没有虚频；但改在溶剂模型下优化，得到的Th点群分子出现了6个虚频（使用SMD溶剂模型优化时则出现了21个虚频）。破坏对称性后继续优化，结构震荡难以收敛，尝试各种手段时偶然出现该报错，于是用GaussView打开输出文件观看溶剂化表面，发现取的空腔明显偏小（SMD则更小）。重新用高对称性的结构加上surface=SAS在溶剂模型优化，优化后Th点群保持并且没有虚频。也顺便试了下加上关键词surface=SES优化，结果在第一步优化的SCF迭代过程中就提示“Density matrix breaks symmetry”。

GoldenBaby · 发表于 Post on 2022-7-24 12:07:46

量化小菜鸡发表于 2021-5-11 16:22
你比较能量和自由能矫正了吗，我也遇到这种问题。我最后优化清一色全改IEFPCM去优化了，最后SMD做单点。

好像IEFPCM也会出现这种情况，目前个人觉得比较好的解决方案还是用能算的模型跳过不行的结构再接着切回来吧

laplacesuanzi · 发表于 Post on 2022-10-2 22:00:12

本帖最后由 laplacesuanzi 于 2022-10-2 22:15 编辑

我最近也遇到了相同的问题，一模一样的报错。关键词最开始为“#p opt freq b3lyp scrf=(solvent=generic,read) em=gd3bj tzvp”，通过参考http://sobereva.com/61，我在关键词中又加上了“SCF=noincfock”，这个问题就解决了。感谢卢天老师。

Winnie1998 · 发表于 Post on 2022-10-16 11:34:19

各位老师好，我在结构优化的时候出现了上述报错，结构优化已经进行到第58帧，请问是要选取最后一帧的结构修改关键词后继续优化吗？还是要重新开始？

JCenter · 发表于 Post on 2023-8-10 11:11:24

老师，我有一个疑问。当我进行结构优化时出现这个报错时，我在使用surface=sas完成结构优化时，还需要使用ses表面计算单点吗。我直接用sas优化的结构在更高的计算级别单点能，使用默认的表面不可以吗。

S221100847 · 发表于 Post on 2023-12-28 17:21:37

您好，想请教下，我这样写这个功能可以被使用吗

zhang123 · 发表于 Post on 2024-5-9 10:00:19

您好，请问您的问题解决了嘛？我是用gaussian 16B.01算的单点能，输入文件如下：
%oldchk=Co_P_2_1_gas-sp.chk
%chk=Co_P_2_1_solv-sp.chk
#p ub3lyp/6-31+g* guess=read geom=allcheck scrf=g03defaults，也出现了一样问题，我看了一下G03默认的就是SES表面，不知道为啥还会报错呢？求大佬指点一下

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

[Gaussian/gview] Inv3 failed in PCMMkU (L502/L508 报错) 的成因与解决方法