计算化学公社

标题: 过渡金属镍体系跑IRC出现难解决的问题（已经仔细阅读过sob老师的IRC帖子） [打印本页]

作者
Author: Wangying001 时间: 2022-9-26 17:18
标题: 过渡金属镍体系跑IRC出现难解决的问题（已经仔细阅读过sob老师的IRC帖子）
本帖最后由 Wangying001 于 2022-9-27 14:06 编辑

如题，我在寻找过渡态需要跑IRC验证是否找对
（1）第一次使用gaussian默认关键词报错（IRC=calcfc)，发现校正步不收敛。
（2）后换用了GS2算法IRC=(GS2,calcfc)后还是报错校正步不收敛，此时我只能采用不太精准的IRC=(LQA,calcfc)，此时报错信息变了，提示我SCF不收敛。
（3）后我增加关键词SCF=vshift=300 ，可以得到IRC曲线，但由于该曲线不光滑不完整，我想得到一条完整、高质量的IRC曲线。故增加关键词IRC=(calcfc,LQA,maxpoint=80,stepsize=10) SCF=vshift=300
（4）此时gaussian又继续报错，报错又变成了原来的SCF不收敛。信息为Convergence failure -- run terminated.
Error termination via Lnk1e in /public1/home/scb0643/g16/l502.exe at Mon Sep 26 16:05:48 2022.

我的问题是：
（1）如何在不希望用LQA避免报错的同时牺牲IRC精度的情况下得到一条完整的IRC曲线？关键词还能怎么增加？（我已经用过GS2算法会报错，只有LQA不会）
（2）想搞清楚使用LQA为什么增加步数以后会报错？

作者
Author: Wangying001 时间: 2022-9-26 17:24
我需要计算的体系是含过渡金属镍的
完整关键词是：
%mem=100GB
%nprocshared=64
#p IRC=(calcfc,LQA,maxpoint=60,stepsize=10) b3lyp/genecp 5d empiricaldispersion=gd3 SCF=vshift=300

Title Card Required

0 2
C                -5.09233700 -0.80210900 -0.44763100
C                -3.58241300 -0.46751700 1.25206000
C                -4.35349600 0.60526400 1.80184700
C                -5.54699100 0.96240600 1.12118000
C                -5.91142800 0.27215600 -0.01041000
C                -2.33695500 -0.80554300 1.86997100
C                -3.90108600 1.26169700 2.97581400
H                -6.15343100 1.77982800 1.50281800
H                -6.81090100 0.52486400 -0.56223500
C                -2.72388500 0.87706700 3.58088600
C                -1.93533000 -0.15378500 3.02161200
H                -4.49647700 2.07271200 3.38730900
H                -2.38084100 1.38157300 4.47872900
H                -0.97953400 -0.40860100 3.46408000
C                -1.46262900 -1.80751600 1.25025500
O                -1.02434600 -2.83243800 2.00762400
C                -0.18723200 -2.97858800 -0.04574900
C                -0.21446200 -3.57086000 1.17707700
H                0.22862400 -4.45635000 1.59837000
C                0.45696100 -3.47733400 -1.32589700
C                -0.66559800 -4.03038800 -2.23854000
C                1.18515400 -2.35972400 -2.01615000
C                1.45551900 -4.60565100 -0.98341300
H                -1.18892800 -4.85505800 -1.74278100
H                -1.38371700 -3.24396800 -2.47923200
H                -0.23084000 -4.40328600 -3.17194400
H                2.18625800 -2.09872100 -1.69006600
H                0.67102200 -1.74998100 -2.74976000
H                1.92522000 -4.97018100 -1.90098500
H                2.24470000 -4.24849800 -0.31330500
H                0.95050600 -5.45394200 -0.50817300
Ni             -1.25025100 -0.00684300 -0.98388200
N                -3.98009700 -1.16830200 0.15434300
N                -0.97948500 -1.82835600 0.03399700
C                5.39516800 -0.87239500 -0.47314800
C                4.28781900 -0.38867600 -1.15608800
C                3.33388200 0.42088500 -0.51717300
C                3.52830100 0.73211500 0.82720400
C                4.63892200 0.25356300 1.53051300
C                5.57823900 -0.55353600 0.88110700
H                6.13338100 -1.49933700 -0.96228800
H                4.16019100 -0.64212900 -2.20634400
H                2.79036600 1.33588700 1.34067800
H                4.75099700 0.51118800 2.57693500
O                6.69514600 -1.08076000 1.46637000
C                6.92759400 -0.79807300 2.83636300
H                6.11614700 -1.17566200 3.47334300
H                7.85554600 -1.30947600 3.09822300
H                7.04766400 0.27913700 3.01542400
C                2.16262200 0.93302500 -1.33029500
H                2.52773300 1.62964800 -2.09284800
H                1.69560000 0.10640000 -1.87682200
N                1.12807600 1.65871000 -0.57934700
C                0.21164600 0.90547600 0.17243700
O                0.49987900 0.28535400 1.17242300
Cl             -1.87324500 -0.68518800 -3.00505500
C                0.55802400 2.85905000 -1.19231800
C                -0.59787000 3.26925900 -0.24361800
H                0.11398400 2.61583300 -2.16982300
H                -1.19900300 4.01440300 -0.77683600
C                -1.46479500 2.01630300 -0.05552500
C                -2.25355300 1.61333600 -1.18964000
H                -1.90101300 1.88615700 0.93207300
H                -3.27231500 1.26761300 -1.02936800
H                -2.10927900 2.12862600 -2.13673700
C                1.59447900 3.95327800 -1.41382300
C                1.43542000 4.78888100 -2.44596200
H                2.11497800 5.61856700 -2.61701600
H                0.61806900 4.67412700 -3.15344200
C                2.76115800 4.07017900 -0.46214300
H                3.32242600 4.98913100 -0.65163300
H                2.44435000 4.06556800 0.58342500
H                3.44841000 3.22489300 -0.56953600
H                -5.36859600 -1.37158400 -1.33238600
C                -0.17071200 3.86799800 1.10080300
H                -1.05669400 4.06536200 1.71290500
H                0.47724500 3.18646400 1.65958600
H                0.35465400 4.81552000 0.96035000

