24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 666  |  回复: 11
当前主题已经存档。

huangyc

木虫 (正式写手)

[交流] 【求助】TS优化问题

cNEB后,取出IBRION=1优化,POTIM大致依手册确定。发现结构变化不大,能量变化也非常小(时升时降,如下),但就是某些原子上的力达不到设定的收敛标准,不知如何处理?请大家支招
   1 F= -.29994723E+03 E0= -.29994977E+03  d E =-.299947E+03
   2 F= -.29994680E+03 E0= -.29994934E+03  d E =0.432291E-03
   3 F= -.29994746E+03 E0= -.29995000E+03  d E =-.660075E-03
   4 F= -.29994800E+03 E0= -.29995054E+03  d E =-.544400E-03
   5 F= -.29994890E+03 E0= -.29995144E+03  d E =-.893707E-03
   6 F= -.29994894E+03 E0= -.29995147E+03  d E =-.374641E-04
   7 F= -.29994896E+03 E0= -.29995150E+03  d E =-.254279E-04
   8 F= -.29995099E+03 E0= -.29995353E+03  d E =-.203132E-02
   9 F= -.29995099E+03 E0= -.29995353E+03  d E =-.905478E-06
  10 F= -.29995347E+03 E0= -.29995601E+03  d E =-.247985E-02
  11 F= -.29995425E+03 E0= -.29995679E+03  d E =-.779520E-03
  12 F= -.29995462E+03 E0= -.29995716E+03  d E =-.364252E-03
  13 F= -.29994968E+03 E0= -.29995222E+03  d E =0.493411E-02
  14 F= -.29994952E+03 E0= -.29995206E+03  d E =0.164753E-03
  15 F= -.29994882E+03 E0= -.29995136E+03  d E =0.694011E-03
  16 F= -.29994882E+03 E0= -.29995136E+03  d E =-.359578E-07
  17 F= -.29994854E+03 E0= -.29995107E+03  d E =0.280731E-03
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guitar2033

至尊木虫 (职业作家)

★ ★
huangyc(金币+1):谢谢参与
wuchenwf(金币+1,VIP+0):谢谢回帖交流 10-10 22:09
是不是你的初始结构离稳定结构较远
可以加大POTIM,或换一种IBRION?
2楼2009-09-08 09:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlshi

金虫 (正式写手)

★ ★
huangyc(金币+1):谢谢参与
wuchenwf(金币+1,VIP+0):谢谢回帖交流 10-10 22:09
用IBRION=2
引用回帖:
Originally posted by huangyc at 2009-9-7 19:38:
cNEB后,取出IBRION=1优化,POTIM大致依手册确定。发现结构变化不大,能量变化也非常小(时升时降,如下),但就是某些原子上的力达不到设定的收敛标准,不知如何处理?请大家支招
   1 F= -.29994723E+03 E0= - ...

3楼2009-09-08 09:44:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)

用2的话就会很容易跑回is或fs,同样POTIM设的过大也会出现此情况

不知道有没有其它控制参数
4楼2009-09-08 10:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


huangyc(金币+1):谢谢参与
TS为啥还要优化?
5楼2009-09-08 10:52:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenxuan5788

金虫 (正式写手)

★ ★
huangyc(金币+1):谢谢参与
wuchenwf(金币+1,VIP+0):谢谢回帖交流 10-10 22:09
建议你用这个网站提供的程序来计算,我一直在用,挺不好用的。具体设置网站上有,而且出现问题了你可以再改网站论坛上咨询
http://theory.cm.utexas.edu/henkelman/
6楼2009-09-08 11:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)

回LS:我是用cNEB在做。至于其它方法我试过dimer,lanczons等。现在是针对我这个问题应如何解决

回4楼:.......
7楼2009-09-08 11:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
huangyc(金币+1):谢谢参与
qasd(金币+2,VIP+0):xiexie~ 9-9 13:38
huangyc(金币+4,VIP+0):xie 9-11 14:26
试着减小PTOIM,将该参数设置为0.1
此外,如果cNEB如果收敛了的话,可以不对TS再做优化,直接计算TS态的振动频率可判断究竟是不是TS态。
8楼2009-09-08 12:29:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minglei288

金虫 (小有名气)

★ ★
huangyc(金币+1):谢谢参与
qasd(金币+1,VIP+0):xiexie~ 9-9 13:39
IBRION=3 试试吧,再调调初始构型
9楼2009-09-08 13:40:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)

看能量时高时低,可能是POTIM过大或过小。我的POTIM已设为0.1,现在我把它调成0.05了,结果也未出来。cNEB没有收敛,耗不起那么长时间

3我没试过,或许有可能?
10楼2009-09-08 14:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huangyc 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见