24小时热门版块排行榜    

查看: 4972  |  回复: 11

dami251

银虫 (小有名气)

引用回帖:
9楼: Originally posted by gengcy at 2011-11-08 08:51:20:
用b3lyp做非键作用要慎重,如果非要用,一定要做BSSE,不然审稿时可能会怀疑你结果。报错的原因是输入结构的问题,建议做点技巧,比如先固定一部分原子优化,再全部放开优化整体或者,慢慢往上增加分子,具体只能 ...

非常抱歉,最近一直在做实验,所以没顾上回复。我现在尝试计算酒石酸酯和布洛芬的分子的相互作用,即结合能。具体分子形式见附件。感觉这两分子主要是通过氢键作用的。文献中有用B3LYP来计算氢键的,所以我就用B3LYP来算,但是计算顺利完成的很少,总是报错。。所以请教大侠,哪一些关键词适合计算我的体系啊,,苦思中。。两分子复合物的输入文件见下:
%Chk=L-tar-S-ibu.chk
%mem=6000MB
%nproc=8
#p b3lyp/6-31g opt

[No Title]

0   1
C  0          0.0000         0.0000         0.7584
C  0          0.0000         0.0000        -0.7584
H  0          0.1255         1.0199         1.1211
H  0         -0.1866        -1.0103        -1.1214
O  0         -1.2664        -0.5713         1.1940
H  0         -1.4355        -0.3165         2.1173
O  0         -1.0362         0.9257        -1.1936
H  0         -1.2825         0.7275        -2.1134
C  0          1.3198         0.5073        -1.3120
O  0          2.1888         1.1049        -0.6750
O  0          1.3855         0.2854        -2.6615
C  0          2.5802         0.7921        -3.3815
C  0          2.4344         0.3699        -4.8345
H  0          3.4640         0.3576        -2.9133
H  0          2.6154         1.8724        -3.2492
C  0          3.6627         0.7293        -5.6948
H  0          2.2742        -0.7105        -4.8617
H  0          1.5379         0.8334        -5.2560
H  0          3.5290         0.2791        -6.6835
H  0          4.5567         0.2608        -5.2671
C  0          3.9109         2.2390        -5.8675
C  0          5.0767         2.5544        -6.8207
H  0          4.1167         2.7011        -4.8958
H  0          2.9975         2.7134        -6.2460
C  0          5.3312         4.0595        -6.9880
H  0          5.9857         2.0691        -6.4489
H  0          4.8706         2.1076        -7.7996
H  0          6.1598         4.2487        -7.6734
H  0          4.4483         4.5664        -7.3857
H  0          5.5782         4.5263        -6.0311
C  0          1.1069        -0.8798         1.3127
O  0          1.7660        -1.7015         0.6736
O  0          1.2241        -0.6985         2.6648
C  0          2.2087        -1.5403         3.3889
C  0          2.2525        -1.0345         4.8216
H  0          3.1627        -1.4535         2.8724
H  0          1.8753        -2.5774         3.3285
C  0          3.1568        -1.8845         5.7373
H  0          2.5930         0.0044         4.8225
H  0          1.2339        -1.0334         5.2177
H  0          3.0345        -1.5253         6.7639
H  0          2.8050        -2.9229         5.7357
C  0          4.6535        -1.8533         5.3768
C  0          5.5318        -2.6126         6.3866
H  0          4.9885        -0.8109         5.3191
H  0          4.8113        -2.2820         4.3813
C  0          7.0247        -2.5916         6.0267
H  0          5.3896        -2.1789         7.3825
H  0          5.1867        -3.6502         6.4532
H  0          7.6186        -3.1364         6.7634
H  0          7.2017        -3.0514         5.0512
H  0          7.4050        -1.5679         5.9842
C  0          4.1258        -2.8706         0.0243
C  0          2.9189        -2.1587        -0.0387
C  0          2.8963        -0.7637        -0.0868
C  0          4.0915        -0.0313        -0.0728
C  0          5.3013        -0.7371        -0.0160
C  0          5.3186        -2.1315         0.0333
C  0          4.7095        -4.9158         1.4169
C  0          4.1309        -4.3950         0.0823
C  0          4.8691        -5.0185        -1.1225
C  0          4.1038         1.4998        -0.0918
C  0          3.3918         2.1158         1.1334
O  0          4.2120         2.4843        -2.3283
C  0          3.5488         2.0312        -1.3960
O  0          2.1677         1.9594        -1.4608
H  0          1.9832        -2.7049        -0.0522
H  0          1.9503        -0.2449        -0.1521
H  0          6.2368        -0.1914        -0.0117
H  0          6.2714        -2.6436         0.0769
H  0          4.6478        -6.0056         1.4629
H  0          5.7595        -4.6350         1.5270
H  0          4.1626        -4.5068         2.2681
H  0          3.0863        -4.7201         0.0295
H  0          4.7948        -6.1081        -1.0950
H  0          4.4444        -4.6711        -2.0656
H  0          5.9296        -4.7565        -1.1157
H  0          5.1458         1.8200        -0.0917
H  0          3.4652         3.2048         1.1206
H  0          3.8595         1.7499         2.0476
H  0          2.3376         1.8462         1.1499
H  0          1.8534         2.2847        -2.3269

R-布洛芬



L-酒石酸酯



二者复合物

11楼2011-11-17 17:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dami251

银虫 (小有名气)

引用回帖:
8楼: Originally posted by baoyuan6969 at 2011-11-08 00:32:51:
如果是vdW相互作用建议用加色散作用的泛函做,
报错可能和你的输入文件写的有关系

请问您用B3LYP计算包含氢键作用的分子间相互作用力如何呢?我现在就是想模拟手性识别作用。具体分子见我回复上面的内容。如果可以使用该方法的话,该再加些什么限定呢?用B3lyp/6-31G不太行,老报错。。我是学化工的,一下子还真找不到解决办法。急求啊。。。
12楼2011-11-17 19:00:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dami251 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 14/700 2026-11-09 04:20 by l48hFd8Nwae2
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 04:16 by l48hFd8Nwae2
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 04:09 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 03:58 by l48hFd8Nwae2
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见