24小时热门版块排行榜    

查看: 467  |  回复: 5
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 sxllgaoli898 的 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sxllgaoli898

金虫 (小有名气)

[交流] 【求助】gaussian计算问题

计算pd(ch3ch2)4分子,我的输入文件为:
%mem=1000mb
%rwf=loc1,1900mb,loc2,1900mb,loc3,-1
#p B3LYP/Gen pseudo=cards opt freq geom=connectivity

Title Card Required

0 1
C
H                  1              B1
H                  1              B2    2              A1
H                  1              B3    3              A2    2              D1
C                  1              B4    3              A3    4              D2
H                  5              B5    1              A4    3              D3
H                  5              B6    1              A5    3              D4
C                  5              B7    1              A6    3              D5
H                  8              B8    5              A7    1              D6
H                  8              B9    5              A8    1              D7
C                  8             B10    5              A9    1              D8
H                 11             B11    8             A10    5              D9
H                 11             B12    8             A11    5             D10
H                 11             B13    8             A12    5             D11
C                  8             B14    5             A13    1             D12
H                 15             B15    8             A14    5             D13
H                 15             B16    8             A15    5             D14
C                 15             B17    8             A16    5             D15
H                 18             B18   15             A17    8             D16
H                 18             B19   15             A18    8             D17
H                 18             B20   15             A19    8             D18
C                 15             B21    8             A20    5             D19
H                 22             B22   15             A21    8             D20
H                 22             B23   15             A22    8             D21
C                 22             B24   15             A23    8             D22
H                 25             B25   22             A24   15             D23
H                 25             B26   22             A25   15             D24
H                 25             B27   22             A26   15             D25
Pb                 5             B28    1             A27   22             D26

   B1             1.09475311
   B2             1.09475311
   B3             1.09507911
   B4             1.53378455
   B5             1.09718801
   B6             1.09718801
   B7             2.49412400
   B8             1.09718801
   B9             1.09718801
   B10            1.53378455
   B11            1.09475311
   B12            1.09507911
   B13            1.09475311
   B14            2.56177515
   B15            1.09718801
   B16            1.09718801
   B17            1.53378455
   B18            1.09475311
   B19            1.09507911
   B20            1.09475311
   B21            2.49412400
   B22            1.09718801
   B23            1.09718801
   B24            1.53378455
   B25            1.09507911
   B26            1.09475311
   B27            1.09475311
   B28            1.55507389
   A1           107.57318734
   A2           107.22255404
   A3           112.17286837
   A4           108.63972299
   A5           108.63972299
   A6           153.79712127
   A7            86.66079806
   A8            86.66079806
   A9           153.79712127
   A10          112.17286837
   A11          110.22030590
   A12          112.17286837
   A13           60.86978877
   A14           91.39979900
   A15          142.86225141
   A16           96.69098057
   A17          112.17286837
   A18          110.22030590
   A19          112.17286837
   A20           60.86978877
   A21           86.66079806
   A22           86.66079806
   A23          153.79712127
   A24          110.22030590
   A25          112.17286837
   A26          112.17286837
   A27          117.11202077
   D1          -115.04502627
   D2          -121.13676266
   D3            62.50219896
   D4           176.28896073
   D5           -60.60442015
   D6          -127.33832288
   D7           127.33832288
   D8             0.00000000
   D9            60.60442015
   D10          180.00000000
   D11          -60.60442015
   D12          -33.86815478
   D13          157.29981422
   D14           39.71831246
   D15          -93.74849125
   D16           89.95383212
   D17         -150.65058803
   D18          -31.25500818
   D19           71.90670940
   D20         -161.20647766
   D21           93.47016809
   D22          -33.86815478
   D23          180.00000000
   D24          -60.60442015
   D25           60.60442015
   D26          -29.34941197

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

  C 0
  6-31G
  ****
  H 0
  6-31G
  ****
  pb 0
  lanl2dz
  ****

  pb 0
  lanl2dz

当考虑其对称性时,计算结果有4个虚频,偶极矩为0,可不考虑对称性无虚频,但其偶极矩不为0,而理论上它的偶极矩应该是0,碰到这种情况该怎么办?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

不要用QQ问我东西

★ ★
zeoliters(金币+1,VIP+0):多谢回帖交流! 10-13 11:50
sxllgaoli898(金币+1,VIP+0):谢谢你的建议,good luck 10-14 17:05
引用回帖:
Originally posted by sxllgaoli898 at 2009-10-13 16:28:
可是考虑对称性的结果显示有四个虚频,那不是说明此分子不稳定吗?

说明你得到的构型 不是你要的 应该继续寻找稳定的构型
Chemistry[]==[]Chem[]is[]try!!!
4楼2009-10-13 17:31:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

半支烟5591

银虫 (正式写手)

专家


★ ★ ★ ★
mingdong(金币+1,VIP+0):感谢交流! 10-13 17:29
sxllgaoli898(金币+3,VIP+0):谢谢你的建议,good luck 10-14 17:05
当然考虑对称性时那说明这个分子是对称的,反之,不对称但是优化稳定点,偶极不为0,那理论上应该选用加对称性的计算!
2楼2009-10-13 16:10:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxllgaoli898

金虫 (小有名气)

可是考虑对称性的结果显示有四个虚频,那不是说明此分子不稳定吗?
3楼2009-10-13 16:28:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)

★ ★
zeoliters(金币+1,VIP+0):多谢回帖交流! 10-13 11:50
sxllgaoli898(金币+1,VIP+0):谢谢你的建议,good luck 10-14 17:05
方法不合适和限制的对称性不合适都可能导致这样的错误
一个人做点好事并不难难的是一辈子做好事!!!
5楼2009-10-13 17:50:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见