H C N O Cl 0
6-31g(d,p)
****
Ni 0
LanL2dz
****

Ni 0
LanL2dz

作者
Author: chands 时间: 2022-9-26 17:45
是Ni(II)吗，如果是Ni(II)，总觉得电荷数和自旋多重度不配。

作者
Author: Wangying001 时间: 2022-9-26 17:50

chands 发表于 2022-9-26 17:45
是Ni(II)吗，如果是Ni(II)，总觉得电荷数和自旋多重度不配。

对的，是NI(Ⅱ)体系，但是和报错信息应该无关。电荷数和自旋多重度我是gussian默认的。

作者
Author: chands 时间: 2022-9-26 18:04

Wangying001 发表于 2022-9-26 17:50
对的，是NI(Ⅱ)体系，但是和报错信息应该无关。电荷数和自旋多重度我是gussian默认的。

含羰基的配体-1价，氯离子-1价，含恶唑的配体0价，Ni(II)，是这样吗？如果是这样，数一数Ni(II)的3d轨道的电子。

作者
Author: fatpig 时间: 2022-9-26 18:34
既然过渡态已经算出来了，为何还要calcfc？直接readfc即可...
另，直接guess=read geom=check即可...

作者
Author: Wangying001 时间: 2022-9-26 18:51

fatpig 发表于 2022-9-26 18:34
既然过渡态已经算出来了，为何还要calcfc？直接readfc即可...
另，直接guess=read geom=check即可...

直接读取过渡态产生的波函数看是否会报错吗？

作者
Author: Wangying001 时间: 2022-9-26 19:14

chands 发表于 2022-9-26 18:04
含羰基的配体-1价，氯离子-1价，含恶唑的配体0价，Ni(II)，是这样吗？如果是这样，数一数Ni(II)的3d轨道 ...

我按照您说的重新设置电荷与自旋多重度为0 1以后报错
The combination of multiplicity 1 and 315 electrons is impossible.
Error termination via Lnk1e in /public1/home/scb0643/g16/l301.exe at Mon Sep 26 18:30:18 2022.

作者
Author: chands 时间: 2022-9-26 19:42
本帖最后由 chands 于 2022-9-26 19:43 编辑

Wangying001 发表于 2022-9-26 19:14
我按照您说的重新设置电荷与自旋多重度为0 1以后报错
The combination of multiplicity 1 and 315 el ...

检查一下那个含恶唑的配体，好像端位甲基缺了个H

作者
Author: sobereva 时间: 2022-9-27 02:41
首先注意前面的人说的电荷、自旋多重度问题，这个要错了其它都是白说。

记住Gaussian怎么拼

减小IRC步长、使用recalc=5，有助于改善IRC质量，以及增加令SCF能够收敛的状态延续下去的概率，此文明确说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

写stepsize=10毫无意义，本来就是默认的

写5d也完全多余，用gen、genecp的时候默认就是5d

作者
Author: fatpig 时间: 2022-9-27 08:50

Wangying001 发表于 2022-9-26 18:51
直接读取过渡态产生的波函数看是否会报错吗？

不是看报错，而是一定要读取吧...

作者
Author: Wangying001 时间: 2022-9-27 14:05

sobereva 发表于 2022-9-27 02:41
首先注意前面的人说的电荷、自旋多重度问题，这个要错了其它都是白说。

记住Gaussian怎么拼

谢谢sob老师的耐心解答！

作者
Author: sai77 时间: 2022-9-27 18:02
配体上的叔丁基少了一个氢，结构错了，前面全白算= =
镍最外层10个电子，二价镍8个电子，都是偶数，没有未成对电子，中性体系，自旋多重度也不可能是2啊

作者
Author: Wangying001 时间: 2022-9-27 19:11

sai77 发表于 2022-9-27 18:02
配体上的叔丁基少了一个氢，结构错了，前面全白算= =
镍最外层10个电子，二价镍8个电子，都是偶数，没有未 ...

是的，已经发现结构有问题了！谢谢您的解答！

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