24小时热门版块排行榜    

CyRhmU.jpeg
查看: 429  |  回复: 5
当前主题已经存档。

fu8777

木虫 (正式写手)

[交流] 【求助】g03优化亚油酸甲酯,出现问题

高手看看。用hyper chem画了一亚油酸甲酯的分子式,当时进行了半经验优化结构,但是1-4个原子的震动频率是负数。当然也就得到内坐标,稍微写了写,想用高斯算算,可是不能运算。

输入文件如下:

%chk=yayousuanjiazhi
%mem=1Gb
%nproc=2
%rwf=yayousuanjiazhi
# opt am1

opt

0 1
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
C                  6              B6    5              A5    4              D4
C                  7              B7    6              A6    5              D5
C                  8              B8    7              A7    6              D6
C                  9              B9    8              A8    7              D7
C                 10             B10    9              A9    8              D8
C                 11             B11   10             A10    9              D9
C                 12             B12   11             A11   10             D10
C                 13             B13   12             A12   11             D11
C                 14             B14   13             A13   12             D12
C                 15             B15   14             A14   13             D13
C                 16             B16   15             A15   14             D14
C                 17             B17   16             A16   15             D15
O                 18             B18   17             A17   16             D16
O                 18             B19   17             A18   16             D17
C                 20             B20   18             A19   17             D18
H                  1             B21    2             A20    3             D19
H                  1             B22    2             A21    3             D20
H                  1             B23    2             A22    3             D21
H                  2             B24    1             A23    3             D22
H                  2             B25    1             A24    3             D23
H                  3             B26    2             A25    1             D24
H                  3             B27    2             A26    1             D25
H                  4             B28    3             A27    2             D26
H                  4             B29    3             A28    2             D27
H                  5             B30    4             A29    3             D28
H                  5             B31    4             A30    3             D29
H                  6             B32    5             A31    4             D30
H                  7             B33    6             A32    5             D31
H                  8             B34    7             A33    6             D32
H                  8             B35    7             A34    6             D33
H                  9             B36    8             A35    7             D34
H                 10             B37    9             A36    8             D35
H                 11             B38   10             A37    9             D36
H                 11             B39   10             A38    9             D37
H                 12             B40   11             A39   10             D38
H                 12             B41   11             A40   10             D39
H                 13             B42   12             A41   11             D40
H                 13             B43   12             A42   11             D41
H                 14             B44   13             A43   12             D42
H                 14             B45   13             A44   12             D43
H                 15             B46   14             A45   13             D44
H                 15             B47   14             A46   13             D45
H                 16             B48   15             A47   14             D46
H                 16             B49   15             A48   14             D47
H                 17             B50   16             A49   15             D48
H                 17             B51   16             A50   15             D49
H                 21             B52   20             A51   18             D50
H                 21             B53   20             A52   18             D51
H                 21             B54   20             A53   18             D52

   B1             1.51211937
   B2             1.52001743
   B3             1.52017400
   B4             1.52154954
   B5             1.48765083
   B6             1.33483220
   B7             1.48959223
   B8             1.48837159
   B9             1.33443696
   B10            1.48861949
   B11            1.52122221
   B12            1.52051110
   B13            1.52027662
   B14            1.51957264
   B15            1.52223947
   B16            1.51910270
   B17            1.51138678
   B18            1.21702958
   B19            1.36666089
   B20            1.41221139
   B21            1.07002379
   B22            1.07000370
   B23            0.00000066
   B24            1.07000000
   B25            1.07000000
   B26            1.07000000
   B27            1.07000000
   B28            1.07000000
   B29            1.07000000
   B30            1.07000000
   B31            1.07000000
   B32            1.07000000
   B33            1.07000000
   B34            1.07000000
   B35            1.07000000
   B36            1.07000000
   B37            1.07000000
   B38            1.07000000
   B39            1.07000000
   B40            1.07000000
   B41            1.07000000
   B42            1.07000000
   B43            1.07000000
   B44            1.07000000
   B45            1.07000000
   B46            1.07000000
   B47            1.07000000
   B48            1.07000000
   B49            1.07000000
   B50            1.07000000
   B51            1.07000000
   B52            1.06876843
   B53            1.06917234
   B54            0.00009065
   A1           111.40834712
   A2           111.55530782
   A3           111.17173632
   A4           111.99970174
   A5           122.12124616
   A6           122.04207636
   A7           113.09543233
   A8           122.23389801
   A9           121.94090358
   A10          112.16605249
   A11          111.17812540
   A12          111.29181630
   A13          111.30483510
   A14          111.22813453
   A15          110.51283300
   A16          115.53529301
   A17          125.84788420
   A18          115.80218771
   A19          118.63412144
   A20           89.99992041
   A21          144.73560059
   A22           54.73737557
   A23          106.36365902
   A24          106.36365902
   A25          106.33201385
   A26          106.33201385
   A27          106.41456123
   A28          106.41456123
   A29          106.23618556
   A30          106.23618556
   A31          118.93937692
   A32          118.97896182
   A33          105.99903496
   A34          105.99903496
   A35          118.88305100
   A36          119.02954821
   A37          106.20026139
   A38          106.20026139
   A39          106.41318751
   A40          106.41318751
   A41          106.38873567
   A42          106.38873567
   A43          106.38593482
   A44          106.38593482
   A45          106.40243357
   A46          106.40243357
   A47          106.55600164
   A48          106.55600164
   A49          105.46662471
   A50          105.46662471
   A51           90.13351381
   A52          144.62898046
   A53           54.54212529
   D1           180.00000000
   D2           180.00000000
   D3           180.00000000
   D4           180.00000000
   D5           180.00000000
   D6           180.00000000
   D7           180.00000000
   D8           180.00000000
   D9           180.00000000
   D10          180.00000000
   D11          180.00000000
   D12          180.00000000
   D13          180.00000000
   D14          180.00000000
   D15          180.00000000
   D16          180.00000000
   D17            0.00000000
   D18          180.00000000
   D19            0.00000000
   D20          180.00000000
   D21          180.00000000
   D22          115.49982441
   D23         -115.49982441
   D24           64.48072302
   D25          -64.48072302
   D26          -64.53158497
   D27           64.53158497
   D28           64.42216163
   D29          -64.42216163
   D30            0.00000000
   D31            0.00000000
   D32          -64.27944512
   D33           64.27944512
   D34            0.00000000
   D35            0.00000000
   D36           64.40034121
   D37          -64.40034121
   D38          -64.53073538
   D39           64.53073538
   D40           64.51563080
   D41          -64.51563080
   D42          -64.51390281
   D43           64.51390281
   D44           64.52408818
   D45          -64.52408818
   D46          -64.61963897
   D47           64.61963897
   D48           63.97022482
   D49          -63.97022482
   D50            0.00000000
   D51          180.00000000
   D52          180.00000000
