24小时热门版块排行榜    

查看: 2341  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见