24小时热门版块排行榜    

查看: 1879  |  回复: 20

chemistryxrw

金虫 (小有名气)

引用回帖:
Originally posted by fatpig8832 at 2011-05-08 14:11:21:
算td的时候是否用了 guess=read 来读取 opt 的 chk 文件中的初猜?

如果直接把构型放进去用td算,有可能会不一样的...

(⊙o⊙)…。。。这个问题,我倒只是用out文件重新生成个gjf文件然后改命令算td的。。。
动力学
11楼2011-05-08 14:21:56
已阅   回复此楼   关注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的回帖

鱼妃

金虫 (著名写手)

挖坑不填

引用回帖:
Originally posted by fatpig8832 at 2011-05-08 14:11:21:
算td的时候是否用了 guess=read 来读取 opt 的 chk 文件中的初猜?

如果直接把构型放进去用td算,有可能会不一样的...

是否使用了guess(read)有区别吗?不管是否使用了guess(read),楼主的初猜都是从401到502的初猜,楼主的fock矩阵都是通过502迭代得到的。怎么就不一样呢?除非波函数极端不稳定呢,否则怎么会不一样呢?如果楼主的操作无误,偶继续表示亚历山大
己所不欲,勿施于人;己所欲者,亦勿施人。无论穷达,但守本分;八风不动,独善自身。
13楼2011-05-08 17:00:33
已阅   回复此楼   关注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的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

引用回帖:
Originally posted by fatpig8832 at 2011-05-08 17:28:57:
据我估计,楼主的opt的初始构型和结束时的构型,这两个的构型可能相差较远,这样在opt时,第一步的初猜跟最后一步可能是很不相同的...就是说,opt的结果的轨道,是要经过一番乱搞狂算才能弄出来的...

而做t ...

同一个几何,同一理论水平,就因为一个是TD,一个是OPT就不一样了?
前途光明,出路难觅!
15楼2011-05-09 08:32:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pw91

铜虫 (小有名气)

引用回帖:
Originally posted by lihb734 at 2011-05-09 08:32:51:
同一个几何,同一理论水平,就因为一个是TD,一个是OPT就不一样了?

不是没有可能啊
如果初猜的波函数太差 反正波函数也是SCF得到的 说不定就猜到别的上面去了呢 不过我总觉得概率较小
其实我还是不太相信楼主自称的“操作无误”
16楼2011-05-09 08:38:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

九仙居士

捐助贵宾 (正式写手)

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

建议试试把优化后的坐标取出来进行td计算
17楼2011-05-09 10:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nybz1986lizi

铁虫 (小有名气)

【答案】应助回帖

我是用opt的chk文件生成td的gjf文件的!然后直接改关键词得到的能量并没有差别,附上图形!希望对楼住有所帮助!!!

18楼2011-05-09 14:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemistryxrw

金虫 (小有名气)

引用回帖:
Originally posted by nybz1986lizi at 2011-05-09 14:53:57:
我是用opt的chk文件生成td的gjf文件的!然后直接改关键词得到的能量并没有差别,附上图形!希望对楼住有所帮助!!!

一个是opt得到的chk文件生成的MO,另一个是td的chk文件生成的MO吗?
动力学
19楼2011-05-09 16:43:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

但个

铜虫 (初入文坛)

引用回帖:
711630楼: Originally posted by yjcmwgk at 2011-05-07 22:05:28:
楼主斩钉截铁的回答 让我对自己大失信心 难道高斯的td有什么特殊之处吗?难道高斯公司那帮神人修改里td理论吗? 偶表示压力很大
赶紧计算了一下 发现楼主肯定是自己搞错了 斌兄 咱可以撤了 哈哈
e6/a3/448368_ ...

斑竹,直接打开log或者out文件另存为gif和直接用查找log或者out文件里面优化后的文件的坐标(查找checkpoint,会快一些的,)有一点点不一样的,结算后会有差别的
20楼2012-04-26 10:51:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chemistryxrw 的主题更新
信息提示
请填处理意见