24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1550  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 太妃糖么么 的主题更新
信息提示
请填处理意见