然后换了方法一就好了。记不得在哪查的，说原因是直接soc自洽是先猜一个电荷密度，这样难以收敛。
需要强调的是由于我们的体系简单,没有磁性不用设置初始磁矩;没有考虑SOC效应;没有考虑+U;且采用的是PBE-level的泛函,没有用meta-GGA或是HSE的泛函,且最重要的是初始结构足够合理(不会有奇怪的原子位置),因此程序很容易就优化成功了。正如前面实操3-2中的自洽计算收敛性问题,结构优化也容易出现不收敛的情况,即...
为了研究材料体系中的相互作用,我们需要研究一些假体系,例如,明明是有自旋轨道耦合(Spin-orbit coupling,SOC)相互作用的Bi2Se3我们会先不加SOC再加SOC,研究自旋轨道耦合对能带反转的作用。在有些体系中我们还会用到调节自旋轨道耦合强度的来研究自旋轨道耦合强度对体系性质的...
ISTART=1ICHARG=1要读取WAVCAR最好选取1而不用11，对结果影响是比较大的 LSORBIT=.TRUE.LMAXMIX=6计算soc的时候一定要有此参数，d电子4，f电子6MAGMOM=6*000400400400-400-400-400300300300300-300-300-300-372*0SAXIS=111磁矩的方向(见vasp手册，有两种设置方法)NBANDS=408能带数是线性计算的二倍ISYM=0计算soc...
过渡态计算的收敛准则只与力相关,即设置EDIFFG为合理的值即可,在前期试算过程中设置的值为0.02 eV/Å; NEB相关的关键词含义: LCLIMB=.TRUE. #通知VASP使用vtst计算; ICAHIN=0 #通知VASP使用vtst计算; IMAGES=7 #之前添加的图像数量N image number; ...
作者您好,我也出现了同样的问题,计算SOC时也会出现 WARNING: Sub-Space-Matrix is not hermitian in DAV 7 -1.459555082343374E-005之后加入IALGO=48后会出现警告WARNING in EDDRMM: call to ZHEGV failed,但是最后会收敛,请问这个警告会对结果有影响吗?或者您有没有更好的解决办法呢?期待您的回复,谢谢! 2022-10...
但是把CONTCAR拷贝成POSCAR,然后再计算能带时,就不收敛了,跑了NELM步之后自己结束。 我的INCAR如下: SYSTEM=0326 ENCUT=400 PREC=medium GGA=91 VOSKOWN=1 ISPIN=2 ICHARGE=11 POTIM=1 ISMEAR=0 LELF=.TRUE LORBIT=11 MAXMIN=30 ALGO=Fast NELM=120 ...
HSE+SOC ... 使用MS计算分波...精华评论 490105179 一共多少原子,运算多长时间? OhMyGod_2015 2楼: Originally posted by 490105179 at 2015-05-28 10:10:03 一共多少原子,运算多长时间? 一共18个原子,算了5分钟。请问就是这种不收敛的情况要怎么修改参数?新手,很多不懂。 490105179 3楼: Originally po...
LDAUTYPE = 1 LDAUL = 2 -1 LDAUU = 4 0 LDAUJ = 1 0 LORBIT=11 算得更准 19如果要用wannier算磁性体系。Spin要分开,这是应该用vasp+wannier+soc的版本跑,同时INCAR设置成spin=2的模式,但是LSORBIT不开。跑完会出来up和down,再分别对他们局域化。