应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】TS过渡态总是L301错误，我要哭了~~~
131/1返回列表
查看: 2283  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lixiaocat

银虫 (小有名气)

[交流] 【求助】TS过渡态总是L301错误，我要哭了~~~

要做一个羟基自由基夺取2,4,5三氯酚上羟基氢的过渡态
这是我的输入文件
%chk=ts1-2
%mem=150MW
%nproc=8
#p opt=(calcfc,ts,noeigentest,maxcycl=200) freq mp2/6-311+g(d,p)
scrf=(solvent=water) guess=read geom=connectivity optcyc=200
iop(1/11=1,2/16=1,5/13=1)

ts1-2

0 2
C
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
C                  4              B4    3              A3    2              D2
C                  5              B5    4              A4    3              D3
H                  2              B6    1              A5    6              D4
H                  5              B7    4              A6    3              D5
O                  4              B8    3              A7    2              D6
H                  9              B9    4              A8    3              D7
Cl                 3             B10    2              A9    1              D8
Cl                 1             B11    2             A10    3              D9
Cl                 6             B12    5             A11    4             D10
O                  9             B13    4             A12    3             D11
H                 14             B14    9             A13    4             D12

   B1             1.38960122
   B2             1.38958796
   B3             1.40463410
   B4             1.39862108
   B5             1.38770592
   B6             1.08147126
   B7             1.08178453
   B8             1.34262550
   B9             1.05400000
   B10            1.72823000
   B11            1.74570346
   B12            1.74434560
   B13            2.35239415
   B14            0.97370785
   A1           120.15822454
   A2           121.54171910
   A3           117.47551141
   A4           121.38215925
   A5           119.79335544
   A6           118.22078710
   A7           124.82662770
   A8           128.00000000
   A9           118.66745836
   A10          118.87027986
   A11          118.52464701
   A12          121.31114852
   A13          100.52534923
   D1             1.57361244
   D2            -1.60846309
   D3             0.48428673
   D4          -179.22255337
   D5          -179.93166928
   D6           178.77857492
   D7           -18.95612259
   D8          -176.86511936
   D9           179.38358497
   D10         -179.87723870
   D11          -31.42153557
   D12          -52.47926078

1 2 1.5 6 1.5 12 1.0
2 3 1.5 7 1.0
3 4 1.5 11 1.0
4 5 1.5 9 1.5
5 6 1.5 8 1.0
6 13 1.0
7
8
9 10 1.0
10 14 1.0
11
12
13
14 15 1.0
15





报错是
United Atom Topological Model (UA0  parameters set).
UA0: Hydrogen   10 has 2 bonds. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.
Error termination via Lnk1e in /public/soft/Gauss03/g03/l301.exe at Thu Apr 14 10:33:26 2011.
Job cpu time:  0 days  0 hours  0 minutes  8.2 seconds.



坐等救命
回复此楼

» 收录本帖的淘帖专辑推荐

lywiailyw

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2011-04-14 10:53:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)
★
lixiaocat(金币+1):谢谢参与
引用回帖:
Originally posted by lixiaocat at 2011-04-14 10:53:03:
要做一个羟基自由基夺取2,4,5三氯酚上羟基氢的过渡态
这是我的输入文件
%chk=ts1-2
%mem=150MW
%nproc=8
#p opt=(calcfc,ts,noeigentest,maxcycl=200) freq mp2/6-311+g(d,p)
scrf=(solvent=water) guess= ...

为什么不用用qts3或qts2啊
赞一下(1人) 回复此楼
2楼2011-04-14 12:44:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)
★ ★ ★
lixiaocat(金币+1):谢谢参与
375642546(金币+2): 2011-04-14 16:27:37
楼主你...

Hydrogen   10 has 2 bonds.

这么大的字你莫非没看到...

简单一点，把 geom=connectivity 删去

再把
1 2 1.5 6 1.5 12 1.0
2 3 1.5 7 1.0
3 4 1.5 11 1.0
4 5 1.5 9 1.5
5 6 1.5 8 1.0
6 13 1.0
7
8
9 10 1.0
10 14 1.0
11
12
13
14 15 1.0
15

全部删去

再试试...
赞一下(2人) 回复此楼
3楼2011-04-14 16:22:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangbohejin

木虫之王 (文学泰斗)
★
lixiaocat(金币+1):谢谢参与
用qts2--
回复此楼
4楼2011-04-14 19:26:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjww0001

金虫 (小有名气)
★
lixiaocat(金币+1):谢谢参与
不要用gv直接创建输入文件，那是垃圾，没有用的。
赞一下(1人) 回复此楼
5楼2011-04-14 19:52:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

曼联无敌

木虫 (著名写手)
★ ★
lixiaocat(金币+1):谢谢参与
375642546(金币+1): 2011-04-15 08:51:45
L301是生成基组信息的步骤，应该是命令行的问题，报错也说了H10有两个键，显然不对，建议楼主再检查一下体系的结构。
赞一下(2人) 回复此楼
6楼2011-04-15 08:16:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jin912

木虫 (著名写手)
★ ★
lixiaocat(金币+1):谢谢参与
375642546(金币+1): 2011-04-15 08:51:29
Hydrogen   10 has 2 bonds
你的初始结构都不正确
还有就是建议用qst2，你用ts，而且还有溶剂化~~~
赞一下(2人) 回复此楼
7楼2011-04-15 08:46:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhongrl898

