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

24小时热门版块排行榜    

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

pengyj

至尊木虫 (著名写手)

【答案】应助回帖

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

6-31G*不算太大的基组，从文献的最终结构来优化再来对比下看看。
赞一下 回复此楼
4楼2013-01-01 23:09:55
已阅   回复此楼   关注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的回帖

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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