24小时热门版块排行榜    

查看: 2471  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:31 by lXgU0RvNTDC9
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:22 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 18:55 by lXgU0RvNTDC9
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见