24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1719  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 太妃糖么么 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211电子信息347求调剂 +3 554916 2026-04-03 4/200 2026-04-09 22:24 by 554916
[考研] 一志愿中南大学物理学,英一66,求调剂 +3 长烟旖旎 2026-04-08 4/200 2026-04-09 21:27 by wutongshun
[考研] 0703调剂 +18 拾玖壹 2026-04-04 20/1000 2026-04-09 20:53 by zhouxiaoyu
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 285求调剂 +12 AZMK 2026-04-05 18/900 2026-04-08 20:43 by 逆水乘风
[考研] 0703调剂,一志愿天津大学319分 +23 haaaabcd 2026-04-05 26/1300 2026-04-08 16:19 by luoyongfeng
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +20 努力奋斗112 2026-04-07 21/1050 2026-04-08 14:54 by screening
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 326分,一志愿沪9,求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 调剂求助(生物与医药) +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 13/650 2026-04-06 07:24 by houyaoxu
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
信息提示
请填处理意见