24小时热门版块排行榜    

查看: 1663  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 太妃糖么么 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +9 大米饭! 2026-03-15 12/600 2026-03-16 21:59 by 大米饭!
[文学芳草园] 伙伴们,祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 本人考085602 化学工程 专硕 +12 不知道叫什么! 2026-03-15 14/700 2026-03-16 16:45 by 我的船我的海
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见