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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 量化新手直通车 » 求大神帮我看看，我的计算思路有没有问题？
51/1返回列表
查看: 1385  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dokinght

新虫 (初入文坛)

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

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

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

未命名.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-12-31 10:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)
引用回帖:
4楼: Originally posted by pengyj at 2013-01-01 10:09:55
无虚频的话就可以排除局部极小了吧？

6-31G*不算太大的基组，从文献的最终结构来优化再来对比下看看。

局部极小是没有虚频的，鞍点才有。
赞一下 回复此楼
5楼2013-01-03 00:01: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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-29 23:14 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:35 by 4qO8zhfjr2AX
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:22 by 4qO8zhfjr2AX
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-29 21:22 by 4qO8zhfjr2AX
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:11 by 4qO8zhfjr2AX
[基金申请] 为什么到现在没收到通知？ +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
关闭
请填处理意见
关闭