求助，用文献相同的参数计算晶格常数，结果差别很大，能否帮我看一下设置是否正确？

AustraliaSucks

采用MS2019的castep模块计算CoCrFeNi高熵合金的晶格参数，方法采用虚晶近似。参数采用文献“An investigation concerning generalized stacking fault behavior of AlCoxCrFeNi (0.25 ≤ x ≤ 2) high entropy alloys: Insights from first-principles study”的参数设置，具体参数如下：
“The Pardew-Burke-Ernzerhof (PBE) [48] in the class of the generalized gradient approximation (GGA) was used as the exchange-correlation functional. Notably, geometry optimization (minimum energy levels of atoms at different volumes) is a necessary condition for simulation. Therefore, the convergence parameters have been selected when the relaxation of structure was terminated. In this case, the convergence parameters were set as follows: total energy tolerance 10−5 eV/atom, force tolerance 0.03 eV/Å, maximum stress 0.05 GPa, and maximum displacement 0.001 Å. After adequate convergence tests, the Monkhorst-Pack [49] scheme k-points was set as 8 × 8 × 8 in the Brillouin zone. Moreover, ultrasoft pseudopotentials were applied and also the plane wave energy cut-off was eventually set as 560 eV.  ”
文献得到的晶格参数的结果为：
HEAs (in fcc phase)        Expt.                    GGA-PBE        REa (%)        GGA-PRBE        RE (%)        GGA-PW91         RE (%)        LDA        RE (%)
CoCrFeNi                         3.575 [52]        3.589        0.39                  3.601        0.72                3.548        0.75               3.484        2.54

我用了相同的参数进行了计算，参数设置如下图所示：


计算结果如下图所示：


我计算得到的结果远大于该文献计算的结果，能否请各位帮我看一下我的参数的设置是否哪里存在问题？非常感谢大家的帮助！

卡开发发

Co Cr Fe Ni可能有局域磁矩，为啥不考虑自旋极化？

AustraliaSucks

卡开发发 发表于 2022-1-25 17:08
Co Cr Fe Ni可能有局域磁矩，为啥不考虑自旋极化？

我极化其实也加了，结果也是差的很多。另外我也发邮件问过作者，他说他没加，他建议我改动的那些，我改动了还是不一样的结果

卡开发发

AustraliaSucks 发表于 2022-1-25 19:16
我极化其实也加了，结果也是差的很多。另外我也发邮件问过作者，他说他没加，他建议我改动的那些，我改动 ...

开自旋极化同时需要设置初始的formal spin。如果还不行，那你可以问问他有没有原始数据，能否核对一下计算参数。另外对称性等影响也应该尝试，原胞和单胞有时候算出来是可能存在一些差异。

AustraliaSucks

卡开发发 发表于 2022-1-25 19:41
开自旋极化同时需要设置初始的formal spin。如果还不行，那你可以问问他有没有原始数据，能否核对一下计 ...

问啦，他就说改精度就好了。我把截断能提高到1500eV左右，算出来的结果才差不多。但是MS ultrafine推荐的截断能才300多……当然如果能量收敛的话，我测试了一下2000eV多才差不多收敛。这种情况请问有遇到过吗？非常感谢

卡开发发

AustraliaSucks 发表于 2022-1-25 23:23
问啦，他就说改精度就好了。我把截断能提高到1500eV左右，算出来的结果才差不多。但是MS ultrafine推荐的 ...

那文章为啥说使用560eV就能行？MS推荐的值一般都偏低，按道理确实要做收敛性测试。不过2000还没收敛可能是赝势自身的问题，但是因为你用了VCA的方式，这个应该或多或少还是会有影响。

AustraliaSucks

卡开发发 发表于 2022-1-26 03:04
那文章为啥说使用560eV就能行？MS推荐的值一般都偏低，按道理确实要做收敛性测试。不过2000还没收敛可能 ...

我也母鸡了，所以我来求助了，我真的是搞得快崩溃了。

卡开发发

AustraliaSucks 发表于 2022-1-26 09:26
我也母鸡了，所以我来求助了，我真的是搞得快崩溃了。

自旋极化和初始磁矩还有原胞、单胞转换的问题你可以做些尝试，我现在手头没有空余服务器帮你一块尝试，有什么问题咱们可以一块讨论。