24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1712  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 太妃糖么么 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +9 chenxrlkx 2026-04-05 9/450 2026-04-06 09:29 by 学员tURuqU
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 086000生物与医药 初试274求调剂 +6 小叮当来了 2026-03-30 7/350 2026-04-05 20:30 by lys0704
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 290求调剂085701 +10 1314捧花 2026-04-02 10/500 2026-04-05 10:19 by Sealedwind
[考研] 材料科学与工程调剂 +19 深V宿舍吧 2026-03-30 20/1000 2026-04-04 22:13 by hemengdong
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 一志愿a区211,085601-307分求调剂 +13 党嘉豪 2026-03-31 26/1300 2026-04-03 08:33 by 495374996
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
信息提示
请填处理意见