24小时热门版块排行榜    

查看: 1833  |  回复: 13

yangxinyeah

铜虫 (小有名气)

[交流] 【求助】含有赝式基组时应该如何用geom=chenck命令已有6人参与

我想问问各位高手,如果我的输入文件中有赝式基组,当用geom=check命令时应该怎样用呢?之前试过直接在命令行加入geom=allcheck guess=read,结果总报错。求解啊
回复此楼

» 猜你喜欢

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

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

manson1998

金虫 (著名写手)

OLED达人

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流!! 2010-10-26 14:29:35
你想要做什么?是接着算么?还是要力常数?
惺惺相惜情不自禁
2楼2010-10-24 23:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxinyeah

铜虫 (小有名气)

我就是用低基组优化完成后接着用高基组优化。。谢谢解答啊
3楼2010-10-26 08:49:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+1): 2010-10-26 09:12:34
你把低基组优化完的结构重新做个输入文件不就可以了?
惺惺相惜情不自禁
4楼2010-10-26 08:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by yangxinyeah at 2010-10-26 08:49:14:
我就是用低基组优化完成后接着用高基组优化。。谢谢解答啊

你把输入文件发上来才好判断究竟哪里出了问题!
Nothing_Is_Impossible!
5楼2010-10-26 09:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxinyeah

铜虫 (小有名气)

%chk=/home/xy4/Cl-Pd-C20H18N5-1.chk
%mem=3000mb
%nproc=8
# HF/gen opt=(calcfc,z-matrix,maxcyc=265) freq=noraman test
# pseudo=read

opt Pd for SDD, C N H Cl for 3-21G

0 1
C
C                  1              B1
N                  2              B2    1              A1
H                  1              B3    2              A2    3              D1
H                  2              B4    1              A3    3              D2
C                  3              B5    2              A4    1              D3
C                  6              B6    3              A5    2              D4
C                  6              B7    3              A6    2              D5
C                  7              B8    6              A7    3              D6
C                  8              B9    6              A8    3              D7
H                  8             B10    6              A9    3              D8
C                  9             B11    7             A10    6              D9
H                  9             B12    7             A11    6             D10
H                 10             B13    8             A12    6             D11
H                 12             B14    9             A13    7             D12
C                  7             B15    6             A14    3             D13
C                 16             B16    7             A15    6             D14
C                 16             B17    7             A16    6             D15
C                 17             B18   16             A17    7             D16
H                 17             B19   16             A18    7             D17
C                 18             B20   16             A19    7             D18
C                 21             B21   18             A20   16             D19
H                 19             B22   17             A21   16             D20
H                 21             B23   18             A22   16             D21
H                 22             B24   21             A23   18             D22
N                 18             B25   16             A24    7             D23
C                 26             B26   18             A25   16             D24
C                 27             B27   26             A26   18             D25
H                 27             B28   26             A27   18             D26
H                 28             B29   27             A28   26             D27
N                  1             B30    2             A29    3             D28
N                 28             B31   27             A30   26             D29
C                 32             B32   28             A31   27             D30
H                 33             B33   32             A32   28             D31
H                 33             B34   32             A33   28             D32
H                 33             B35   32             A34   28             D33
C                 31             B36    1             A35    2             D34
H                 37             B37   31             A36    1             D35
H                 37             B38   31             A37    1             D36
H                 37             B39   31             A38    1             D37
C                 31             B40    1             A39    2             D38
C                 32             B41   28             A40   27             D39
N                 16             B42    7             A41    6             D40
Pd                43             B43   16             A42    7             D41
Cl                44             B44   41             A43   16             D42

   B1             1.33071875
   B2             1.40184424
   B3             1.06269931
   B4             1.06120074
   B5             1.42604021
   B6             1.39626219
   B7             1.38783285
   B8             1.40242998
   B9             1.37684934
   B10            1.07090926
   B11            1.37485112
   B12            1.07055354
   B13            1.07090245
   B14            1.07212847
   B15            2.41348888
   B16            1.40560282
   B17            1.39371780
   B18            1.37233724
   B19            1.07064165
   B20            1.39150036
   B21            1.37395748
   B22            1.07206210
   B23            1.07144328
   B24            1.07086457
   B25            1.43072598
   B26            1.40764977
   B27            1.32961044
   B28            1.06123871
   B29            1.06291405
   B30            1.39445277
   B31            1.39907848
   B32            1.47159983
   B33            1.07373672
   B34            1.08136131
   B35            1.08245175
   B36            1.47077157
   B37            1.07993506
   B38            1.08126293
   B39            1.07285213
   B40            1.34569962
   B41            1.34494234
   B42            1.39432440
   B43            2.01971801
   B44            2.41000000
   A1           106.87507202
   A2           130.95468935
   A3           130.82529079
   A4           123.21831131
   A5           120.25439671
   A6           118.72082691
   A7           116.89039738
   A8           121.06539129
   A9           119.30784175
   A10          122.01753442
   A11          118.08957006
   A12          120.21016429
   A13          119.62813466
   A14          142.80938386
   A15           96.74816858
   A16          142.56984140
   A17          122.41516199
   A18          117.65903474
   A19          120.52917170
   A20          121.55443188
   A21          119.71397665
   A22          118.68888789
   A23          120.37908780
   A24          122.33878822
   A25          122.93200632
   A26          106.99462454
   A27          122.22740944
   A28          130.82418545
   A29          106.98840052
   A30          107.05903945
   A31          122.40439573
   A32          109.12518068
   A33          108.59269765
   A34          109.65916404
   A35          121.58142627
   A36          109.76546115
   A37          108.24605361
   A38          109.66743277
   A39          111.01560462
   A40          110.82240716
   A41           30.09605245
   A42          118.93663607
   A43           94.93397627
   D1          -179.05006450
   D2           177.70541947
   D3          -175.57854931
   D4          -147.21040132
   D5            29.08675177
   D6           172.80535424
   D7          -173.24556828
   D8             7.40153121
   D9             1.20486157
   D10         -179.72784234
   D11         -179.98068495
   D12          179.72700861
   D13           16.24704240
   D14           86.81214391
   D15         -119.31358172
   D16          163.15798737
   D17          -17.30848017
   D18         -153.80807595
   D19            3.70183075
   D20          179.82731637
   D21         -175.83924069
   D22          179.40801575
   D23           26.60125564
   D24         -135.04774418
   D25          161.98159170
   D26          -17.62884778
   D27          179.12138597
   D28            0.08610791
   D29           -0.65845181
   D30         -168.92813225
   D31         -160.85779039
   D32          -42.05550416
   D33           77.78451146
   D34          175.48160561
   D35          -84.36323193
   D36           36.02240685
   D37          156.35788925
   D38            0.03669956
   D39            0.00280598
   D40          -56.65792953
   D41         -179.69140968
   D42         -159.87127895

