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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 过渡态优化不收敛，急！！！
91/1返回列表
查看: 2389  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

qingjing99

木虫 (小有名气)

[求助] 过渡态优化不收敛，急！！！

我用g09优化一个过渡态，总是不收敛。优化初始阶段结构挺好，力也逐渐减小（甚至有两个yes出现），但继续优化就会偏离原来的构型，力也突然增大。我取出较好的构型计算力常熟优化或者减小步长过渡态优化都不行。用m05，m06方法都可以得到满意的结构，但用回原来的dft方法还是会有同样的情况。急盼高手出手相救啊，课题就差这一步了。。。

跟体系大有关系（129个原子）吗？scf=qc可行吗

[ Last edited by qingjing99 on 2011-10-11 at 16:01 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-11 15:59:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

【答案】应助回帖

首先看一下你的输入文件吧， 我估计你的输入文件就不是很好的构型。

然后调大优化步长， 比如 iop 调到4
赞一下(1人) 回复此楼
我想做个科学家，可我又喜欢写诗；我想做一个诗人，可我却读的理工科..
2楼2011-10-11 16:47:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingjing99

木虫 (小有名气)
引用回帖:
2楼: Originally posted by jiewei at 2011-10-11 16:47:17:
首先看一下你的输入文件吧， 我估计你的输入文件就不是很好的构型。

然后调大优化步长， 比如 iop 调到4

我用的是QST2方法优化计算，然后取较好的点出来再用ts方法。还有，我用m05，m06方法算好的结构作为初始构型也不行啊。。。
赞一下 回复此楼
3楼2011-10-11 17:51:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

【答案】应助回帖

qingjing99(金币+8): 受益 2011-10-12 10:12:22
引用回帖:
3楼: Originally posted by qingjing99 at 2011-10-11 17:51:55:
我用的是QST2方法优化计算，然后取较好的点出来再用ts方法。还有，我用m05，m06方法算好的结构作为初始构型也不行啊。。。

初始构型的好坏和选用什么方法优化的构型有什么必然联系么？

你方法完全合理就能优化出合理的构型么？


一个有机分子常常有很多种亚稳态，而各个亚稳态之间在优化时不一定跳得过去，
很多时候优化结果未必就是能量最低的点，所以你的输入构型才显得重要。
这和方法无关，仅仅和你输入构型有关。。
赞一下 回复此楼
我想做个科学家，可我又喜欢写诗；我想做一个诗人，可我却读的理工科..
4楼2011-10-11 18:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangth0921

铁虫 (职业作家)

【答案】应助回帖

qingjing99(金币+8): 受益 2011-10-12 10:13:00
可能是DFT方法的问题。例如，很多情况下，最常用的B3LYP就找不到H-原子转移的过渡态，而BH&HLYP或BB1K等就可以。既然你能在m05和m06找到，干脆就用他们不行吗？
赞一下(2人) 回复此楼
5楼2011-10-11 22:41:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vallen

银虫 (正式写手)

【答案】应助回帖

qingjing99(金币+4): 谢谢 2011-10-12 10:14:35
scf=qc是用于SCF不收敛的，你这个不用
还是换方法吧
赞一下(1人) 回复此楼
6楼2011-10-11 23:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingjing99

木虫 (小有名气)
引用回帖:
4楼: Originally posted by jiewei at 2011-10-11 18:03:15:
初始构型的好坏和选用什么方法优化的构型有什么必然联系么？

你方法完全合理就能优化出合理的构型么？


一个有机分子常常有很多种亚稳态，而各个亚稳态之间在优化时不一定跳得过去，
很多时候优化结果未 ...

虽然不太同意您的观点，不过还是谢谢啦。
赞一下 回复此楼
7楼2011-10-12 10:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingjing99

木虫 (小有名气)
引用回帖:
5楼: Originally posted by wangth0921 at 2011-10-11 22:41:06:
可能是DFT方法的问题。例如，很多情况下，最常用的B3LYP就找不到H-原子转移的过渡态，而BH&HLYP或BB1K等就可以。既然你能在m05和m06找到，干脆就用他们不行吗？

其他结构都算好了，就差这一个了，如果换方法，要全部重算啊。。。
赞一下 回复此楼
8楼2011-10-12 10:14:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦话

新虫 (正式写手)
引用回帖:
2402672楼: Originally posted by jiewei at 2011-10-11 16:47:17
首先看一下你的输入文件吧， 我估计你的输入文件就不是很好的构型。

然后调大优化步长， 比如 iop 调到4

我想问一个问题，我算激发态优化时用的如下方法#p hf/6-31+g(d) cis(singlets,nstates=5) scrf(pcm,solvent=CH3CH2OCH2CH3
) opt optcyc=300 maxdisk=20GB
出现了下面的结果
Root      1 :  ********************   Change is ********************
Root      2 :  ********************   Change is -177.750956439561200
Root      3 :  ********************   Change is -109.476836921007000
Root      4 :  ********************   Change is   -6.854283804947429
Root      5 :  ********************   Change is  -58.866471450604870
这是怎么回事呢？请求你的帮助，我都算了一个多月了，还没出来，还不报错。
赞一下 回复此楼
用微笑收集快乐的金字塔
9楼2012-08-16 08:48:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qingjing99 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