计算化学公社

标题: VASP采用虚晶近似（VCA）计算gamma-Fe2O3能量为正值 [打印本页]

作者
Author: I10140317 时间: 2024-5-24 22:48
标题: VASP采用虚晶近似（VCA）计算gamma-Fe2O3能量为正值
各位老师好：
最近要计算一个gamma-Fe2O3相关的体系，该体系比较复杂，里面的Fe分为Fe-oct和Fe-tet，原胞中8个Fe-oct占据数为1，而剩下的13.33个Fe-oct占据了16个位置，平均每个位置的Fe-oct占据数为0.833。最开始尝试扩胞计算，但是体系大并且有磁性，结构和自洽收敛都比较困难。后来了解了一些虚晶近似相关的计算方法，决定试一下该方法，在INCAR增加了VCA = 0.833 1.0 1.0, POSCAR和POTCAR将Fe元素分为了2类，16个Fe-oct对应的权重为0.833，剩余8个Fe-tet对应权重1.00，O的权重为1.0，但是计算得到的能量却是正值，优化倒是也能正常进行，OUTCAR里面得到的NELECT似乎也和Fe21.33O32的化学计量式相对应，得到的结构也没有特别离谱，想咨询一下各位老师这样的结果是正常的还是哪里出了问题？
PS：同时附上了INCAR，POSCAR，POTCAR，OUTCAR以及XDATCAR文件。

作者
Author: xxxhhh 时间: 2024-5-25 11:28
您好，我想问一下您使用虚晶近似能算出来吗

作者
Author: I10140317 时间: 2024-5-25 14:08

xxxhhh 发表于 2024-5-25 11:28
您好，我想问一下您使用虚晶近似能算出来吗

附件就是优化完的，结构倒是肉眼看着不算离谱，但是能量都是正值，不知道是哪里有问题还是就这样。

作者
Author: xxxhhh 时间: 2024-5-25 14:17

I10140317 发表于 2024-5-25 14:08
附件就是优化完的，结构倒是肉眼看着不算离谱，但是能量都是正值，不知道是哪里有问题还是就这样。

我也是再算一个分数占据的体系，也是用了虚晶近似的方法，但优化着优化着，吸附分子就跑了，也不知道咋回事

作者
Author: dmacgyj 时间: 2024-10-17 21:34
同样的问题

作者
Author: zhengyr 时间: 2024-11-18 12:01
请问我VCA计算出来只跑一步然后能量都是NAN是怎么回事啊

欢迎光临计算化学公社 (http://bbs.keinsci.com/)