vasp使用VCA弛豫结构原子受力收敛问题

Shana

各位老师好，
最近尝试在vasp中使用VCA计算固溶体性质，在INCAR中设置的原子受力收敛精度为0.01 eV/A，vasp在结构弛豫reached required accuracy之后，小弟查看OUTCAR，发现TOTAL-FORCE却远大于收敛精度。小弟想咨询两个问题：
1. OUTCAR中的TOTAL-FORCE为什么在大于收敛精度的情况下，vasp仍然会reached required accuracy？
2. 这种情形下的结构，能否认为达到了收敛精度，可以进行进一步的性质计算？

卡开发发

VASP当中每个原子的力F=(Fx,Fy,Fz)，|F|=sqrt(Fx^2+Fy^2+Fz^2)，EDIFFG收敛标准是按照原子最大力|F|max的情况判断的，这个在其他帖子有相关讨论。

Shana

卡开发发 发表于 2023-9-26 11:02
VASP当中每个原子的力F=(Fx,Fy,Fz)，|F|=sqrt(Fx^2+Fy^2+Fz^2)，EDIFFG收敛标准是按照原子最大力|F|max的情 ...

卡老师好。|F|肯定是大于|Fx|的，以我的OUTCAR中的受力来看，第一个原子的受力|Fx|已经远大于EDIFFG了，但是vasp仍然认为达到了收敛标准。这是因为VCA的原因吗？

卡开发发

Shana 发表于 2023-9-26 12:53
卡老师好。|F|肯定是大于|Fx|的，以我的OUTCAR中的受力来看，第一个原子的受力|Fx|已经远大于EDIFFG了， ...

你可以传一下完整文件，按道理不会在没达到标准而给出达到标准的提示。

Shana

卡开发发 发表于 2023-9-26 13:16
你可以传一下完整文件，按道理不会在没达到标准而给出达到标准的提示。

卡老师，这是我做VCA的INCAR，POSCAR和OUTCAR

卡开发发

Shana 发表于 2023-9-26 13:22
卡老师，这是我做VCA的INCAR，POSCAR和OUTCAR

VCA的力判据不是按照POSCAR当中每个原子的。

Shana

卡开发发 发表于 2023-9-26 13:48
VCA的力判据不是按照POSCAR当中每个原子的。

好的，谢谢卡老师。现在敢放心的继续去算了

卡开发发

Shana 发表于 2023-9-26 13:50
好的，谢谢卡老师。现在敢放心的继续去算了

问题应该不大，只是VCA自身可靠性需要你自己去把握。你的VASP有VASP-VTST的插件，你也可以参照你的OUTCAR中FORCES: max atom, RMS字段。VCA下力的处理可以参考源码的force.F中的VCA_FORCE这个subroutine，和一般非VCA的情况是有所不同的。

Shana

卡开发发 发表于 2023-9-26 13:48
VCA的力判据不是按照POSCAR当中每个原子的。

卡老师，再问您一个问题。小弟最近也在关注磁性的计算，以及居里温度的模拟，然后在调研之后，遇到了两个问题。
一是，关于磁性结构的计算问题。目前看到的文章在计算磁性结构方面有一个很大的区别。一类文章在计算磁性结构的时候，会以铁磁态、反铁磁态进行结构弛豫，进而得到最低能量态及其对应的晶格常数；另一类文章会先进行弛豫（这个弛豫好像是无磁性的），然后静态自洽计算铁磁态、反铁磁态，找到最低的磁性态。第二类方法我遇到的更多，虽然说磁性对结构的影响可能会比较弱。但是我更倾向于第一类方法。
二是，关于海森堡模型的问题。海森堡模型在计算磁交换能和磁晶各项异能的时候，需要考虑铁磁和反铁磁的能量差。如果采用第一类方法计算磁性结构的话，铁磁态的结构和反铁磁结构是不同的，那么这种方法计算的交换能和各向异能可能是存在问题的。第二类方法，固定了结构，计算出的交换能和各向异能可能更复合模型一点。
小弟现在处于迷惑之中。无法区分这两类方法的优劣。

卡开发发

Shana 发表于 2023-9-26 14:00
卡老师，再问您一个问题。小弟最近也在关注磁性的计算，以及居里温度的模拟，然后在调研之后，遇到了两个 ...

磁性问题我了解有限，仅作为讨论：
1、只讨论不同磁性所对应的稳定结构，应该各自优化比较合适。
2、对Heissenberg模型来说，如果基态是铁磁性，然后需要通过计算这个结构的几种反铁磁的状态来联立解耦合常数J，按我理解此时应当保持在同一结构上，最后再利用平均场或Monte Carlo去估计Curie温度之类的性质。

Shana

卡开发发 发表于 2023-9-26 14:22
磁性问题我了解有限，仅作为讨论：
1、只讨论不同磁性所对应的稳定结构，应该各自优化比较合适。
2、对 ...

好的，谢谢卡老师。有了新的理解了