24小时热门版块排行榜    

查看: 381  |  回复: 1

fangyongxinxi

新虫 (正式写手)

[求助] 关于优化过程的进一步询问。

没有测试就没有发言权。关于优化过程的进一步询问.自己的INCAR 文件如下,是个二元体系,在优化结构的时候做了一些测试。
测试目的是如何选取ISIF,及如何实施优化过程。(看以前的帖子说要在优化过程中设置循环,找到能量最低点,然后再进行进一步优化,这种情况自己也做了测试,同样存在问题,以后再讨论。而且这种方法只适用于,体系只有一个内参数的情形。)。
思路:
1.直接用ISIF=3对INCAR文件进行数次优化,用优化后的CONCAR,当成下次优化的POSCAR,直到晶格参数变化在很小范围内,则算优化成功。
2.先用ISIF=2优化,在用isif = 3 优化;
3.先用ISIF=4优化,再用ISIF=3优化;
结果:
1. 每次优化,晶格参数都变化比较大。
第一次优化后, 以第一个参数为例(5.7516),变成5.6274,再优化一次,变成5.5636,再优化一次,变成 5.5283,再优化一次,变成 5.5115. 悲剧,晶格参数越来越小,体系能量随优化的进行,逐渐升高。体系稳定的有两个标准,受力和体系能量。但这个结果如何解释 ?
2 和3 优化出来的晶体参数一致,是否以后用着种方法就好了 ?


NiTi B2
   1.00000000000000
     5.7516         0.0000000000000000    0.0000000000000000
     0.0000000000000000     4.06762   0.0000000000000000
     0.0000000000000000    0.0000000000000000    4.06762
   Ni   Ti
   4   4
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5                   0.75                0.25
  0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5                   0.25                0.75
  0.25000000000000000  0.2500000000000000  0.2500000000000000
  0.75                0.25                0.25
  0.25000000000000000  0.7500000000000000  0.7500000000000000
  0.75                 0.75                0.75
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saikun

木虫 (正式写手)


gzqdyouxia(金币+1): 谢谢提示 2011-06-11 22:35:03
你看看pressure, 等于 +-0.0x, max force < 0.0x eV/A 就对了。
2楼2011-06-11 14:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fangyongxinxi 的主题更新
信息提示
请填处理意见