24小时热门版块排行榜    

查看: 385  |  回复: 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 的主题更新
信息提示
请填处理意见