24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2139  |  回复: 8

qingjing99

木虫 (小有名气)

[求助] 过渡态优化不收敛,急!!!

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

跟体系大有关系(129个原子)吗?scf=qc可行吗

[ Last edited by qingjing99 on 2011-10-11 at 16:01 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

【答案】应助回帖

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

然后调大优化步长, 比如 iop 调到4
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
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找到,干脆就用他们不行吗?
5楼2011-10-11 22:41:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vallen

银虫 (正式写手)

【答案】应助回帖

qingjing99(金币+4): 谢谢 2011-10-12 10:14:35
scf=qc是用于SCF不收敛的,你这个不用
还是换方法吧
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 288求调剂 +13 代fish 2026-04-09 14/700 2026-04-10 23:23 by babysonlkd
[考研] 297求调剂 +11 ORCHID1 2026-04-10 12/600 2026-04-10 21:07 by cyh—315
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 312求调剂 +4 Say Never 2026-04-04 4/200 2026-04-08 08:41 by barlinike
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 第一志愿东南大学物理313,有科研竞赛获奖经历,希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 083200 333求调剂 +3 十二!! 2026-04-04 3/150 2026-04-05 08:28 by barlinike
信息提示
请填处理意见