24小时热门版块排行榜    

查看: 1994  |  回复: 12

太妃糖么么

金虫 (正式写手)

[求助] tddft做优化,convergence failed

tddft做优化,采用优化基态构型,但是出来说的是convergence failed,看了一下是能量再循环过程中震荡了,想问一下有什么好的解决方法吗?
%NprocShared=2
%mem=2000MW
%chk=HP1-0-OPP-OPP-RH-td-opt.chk
#p ub3lyp/6-311++g(d,p) td(nstate=6,singlets,root=2) pop=full gfinput
opt

HP1-0-OPP-OPP-RH-td-opt

C                    -1.18613   0.43894   0.33366
H                    -1.05065   0.38844   1.42555
H                    -2.04181   1.09039   0.14402
C                     0.06722   1.02658  -0.31817
H                    -0.03952   0.97429  -1.40663
H                     0.15947   2.08658  -0.05252
C                     1.36256   0.31799   0.09824
H                     2.25631   0.733    -0.39835
O                    -1.5415   -0.83438  -0.18961
H                    -0.76285  -1.40157  -0.11998
H                     1.54967   0.49319   1.18294
O                     1.34995  -1.0488    0.00743
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

枪下游魂

木虫 (著名写手)

★ ★
gmy1990: 金币+2 2013-10-25 20:38:24
引用回帖:
3楼: Originally posted by 太妃糖么么 at 2013-10-22 11:50:35
还是不行啊。我已经设过了opt(maxstep=10),你说的的是基态的用小精度算一下吗?用tddft优化是不是一定可以得到构型啊?...

激发态的构型用小方法优化后,读入优化好的chk用你的方法继续优化。
你做的是opt,优化后肯定可以得到构型。
可设maxstep=1试试。
或者不在基态稳定结构的基础上算激发态,有时是会收敛。
你既然说了能量震荡,那应该是l9999错误?
另外你确定你要研究的是第二激发态吧?
4楼2013-10-22 12:03:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

★ ★
gmy1990: 金币+2 2013-10-25 20:39:14
很多人认为优化一直震荡应该用小步长来优化,其实再下不这么认为,一直在震荡表面你这个结构处于一个具有较小能垒的“谷底”中,而这个谷底并不是全局最优点。你需要做的是加很大的步长让结构跳出这个能垒进行优化,或者微调分子结构进行优化。 激发态优化中很多时候是H原子变化,或者柔性角的转到,你调结构的时候微调这两个地方就行
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
10楼2013-10-25 10:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-10-22 09:45:47
太妃糖么么: 金币+10, 有帮助 2013-10-22 16:54:08
用opt(maxstep=x)减小步长试试看?
或先用小基组低精度方法算一下,再读入构型用你本来的方式算算看?
2楼2013-10-22 09:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

太妃糖么么

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2013-10-22 09:30:16
用opt(maxstep=x)减小步长试试看?
或先用小基组低精度方法算一下,再读入构型用你本来的方式算算看?

还是不行啊。我已经设过了opt(maxstep=10),你说的的是基态的用小精度算一下吗?用tddft优化是不是一定可以得到构型啊?
3楼2013-10-22 11:50:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu0413

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
看看算到最后一步的时候 第二激发态跟其他激发态之间有没有交叉
如果有交叉的话 那个位置的能量不可导 也会报错的
5楼2013-10-22 12:23:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

太妃糖么么

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 枪下游魂 at 2013-10-22 12:03:21
激发态的构型用小方法优化后,读入优化好的chk用你的方法继续优化。
你做的是opt,优化后肯定可以得到构型。
可设maxstep=1试试。
或者不在基态稳定结构的基础上算激发态,有时是会收敛。
你既然说了能量震荡, ...

是的,我想问问,maxstep是不是等于iop(1/8=N)
6楼2013-10-22 15:55:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

太妃糖么么

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 枪下游魂 at 2013-10-22 12:03:21
激发态的构型用小方法优化后,读入优化好的chk用你的方法继续优化。
你做的是opt,优化后肯定可以得到构型。
可设maxstep=1试试。
或者不在基态稳定结构的基础上算激发态,有时是会收敛。
你既然说了能量震荡, ...

不是l9999错误,只是写的not converged
7楼2013-10-22 15:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

引用回帖:
6楼: Originally posted by 太妃糖么么 at 2013-10-22 15:55:05
是的,我想问问,maxstep是不是等于iop(1/8=N)...

从说明上来看是一样的
8楼2013-10-22 16:25:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

引用回帖:
7楼: Originally posted by 太妃糖么么 at 2013-10-22 15:56:16
不是l9999错误,只是写的not converged...

可能SCF not converged
甲烷无敌
9楼2013-10-25 08:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 太妃糖么么 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见