在考虑SOC 后需要关闭对称性,也就是设置 ISYM =0 ,但此时很容易碰到一个报错:WARNING: Sub-Space-Matrix is not hermitian in DAV。下面是我的经验: 方法一: 增加LMAXMIX=4 。 方法二: 增加LMAXMIN=4 和 GGA_COMPAT = False GGA_COMPAT 这个参数很重要,一般容易忽略。 官网解释:We recommend setting GGA_...
通常,在使用自旋-轨道耦合时要格外小心:能量差很小,k点收敛很慢,而且所需的计算机时间可能是无限的。 建议将GGA_COMPAT设置为FALSE。因为这提高了GGA计算的数值精度。 参考:Category:Noncollinear magnetism - Vaspwiki 编辑于 2022-11-27 10:12・IP 属地江苏 ...
2、一般来说,使用自旋-轨道耦合时要非常小心:能量变化和差异很小,K点收敛繁琐而缓慢,因此所需要的机时是不确定的。 3、建议在VASP.4.6中对非共线计算设置GGA_COMPAT = .FALSE.,因为这样可以提高GGA计算的数值精度。 在计算成功后,VASP将以下结果写入OUTCAR中:Spin-Orbit-Coupling matrix elements Ion:1E_soc:...
GGA_COMPAT通常需要,应该被设置,因为磁各向异性能量通常需要精确到亚meV数量级。 当计算自旋轨道耦合,特别是磁各向异性时通常需要非常小心:能量差异非常小,k点的收敛冗长而且缓慢,需要耗费大量的计算时间。此外,这一特征--尽管长期存在于VASP中--在最新的版本中依然存在,你可以尝试频繁地升级发现这一点。不敢保证,你...
3、建议在VASP.4.6中对非共线计算设置GGA_COMPAT = .FALSE.,因为这样可以提高GGA计算的数值精度。 在计算成功后,VASP将以下结果写入OUTCAR中:Spin-Orbit-Coupling matrix elements Ion:1E_soc:-0.0984080l=10.0000000-0.0134381-0.0134381-0.01343810.0000000-0.0134381-0.0134381-0.01343810.0000000l=20.0000000-0.00050720...
GGA_COMPAT=.FALSE. ! it improves the numerical precision of GGA for non collinear calculations LORBMOM=.TRUE. !计算轨道磁矩 继续计算,VASP会读取WA.VECAR 和 CHGCAR将自旋量子化方向(磁场方向)平行于SAXIS方向。 最后可以比较各个方向磁矩时能量的不同。
SOC得到了vasp,得到了mvvasvaspvasp.MPI.SOC.nebMAE (磁各向异性能) -非共线磁矩计算系统=Fe/gral real=自动ALGO=FastIALGO=48ISYM=0ISTART=1ICHARG=11ENCUT=500NPAR=2ISMEAR=0; SIGMA=0.2GGA=91; VOSKOWN=1GGA_COMPAT=.FALSE。ISPIN=2#MAGMOM=1*5 2*0 1*4LORBIT=11LNONCOLLINEAR=.TRUELSORBIT=....
3、建议在VASP.4.6中对非共线计算设置GGA_COMPAT = .FALSE.,因为这样可以提高GGA计算的数值精度。 在计算成功后,VASP将以下结果写入OUTCAR中:Spin-Orbit-Coupling matrix elements Ion: 1 E_soc: -0.0984080 l= 1 0.0000000 -0.0134381 -0.0134381 -0.0134381 0.0000000 -0.0134381 -0.0134381 -0.0134381 ...
GGA_COMPAT=.FALSE.!在梯度场中应用球面截断能 VASP读入WAVECAR和CHGCAR文件,将自旋量子轴对齐SAXIS矢量,这意味着现在磁场平行于SAXIS矢量,执行非线性计算。通过比较不同方向的能量,可以确定磁各向异性。请记住,原则上,在VASP中一个完全地自洽计算(ICHARG=1)也是有可能的,但是这种状况将会允许自旋波函数从它们的初始值...
(31) LSORBIT:自旋轨道耦合计算,默认为False. 如果设置LSORBIT=.TRUE.,则VASP会自动设置LNONCOLLINEAR=.TRUE.同时注意设置MAG-MOM SAXIS和GGA COMPAT(=False)等参数。 (32) I_CONSTRAINED_M:通过引入惩罚项以限制磁矩方向和大小。I_CONSTRAINED_M=1只限制磁矩方向,I_CONSTRAINED_M=2同时限制磁矩大小和方向。