24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1938  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

bjwang

铁杆木虫 (正式写手)

[求助] siesta结构优化收敛标准

用sieata优化一个结构,优化部分设置如下:
MD.TypeOfRun          cg
MD.NumCGsteps         150
MD.MaxCGDispl         0.1  Ang
MD.MaxForceTol        0.02 eV/Ang
MD.VariableCell       F
scf部分:
MaxSCFIterations     300      # Maximum number of SCF iter
DM.MixingWeight      0.25     # New DM amount for next SCF cycle
DM.NumberPulay       10
DM.Tolerance         1.d-4    # Tolerance in maximum difference between input and output DM
DM.UseSaveDM         T        # to use continuation files
DM.MixSCF1           F
DM.PulayOnFile       F        # Store in memory ('F') or in files ('T')
SolutionMethod       diagon   # OrderN or Diagon
ElectronicTemperature  100.0 K  #Temp. for Fermi smearing

但58步正常结束后发现:
1.各个原子受力里面有0.03左右的:
siesta: Atomic forces (eV/Ang):
siesta:     27    0.022356    0.023712    0.000000
2.总力也不小于0.02呀
siesta:    Tot   -0.041567    0.000155    0.000000
问题:优化收敛标准不是是看输出文件中的siesta: Atomic forces (eV/Ang):项吗?
如果是,为什么会在最大步数达到前,还没达到收敛标准就正常结束呢?
谢谢!
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

It is a shame for the soul to be first to give way in this life,when thy body does not give way.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-12-08 18:08:28
bjwang: 金币+5, 有帮助, 谢谢交流 2013-01-03 11:09:45
可能你对某些维度进行了限制。这个时候,只看不受限制的自由度的受力。
弘德明志博学笃行
2楼2012-12-08 17:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 鼓励交流。。 2012-12-24 21:40:26
引用回帖:
4楼: Originally posted by yangxuezhang at 2012-12-24 15:14:22
麻烦问一下,关于收敛,有人告诉我dDmax这项如果一直比较大,就说明不收敛,请问这个说法对不对?...

当然对了。另外这个“一直比较大”要想对你目标精度而言的。
弘德明志博学笃行
5楼2012-12-24 15:21:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bjwang 的主题更新
信息提示
请填处理意见