Gaussian 常见报错及解决方法；新手求助报错时的注意事项

clvn · 发表于 Post on 2017-4-16 11:15:17

老师，我看了您最后部分的通过做柔性扫描找过渡态，我有几个问题，还想请老师解决一下：
1.我用了您文中的opt=(calcfc,ts,noeigentest) freq这条命令，如您文中所述，最后的结构报错L9999,不收敛，但我不明白的是您在文中写到：“描时最左侧一点时会遇到了L9999 的 opt 未收敛错误“，而我做完opt=(calcfc,ts,noeigentest) freq只有一个结构，还是不收敛的，所以，我更加不明白老师文章里写的”此时打开输出文件可见路径上的最高点已经找到了“
2.您在文中说，近4步优化就找到了正确的过渡态，过渡态 C-H 键键长 1.675 A。我不是很懂您的4步是怎么看出来的
以上是我的问题，还请老师百忙之中，解答一下我的困惑，谢谢！

liyuanhe211 · 发表于 Post on 2017-4-16 11:54:10

clvn 发表于 2017-4-16 11:15
老师，我看了您最后部分的通过做柔性扫描找过渡态，我有几个问题，还想请老师解决一下：
1.我用了您文中的 ...

柔性扫描根本不是opt=TS，看教程
另外最好看前半部分该如何提问。

clvn · 发表于 Post on 2017-4-16 14:37:48

liyuanhe211 发表于 2017-4-16 11:54
柔性扫描根本不是opt=TS，看教程
另外最好看前半部分该如何提问。

老师，您好，通过您的文章，我知道柔性扫描，不是opt=ts,我现在就是通过柔性扫描找到了势垒点，然后通过opt=(calcfc,ts,noeigentest) freq去找过渡态，
但是我仔细地阅读了您的文章，还是不太懂您的“4步优化就找到了正确的过渡态，过渡态 C-H 键键长 1.675 A”以及不用太在意L9999报错说结构不收敛，此时打开输出文件可见路径上的最高点已经找到了”
还请老师指点一下，谢谢！

liyuanhe211 · 发表于 Post on 2017-4-16 17:16:56

clvn 发表于 2017-4-16 14:37
老师，您好，通过您的文章，我知道柔性扫描，不是opt=ts,我现在就是通过柔性扫描找到了势垒点，然后通过o ...

文中说不必理会L9999不是opt=ts时不必理会，而是在opt=modredundant阶段时、若已找到最高点，则不必考虑。

如果还不明白，压缩后上传你的柔性扫描及opt=ts两个输出文件。

clvn · 发表于 Post on 2017-4-18 16:57:33

liyuanhe211 发表于 2017-4-16 17:16
文中说不必理会L9999不是opt=ts时不必理会，而是在opt=modredundant阶段时、若已找到最高点，则不必考虑 ...

老师，我计算的是C6H14=C2H5+C4H9,通过做柔性扫描，发现，在2A左右是势垒，但是势垒附近变化平缓，然后我按照您文章后的操作，进行过渡态搜索，结果就报错了，L9999的错误，之前试过Sob老师博文里的消多于虚频的方法，结果不是很理想，还请老师百忙之中指点一下，谢谢！

liyuanhe211 · 发表于 Post on 2017-4-18 17:27:13

clvn 发表于 2017-4-18 16:57
老师，我计算的是C6H14=C2H5+C4H9,通过做柔性扫描，发现，在2A左右是势垒，但是势垒附近变化平缓，然后我 ...

你没有提供柔性扫描的文件，不知做的是否恰当，但2A这个距离经验上是比较合适的。这种过程的势垒确实一般不大，也常没有过渡态
另外既然在2A左右为扫描曲线的最高点，你为什么要将过渡态初猜中的键长定为3.5A？
L9999错误和“消虚频”没关系
你的优化振荡了，用本帖解决L9999振荡的方法解决。

clvn · 发表于 Post on 2017-4-18 17:38:36

liyuanhe211 发表于 2017-4-18 17:27
你没有提供柔性扫描的文件，不知做的是否恰当，但2A这个距离经验上是比较合适的。这种过程的势垒确实一般 ...

抱歉，老师，是我自己表述不当，我的意思是键长增加了2A的位置，到达3.5A埃才出现势垒高度
老师觉得这个是合理的吗
附上的表格是自己整理的数据，以及柔性扫描文件
震荡和不收敛我都出现了问题，不过我试过了那些方法，似乎对我的这个不太管用……请老师指点一下，谢谢！

liyuanhe211 · 发表于 Post on 2017-4-18 18:14:08

