应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » RB3LYP和B3LYP有何区别
122/2返回列表上一页12
查看: 4854  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 的主题更新
122/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯，如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