g09 scan出错提问

小菲丽

我用g09走IRC reverse，只跑出来6个点，然后我想固定7号碳与10号碳的键长，优化过渡态剩下的结构，出错如下，我的输入应该怎么改一下？
%mem=800MW
%nprocshared=4
Will use up to    4 processors via shared memory.
------------------------------------------------------------
# scan Opt=ModRedundant b3lyp/6-311g(2d,2p) scf=maxcycle=500
------------------------------------------------------------
CompJT:  unrecognized IType=10.
Error termination via Lnk1e in /home/ying/g09/l1.exe at Sun Nov 19 09:58:11 2017.
Job cpu time:       0 days  0 hours  0 minutes  0.1 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      

输入文件如下
%mem=800MW
%nprocshared=4
# scan Opt=ModRedundant b3lyp/6-311g(2d,2p) scf=maxcycle=500

Title Card Required

0 1                                                            
H                 -2.232830   -1.368602    0.854839            
N                 -3.258157   -0.340854   -1.072977            
C                 -4.546973   -0.905615   -0.665618            
C                 -4.585058   -1.244052    0.893363            
O                 -5.715772   -1.351071    1.337085            
H                 -5.319735   -0.158561   -0.821234            
C                 -0.753200   -0.086535    0.455305            
O                 -1.333296   -1.327449    0.338036            
C                 -2.758357    0.805964   -0.719802            
C                 -1.287867    0.927209   -0.776810            
C                 -3.605333    1.836191   -0.041706            
C                 -0.763371    2.372605   -0.735816            
H                 -0.911383    0.420199   -1.665187            
C                 -3.037781    3.258348   -0.129603            
H                 -4.625751    1.793548   -0.420230            
H                 -3.669313    1.522017    1.008564            
C                 -1.549220    3.282837    0.213220            
H                 -0.818227    2.789527   -1.745040            
H                  0.292396    2.361150   -0.472571            
H                 -3.184561    3.648313   -1.140185            
H                 -3.605103    3.904081    0.540349            
H                 -1.407595    2.964012    1.249477            
H                 -1.163287    4.300628    0.146422            
H                 -1.027275    0.402835    1.394153            
C                  0.752544   -0.172147    0.344597            
C                  1.567811    0.688453    1.076265            
C                  1.356137   -1.098047   -0.507812            
C                  2.952525    0.653921    0.948767            
H                  1.123950    1.390385    1.770431            
C                  2.737537   -1.144633   -0.647717            
H                  0.738796   -1.803704   -1.044925            
C                  3.524689   -0.259989    0.078299            
H                  3.572945    1.323235    1.525440            
H                  3.196510   -1.866160   -1.306993            
O                 -3.479238   -1.375027    1.498230            
C                 -4.847602   -2.152615   -1.487708            
H                 -5.797398   -2.561524   -1.153272            
H                 -4.913763   -1.930809   -2.553531            
H                 -4.081092   -2.914295   -1.332371            
Br                 5.433940   -0.314534   -0.110998            
H                 -2.561215   -1.035582   -1.326867            
                                                                                                                           
7 10 1.9 s 5 -0.1

Mikasa

1. 没搞明白你要干嘛，什么叫“优化过渡态剩下的结构”？
2. IRC问题和扫描有什么关系，IRC跑不完看下边PPT
3. scan是刚性扫描，opt=modredundant是柔性扫描，这两个关键词同时用是怎么回事。
4. 你都没说你是什么错误。

小菲丽

Mikasa 发表于 2017-11-19 10:21
1. 没搞明白你要干嘛，什么叫“优化过渡态剩下的结构”？
2. IRC问题和扫描有什么关系，IRC跑不完看下边PP ...

谢谢大神，我想柔性扫描7和10之间的距离，我的输入中不应该加关键词SCAN对吧？

Mikasa

小菲丽 发表于 2017-11-19 10:38
谢谢大神，我想柔性扫描7和10之间的距离，我的输入中不应该加关键词SCAN对吧？

对

liyuanhe211

scan 和 opt 冲突。想好要柔性扫描还是刚性扫描
6-311G(2d,2p)是极不匹配基组，去sobereva.com看量化基组的选择博文
scf=maxcyc=500垃圾关键词

小范范1989

liyuanhe211 发表于 2017-11-19 15:32
scan 和 opt 冲突。想好要柔性扫描还是刚性扫描
6-311G(2d,2p)是极不匹配基组，去sobereva.com看量化基组 ...

借楼主的帖子，咨询一下李老师。
我一个文章，计算都是6-31g*。然后审稿人这么说：The authors only tested and used double-zeta quality basis functions. They should provide a benchmark for basis sets including triple-zeta quality basis sets.
是不是建议要用6-311g*或者6-311g**等计算对比一下？谢谢李老师指点

liyuanhe211

小范范1989 发表于 2017-11-19 15:52
借楼主的帖子，咨询一下李老师。
我一个文章，计算都是6-31g*。然后审稿人这么说：The authors only tes ...

是这个意思。算的动就尽量算一下。

小范范1989

liyuanhe211 发表于 2017-11-19 17:20
是这个意思。算的动就尽量算一下。

恩，好的，谢谢李老师指点。

小菲丽

liyuanhe211 发表于 2017-11-19 15:32
scan 和 opt 冲突。想好要柔性扫描还是刚性扫描
6-311G(2d,2p)是极不匹配基组，去sobereva.com看量化基组 ...

老师，6-311G(2d,2p)是极不匹配基组，这是什么意思？是不能用这个基组计算吗？

sobereva

小菲丽 发表于 2017-11-21 17:04
老师，6-311G(2d,2p)是极不匹配基组，这是什么意思？是不能用这个基组计算吗？

仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336

重原子极化才2d，根本没必要给氢加这么多极化。想要比6-311G**更好的结果，用6-311G(2d,2p)明显不如用TZVP

wwen

Mikasa 发表于 2017-11-19 10:21
1. 没搞明白你要干嘛，什么叫“优化过渡态剩下的结构”？
2. IRC问题和扫描有什么关系，IRC跑不完看下边PP ...

您好，在做柔性扫描时我想固定两个原子之间的键长优化扫描其他原子，在写输入文件时最下面的被扫描变量该如何写？

chrischen1128

wwen 发表于 2020-8-12 19:59
您好，在做柔性扫描时我想固定两个原子之间的键长优化扫描其他原子，在写输入文件时最下面的被扫描变量该 ...

详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

wwen

chrischen1128 发表于 2020-8-12 20:53
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.htm ...

非常感谢您的回答，我明白了