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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 量化新手直通车 » 求大神帮我看看，我的计算思路有没有问题？
61/1返回列表
查看: 1407  |  回复: 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

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +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的回帖

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的回帖

pengyj

至尊木虫 (著名写手)
对，当时考虑错了。

那楼主的应该是局部极小的问题了。
赞一下 回复此楼
6楼2013-01-04 22:09:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dokinght 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 7/350 2026-09-03 18:50 by jyNB2vrbdUpR
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 18:50 by jyNB2vrbdUpR
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 YLHlRHNwYkce 2026-09-02 11/550 2026-09-03 18:42 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 rFEsUBKXRll0 2026-09-02 6/300 2026-09-03 16:49 by T0rGB46095mJ
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 16:38 by T0rGB46095mJ
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 14:48 by T0rGB46095mJ
[基金申请] 国社科又开始会评了，不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗？ +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
关闭
请填处理意见
关闭