24小时热门版块排行榜    

查看: 1390  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dokinght

新虫 (初入文坛)

[求助] 求大神帮我看看,我的计算思路有没有问题?

我根据文献,模拟图中所示的化合物,采用和文献中的同样的计算方法(B3LYP/6-31G*),我的思路是这样的。
1.用chem3d画出此化合物,然后MM2优化一下。
2.输出成高斯03文件,然后用高斯在B3LYP/6-31G*下结构优化。
3.用FREQ查看无虚频。
4.用pop=full 查看LUMO和HOMO。

  问题是,不管是最终算出的能量 还是LUMO和HOMO的值与文献的值都有差距。所以想问大神们,我的思路有没有问题?
如果有问题,那么最佳的优化步骤是怎么样的?
如果没有问题,那么这种差距是计算机本身的原因,还是什么?

未命名.jpg
回复此楼

» 猜你喜欢

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

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

pengyj

至尊木虫 (著名写手)

对,当时考虑错了。

那楼主的应该是局部极小的问题了。
6楼2013-01-04 22:09:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

virtualzx

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-12-31 16:19:01
检查你的收敛和积分精度,如果都没问题,而你得到的能量比文献能量高,可能是优化到了局部极小,尝试直接从文献结构开始优化。
2楼2012-12-31 12:36:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyz7819227

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-12-31 16:19:05
你可以先用CCDC来查一下这个化合物的晶体结构,如果没有的话,那就直接按照文献上的最终结构来画。  直接用GV画,不用chem3D,画完就直接opt+freq优化,HOMO和LUMO在输出文件里面自然能看到。
海纳百川,有容乃大
3楼2012-12-31 13:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengyj

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
无虚频的话就可以排除局部极小了吧?

6-31G*不算太大的基组,从文献的最终结构来优化再来对比下看看。
4楼2013-01-01 23:09:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 10:33 by winsaint
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:13 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
信息提示
请填处理意见