计算化学公社

标题: g09 scan出错提问 [打印本页]

作者
Author: 小菲丽 时间: 2017-11-19 10:10
标题: 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

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

作者
Author: 小菲丽 时间: 2017-11-19 10:38

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

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

作者
Author: Mikasa 时间: 2017-11-19 11:03

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

对

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

作者
Author: 小范范1989 时间: 2017-11-19 15:52

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**等计算对比一下？谢谢李老师指点

作者
Author: liyuanhe211 时间: 2017-11-19 17:20

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

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

作者
Author: 小范范1989 时间: 2017-11-19 19:07

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

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

作者
Author: 小菲丽 时间: 2017-11-21 17:04

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

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

作者
Author: sobereva 时间: 2017-11-25 16:09

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

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

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

作者
Author: wwen 时间: 2020-8-12 19:59

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

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

作者
Author: chrischen1128 时间: 2020-8-12 20:53

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

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

作者
Author: wwen 时间: 2020-8-13 09:14
本帖最后由 wwen 于 2020-8-13 19:24 编辑

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

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

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