24小时热门版块排行榜    

查看: 827  |  回复: 1

gscn

新虫 (小有名气)

[求助] QTS2运行时停留在L103 已有1人参与

想用QTS2结构优化,然后IRC算过渡态,但是在结构优化时总是在L103的时候卡住,显示Performing Berny Optimization。
%mem=1000MB
%chk=C:\Users\ASUS\Desktop\HfC\HfCl-N\Hf-NH2 reaction.chk
# opt=(calcfc,qst2,cartesian) freq b3lyp/lanl2dz geom=connectivity

Title Card Required

0 1
C                  1.42078200   -0.00001400    0.00876900
C                  0.69518200   -1.19720800    0.00553600
C                 -0.69519600   -1.19719500   -0.00552800
C                 -1.42075600    0.00001100   -0.00889400
C                 -0.69516600    1.19718400   -0.00569600
C                  0.69520800    1.19718900    0.00562300
H                  1.22558100   -2.14606300    0.01767200
H                 -1.22555500   -2.14606600   -0.01760600
H                 -1.22556200    2.14603500   -0.01773600
H                  1.22556100    2.14606400    0.01783200
N                 -2.82729100    0.00004500   -0.09031600
H                 -3.25532100   -0.82826500    0.29994800
H                 -3.25540100    0.82811100    0.30038900
N                  2.82721500   -0.00000500    0.09035400
H                  3.25541700   -0.82817800   -0.29997100
H                  3.25548800    0.82828400   -0.29965400
Cl                 4.29265444    3.14458792   -1.73386128
Hf                 4.51337999    0.85173409   -2.15423011
Cl                 4.23688213   -0.33463575   -0.15455915
Cl                 2.88405729    0.17838085   -3.69517506
Cl                 6.63992612    0.41860335   -3.03332495

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

Title Card Required

0 1
C                  1.42078200   -0.00001400    0.00876900
C                  0.69518200   -1.19720800    0.00553600
C                 -0.69519600   -1.19719500   -0.00552800
C                 -1.42075600    0.00001100   -0.00889400
C                 -0.69516600    1.19718400   -0.00569600
C                  0.69520800    1.19718900    0.00562300
H                  1.22558100   -2.14606300    0.01767200
H                 -1.22555500   -2.14606600   -0.01760600
H                 -1.22556200    2.14603500   -0.01773600
H                  1.22556100    2.14606400    0.01783200
N                 -2.82729100    0.00004500   -0.09031600
H                 -3.25532100   -0.82826500    0.29994800
H                 -3.25540100    0.82811100    0.30038900
N                  2.82721500   -0.00000500    0.09035400
H                  3.25541700   -0.82817800   -0.29997100
H                  3.25548800    0.82828400   -0.29965400
Cl                 4.29265444    3.14458792   -1.73386128
Hf                 4.51337999    0.85173409   -2.15423011
Cl                 4.23688213   -0.33463575   -0.15455915
Cl                 2.88405729    0.17838085   -3.69517506
Cl                 6.63992612    0.41860335   -3.03332495

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

C H Cl N 0
6-31 G
****

Hf 0
lanl2dz
****
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+8 2020-09-24 13:19:56
gscn: 金币+10 2020-10-12 11:08:32
最大的问题在自定义基组的写法上
1. 关键词部分 应该写 GENECP而不是 lanl2dz   因为你要读取文件结尾的自定义混合基组。就必须要写GenECP要不然就根本不会去读取那些你写的基组。而是对所有的原子都使用LanL2DZ这显然是很不合理的设置。

2. 自定义基组部分的格式也有问题。有多余的空行,以及没有定义赝势部分。 具体自定义基组如何书写见我这个教程1.1节  https://muchong.com/t-14242735-1

3. 6-31G基组 太差 至少换成 6-31G(d)

4. 不知道你有没有设置默认的核心数目。如果就是用单核计算会很慢。即使是阉割的32bit win版高斯也可以用到4核心。
2楼2020-09-23 22:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gscn 的主题更新
信息提示
请填处理意见