[求助] 对Fe110面优化不收敛

» 猜你喜欢

大家能不能给点意见啊已经有30人回复
解某透露，生命学部的项目质量在评审中明显高于医学部已经有5人回复
AI派发是根据关键词、基金题目、研究者研究方向匹配？已经有11人回复
特聘教授，算高级职称吗？青年基金还没结题的情况下还能报面上吗？已经有18人回复
函评这一关首先干掉一大半，今年卷上天了已经有21人回复
领导好安排自己做杂事。在成果上给领导挂名，请其不要经常安排自己做杂事已经有8人回复
盲审已经有3人回复
校级优秀本科毕业论文，比例多少？看指导教师资历吗？已经有11人回复
请教下各位，南京航空航天大学天目湖校区如何？已经有5人回复
24-25申博已经有8人回复

1楼2018-08-31 21:07:35

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

in文件没上传成功？

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
附件 1 : Fe.relax.in

2018-08-31 21:12:15, 5.03 K

2楼2018-08-31 21:12:49

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

引用回帖:

2楼: Originally posted by guiwei2018 at 2018-08-31 21:12:49
in文件没上传成功？

这是俯视图

Fe.relax2.png

3楼2018-08-31 21:15:14

张玲海

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 3067.1
散金: 766
红花: 3
帖子: 398
在线: 1574.5小时
虫号: 1419918
注册: 2011-09-27
性别: GG
专业: 理论和计算化学

增大nbnd 数目，提高degauss，减少conv_thr

简单的生活

4楼2018-09-01 11:32:21

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

引用回帖:

4楼: Originally posted by 张玲海 at 2018-09-01 11:32:21
增大nbnd 数目，提高degauss，减少conv_thr

大神，我看到手册里说的是对于金属体系，for a metal, 20% more ，意思是超过电子数目一半的20%吗？比如480个电子，nbnd>240*1.2=288,我取了300不知道够不够。另外我把degauss从0.2增加到0.6，conv_thr减小到1.0*10-4

5楼2018-09-01 21:04:49

张玲海

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 3067.1
散金: 766
红花: 3
帖子: 398
在线: 1574.5小时
虫号: 1419918
注册: 2011-09-27
性别: GG
专业: 理论和计算化学

引用回帖:

5楼: Originally posted by guiwei2018 at 2018-09-01 21:04:49
大神，我看到手册里说的是对于金属体系，for a metal, 20% more ，意思是超过电子数目一半的20%吗？比如480个电子，nbnd>240*1.2=288,我取了300不知道够不够。另外我把degauss从0.2增加到0.6，conv_thr减小 ...

大概nbnd 为 380 ，conv_thr减小到1.0*10-5

简单的生活

6楼2018-09-01 23:54:08

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

引用回帖:

6楼: Originally posted by 张玲海 at 2018-09-01 23:54:08
大概nbnd 为 380 ，conv_thr减小到1.0*10-5...

好的，我再试试！谢谢

7楼2018-09-02 18:31:52

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

引用回帖:

6楼: Originally posted by 张玲海 at 2018-09-01 23:54:08
大概nbnd 为 380 ，conv_thr减小到1.0*10-5...

确实可以收敛。但是现在有两个问题，一是out文件里出现了SCF correction compared to forces is large: reduce conv_thr to get better values，其实我把conv_thr继续减小到1.0*10-6也可以收敛，但也会出现这样的错误。二是我的课题内容是做吸附，计算结合能，这就要求在优化“基底+吸附物”和‘基底’时，用同一套参数，现在是把优化“基底”的参数改了（增加nbnd=380,更换degauss=0.02为0.06），这样得到的结合能为-0.18与文献里-0.14不符合，请问大神，又该怎么办呢？

8楼2018-09-04 21:46:28

张玲海

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 3067.1
散金: 766
红花: 3
帖子: 398
在线: 1574.5小时
虫号: 1419918
注册: 2011-09-27
性别: GG
专业: 理论和计算化学

引用回帖:

8楼: Originally posted by guiwei2018 at 2018-09-04 21:46:28
确实可以收敛。但是现在有两个问题，一是out文件里出现了SCF correction compared to forces is large: reduce conv_thr to get better values，其实我把conv_thr继续减小到1.0*10-6也可以收敛，但也会出现这样的错 ...

degauss=0.06 优化完以后，提取优化以后的结构作为输入文件，然后把改成degauss=0.02 再优化一次。

简单的生活

9楼2018-09-04 21:52:06

guiwei2018

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.2
帖子: 36
在线: 30.2小时
虫号: 8080309
注册: 2018-02-24
性别: GG
专业: 机械摩擦学与表面技术

引用回帖:

9楼: Originally posted by 张玲海 at 2018-09-04 21:52:06
degauss=0.06 优化完以后，提取优化以后的结构作为输入文件，然后把改成degauss=0.02 再优化一次。...

实际情况是，即使一开始degauss=0.02,只需要增加nbnd=380就已经可以收敛了。是不是etot_conv_thr=1.0D-3和forc_conv_thr=1.0D-2设置得太大了？
Total force =    0.001442    Total SCF correction =    0.000399
   SCF correction compared to forces is large: reduce conv_thr to get better values

   BFGS Geometry Optimization

   bfgs converged in 1 scf cycles and 0 bfgs steps
   (criteria: energy <  1.0E-03 Ry, force <  1.0E-02 Ry/Bohr)

   End of BFGS Geometry Optimization

   Final energy = -8502.2661470662 Ry

10楼2018-09-05 10:11:46