回复此楼
芳华尽放,韶华难求
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 11-16 13:11
怎么叫做不能算啊?总得有个说法啊。什么就叫做不能算呢?
抽空给你跑了一下,赫然发现,202就出错了。奇怪了。
再看log文件才知道。。。
你的1号C原子在相同的位置上隐藏了一个H原子(24号)
你的21号C原子在相同的位置上同样隐藏了一个H原子(55号)
原子重叠在一起,H被C覆盖了,没有显示出来。。
计算不出错才怪。
删掉24号和55号H,保证没问题了。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
2楼2009-11-13 22:50:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

yjcmwgk(金币+0,VIP+0):啥意思?扣你的金币,给zhangmt?哈哈哈哈 11-16 13:11
本帖仅楼主可见
3楼2009-11-15 19:44:17
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

fu8777

木虫 (正式写手)

非常谢谢2楼的帮助,可惜我不知道怎么给金币表示感谢。
芳华尽放,韶华难求
4楼2009-11-16 20:57:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fu8777

木虫 (正式写手)

非常谢谢1楼的帮助,可惜我不知道怎么给金币表示感谢
芳华尽放,韶华难求
5楼2009-11-16 20:58:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fu8777

木虫 (正式写手)

引用回帖:
Originally posted by fu8777 at 2009-11-3 20:59:
高手看看。用hyper chem画了一亚油酸甲酯的分子式,当时进行了半经验优化结构,但是1-4个原子的震动频率是负数。当然也就得到内坐标,稍微写了写,想用高斯算算,可是不能运算。

输入文件如下:

%chk=yayou ...

非常谢谢你
芳华尽放,韶华难求
6楼2009-11-16 21:01:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fu8777 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见