C H N Cl 0
3-21G
****
Pd 0
SDD
****

Pd 0
SDD



当用3-21G优化完成后,我把第一步产生的Cl-Pd-C20H18N5-1.chk复制为Cl-Pd-C20H18N5-2.chk,然后用以下的输入文件接着用6-31g优化
%chk=/home/xy4/Cl-Pd-C20H18N5-2.chk
%mem=3000mb
%nproc=8
# HF/gen geom=allcheck guess=read  freq=noraman test
# pseudo=read opt=(calcfc,z-matrix,maxcyc=265)

opt Pd for SDD, C N H Cl for 3-21G

但是总报错。我看别人用在没有赝式基组时用 geom=allcheck guess=read命令都可以算的。难道是由于我有赝式基组??
6楼2010-10-26 09:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流!! 2010-10-26 14:29:18
基组不同 波函数不同,不能用guess=read吧
惺惺相惜情不自禁
7楼2010-10-26 12:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流!! 2010-10-26 14:29:25
第二步里为何用calcfc,而不用rcfc?
8楼2010-10-26 14:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxinyeah

铜虫 (小有名气)

引用回帖:
Originally posted by manson1998 at 2010-10-26 12:13:41:
基组不同 波函数不同,不能用guess=read吧

貌似不是这个原因哦~如果没有赝式,比如对所有原子都用的是HF/3-21G。当接着用HF/6-31G优化时,用guess=read就可以的。
9楼2010-10-26 15:24:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxinyeah

铜虫 (小有名气)

引用回帖:
Originally posted by fatpig8832 at 2010-10-26 14:26:57:
第二步里为何用calcfc,而不用rcfc?

不好意思,之前写错了。第二步就是用的rcfc.命令如下
# B3LYP/gen geom=allcheck guess=read
# opt(rcfc,z-matrix,maxcyc=265) freq=noraman test
# pseudo=read
10楼2010-10-26 15:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangxinyeah 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见