24小时热门版块排行榜    

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

chemistryxrw

金虫 (小有名气)

[求助] opt分子轨道图与td生成的分子轨道图不一样????

1.我用优化结构后的out文件来生成一个gjf文件,用这个gjf文件来生成分子轨道;
2.使用算过td后的chk文件生成fchk文件后,算出的分子轨道;

为什么这两种算出来的分子轨道图不一样呢?
回复此楼
动力学
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

算td的时候是否用了 guess=read 来读取 opt 的 chk 文件中的初猜?

如果直接把构型放进去用td算,有可能会不一样的...
10楼2011-05-08 14:11:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
Originally posted by chemistryxrw at 2011-05-08 14:21:56:
(⊙o⊙)…。。。这个问题,我倒只是用out文件重新生成个gjf文件然后改命令算td的。。。

这样的话...有可能会不一样的...

一般算td都会用 guess=read 来读取初猜吧...

把你opt的chk文件copy成td的chk文件,然后在td的input file里加上guess=read来读取初猜,同时加上geom=allcheck来读取自旋多重度及构型,不知这样算是否能使两者的分子轨道一致?
12楼2011-05-08 14:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
Originally posted by 鱼妃 at 2011-05-08 17:00:33:
是否使用了guess(read)有区别吗?不管是否使用了guess(read),楼主的初猜都是从401到502的初猜,楼主的fock矩阵都是通过502迭代得到的。怎么就不一样呢?除非波函数极端不稳定呢,否则怎么会不一样呢?如果楼主 ...

据我估计,楼主的opt的初始构型和结束时的构型,这两个的构型可能相差较远,这样在opt时,第一步的初猜跟最后一步可能是很不相同的...就是说,opt的结果的轨道,是要经过一番乱搞狂算才能弄出来的...

而做td时,虽然方法和基组相同,但初猜的结构换了个新的,用b3直接猜很可能就猜到另外一个轨道上去了,而不会猜到那个在opt中经过一番计算才得到的轨道。据我的经验,b3在这种情况下猜错的几率还是挺高,用 hf 的话就好得多...

当然也可能是波函数不稳定的问题,因为楼主在opt中并没有加stable...
14楼2011-05-08 17:28:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chemistryxrw 的主题更新
信息提示
请填处理意见