木虫 (小有名气)
★
lixiaocat(金币+1):谢谢参与
构型猜测不正确，建议使用QST2
赞一下(1人) 回复此楼
8楼2011-04-15 09:13:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice625

铜虫 (正式写手)
★
lixiaocat(金币+1):谢谢参与
把
1 2 1.5 6 1.5 12 1.0
2 3 1.5 7 1.0
3 4 1.5 11 1.0
4 5 1.5 9 1.5
5 6 1.5 8 1.0
6 13 1.0
7
8
9 10 1.0
10 14 1.0
11
12
13
14 15 1.0
15
删去，建议做好用CHEMCRAFT软件画图，直接拷坐标替换固定的程序文件就行了
记得坐标下面还要空上几行
赞一下(1人) 回复此楼
9楼2011-04-28 00:31:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangbaowei

木虫 (正式写手)
★
lixiaocat(金币+1):谢谢参与
估计建模有问题
回复此楼
10楼2011-05-03 20:33:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)
★
lixiaocat(金币+1):谢谢参与
应该不是结构的问题，而是溶剂化模型的问题，G03中的PCM默认的radii是UA0，而G09中PCM的radii默认是UFF，这两种radii在处理H的时候有一些区别……
赞一下(1人) 回复此楼
11楼2011-10-12 11:04:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
scrf=(cpcm,solvent=water,read)
......


radii=uff
赞一下(1人) 回复此楼
12楼2011-10-13 18:33:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zszjn

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
要做一个结构优化（考虑溶剂效应）
这是输入文件：
# B3LYP/6-31G** geom=connectivity OPT SCRF=(PCM,READ,Solvent=water) test SCF=Tight

A44

0 1
C
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
N                  1              B4    2              A3    3              D2
H                  2              B5    1              A4    5              D3
H                  3              B6    2              A5    1              D4
H                  4              B7    3              A6    2              D5
C                  1              B8    5              A7    4              D6
C                  9              B9    1              A8    5              D7
C                  5             B10    1              A9    2              D8
H                 11             B11    5             A10    1              D9
H                 11             B12    5             A11    1             D10
H                 11             B13    5             A12    1             D11
C                  9             B14    1             A13    5             D12
H                 15             B15    9             A14    1             D13
C                  9             B16    1             A15    5             D14
N                 17             B17    9             A16    1             D15
O                 18             B18   17             A17    9             D16
O                 10             B19    9             A18    1             D17
C                 17             B20    9             A19    1             D18
H                 21             B21   17             A20    9             D19
H                 21             B22   17             A21    9             D20
H                 21             B23   17             A22    9             D21

   B1             1.40175000
   B2             1.40124782
   B3             1.40184834
   B4             1.39179783
   B5             1.08481903
   B6             1.08481145
   B7             1.08923173
   B8             2.10384340
   B9             1.20662307
   B10            1.47000000
   B11            1.07000000
   B12            1.07000000
   B13            1.07000000
   B14            1.36170000
   B15            1.07000000
   B16            1.38029460
   B17            1.28500000
   B18            1.32587874
   B19            1.43000000
   B20            1.49230000
   B21            1.07000000
   B22            1.07000000
   B23            1.07000000
   A1           109.07834178
   A2           106.74085740
   A3           106.50573988
   A4           126.44075323
   A5           126.81725494
   A6           130.11211410
   A7           155.50317747
   A8           109.11082165
   A9           125.59658900
   A10          109.47120255
   A11          109.47120255
   A12          109.47123134
   A13           41.24974249
   A14          121.82032679
   A15          152.23676038
   A16          111.83860613
   A17          100.64448547
   A18          124.55736225
   A19          133.72548892
   A20          109.47120255
   A21          109.47120255
   A22          109.47123134
   D1            -0.00472875
   D2             0.00000000
   D3           180.00000000
   D4          -180.00000000
   D5          -179.99201401
   D6           147.76327266
   D7          -150.12552641
   D8           179.98442604
   D9           150.38072621
   D10          -89.61925899
   D11           30.38073361
   D12           45.68019108
   D13          179.07632216
   D14           35.11549995
   D15          156.44810647
   D16           16.38370461
   D17           10.26511811
   D18           -8.15428910
   D19          -48.61759650
   D20           71.38241829
   D21         -168.61758910

1 2 2.0 5 1.0 15 1.0
2 3 1.0 6 1.0
3 4 2.0 7 1.0
4 5 1.0 8 1.0
5 11 1.0
6
7
8
9 10 1.0 15 2.0 17 1.0
10 19 1.0 20 2.0
11 12 1.0 13 1.0 14 1.0
12
13
14
15 16 1.0
16
17 18 2.0 21 1.0
18 19 1.0
19
20
21 22 1.0 23 1.0 24 1.0
22
23
24

ALPHA=1.21
TSARE=0.4
报错信息：
United Atom Topological Model (UA0  parameters set).
UA0: Hydrogen    6 has 2 bounds. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.
Error termination via Lnk1e in d:\gauss\l301.exe at Tue May 06 09:50:02 2014.
Job cpu time:  0 days  0 hours  0 minutes  2.0 seconds.
File lengths (MBytes):  RWF=     12 Int=      0 D2E=      0 Chk=      1 Scr=      1
声明一下，在不考虑溶剂的情况下没有错误。
赞一下 回复此楼
13楼2014-05-06 10:00:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lixiaocat 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