clvn 发表于 2017-4-18 17:38
抱歉，老师，是我自己表述不当，我的意思是键长增加了2A的位置，到达3.5A埃才出现势垒高度
老师觉得这个 ...

你的CCSD(T) “单点扫描曲线”上确实有最大值，但在你找过渡态的 wb97xd/cc-pvtz 水平上却没有最大值。一般的TS都寻找的是电子能势能面上的鞍点，故而你在右侧的、wb97xd/cc-pvtz 势能面上做 opt=TS 是显然找不到过渡态的。右侧图应该下的结论是 wb97xd/cc-pvtz 水平下这个过程很可能没有过渡态。

这个例子倒是提醒了DFT和CCSD(T)的区别。看扫描曲线，或许 CCSD(T) 水平下应该能找到过渡态，但CCSD(T)没有解析梯度，过于耗时，如果算的动的话或许可以尝试用CCSD扫描一下看有没有极大值点，如果有的话，用CCSD找过渡态还算可行、但也不简单。双杂化泛函诸如DSD-PBEP86也可以尝试。

这个过程的过渡态即使找到，也肯定是势垒很低的过渡态，其实意义不大

clvn · 发表于 Post on 2017-4-18 18:35:37

liyuanhe211 发表于 2017-4-18 18:14
你的CCSD(T) “单点扫描曲线”上确实有最大值，但在你找过渡态的 wb97xd/cc-pvtz 水平上却没有最大值。一 ...

听了老师的讲解，我还有两个疑惑：
1.我的这类断裂C-C键或者C-H键，用我自己的方法找出的势垒，势垒点的结构能够代表过渡态的结构吗？和搜索过渡态的相比，会有很大差别吗?
2.之前尝试过CCSD(T),但是代价昂贵，所以放弃了，您说的双杂化泛函诸如DSD-PBEP86，我们没有试过，现在就是按照我们之前的扫描，得到的每一帧做freq，虚频多，而且我发现我扫描其他反应物断裂C-C键，也有同样的情况，但我做C-H键断裂，就没怎么出现这样的情况，这是C-C断裂的共性吗？而且我试过Sob老师关于消虚频的文章里方法，但我的虚频还是有，请问我应该怎么办呢？

liyuanhe211 · 发表于 Post on 2017-4-18 19:48:49

clvn 发表于 2017-4-18 18:35
听了老师的讲解，我还有两个疑惑：
1.我的这类断裂C-C键或者C-H键，用我自己的方法找出的势垒，势垒点的 ...

1. “你自己的方法”如果指柔性扫描，然后计算曲线上各几何结构的高精度单点能，根本就不能代表反应的势垒；可能会出现很大区别。
2. 我没动懂你到底想问什么。这和“消虚频”没关系。

tacoiong · 发表于 Post on 2017-4-20 07:45:19

对于L123，Max corrector steps exceded
可以用recorrect=never, 跳过校正步吗?
因为我用了LQA后曲线不平滑...

话说我最近算出来的过渡态做IRC时都出现L123...

clvn · 发表于 Post on 2017-4-20 17:35:52

liyuanhe211 发表于 2017-4-18 19:48
1. “你自己的方法”如果指柔性扫描，然后计算曲线上各几何结构的高精度单点能，根本就不能代表反应的势 ...

谢谢老师指点，我的疑惑就是断键反应很多都没有明显的势垒，即按我自己的拙见就是没有过渡态，那请问做搜索过渡态有何意义呢？
还请老师解答一下，谢谢！

laomalao · 发表于 Post on 2017-4-21 17:44:24

学习一下

EdRhodes · 发表于 Post on 2017-4-22 18:26:36

各位老师们你们好，在这里我想请教一个关于轨迹运算的问题。之前使用Gaussian默认的ADMP条件admp=(maxpoints=50,stepsize=1000)在不改变基组等条件的基础下可以顺利结束运算。如今想改变为admp=(maxpoints=120,stepsize=10000)，计算总是出错。还请各位老师指点一下。

输入条件：# b3lyp/6-31g(d,p) geom=(connectivity,crowd) admp=(maxpoints=120,stepsize=10000)，以下上传结束前的最后一页输出信息。

liyuanhe211 · 发表于 Post on 2017-4-22 18:31:56

EdRhodes 发表于 2017-4-22 18:26
各位老师们你们好，在这里我想请教一个关于轨迹运算的问题。之前使用Gaussian默认的ADMP条件admp=(maxpoint ...

你看了帖子了么

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

[Gaussian/gview] Gaussian 常见报错及解决方法；新手求助报错时的注意事项