需要说明的是我们不能仅仅通过改变VASP的INCAR输入参数来实现这种调节。要调优VASP中的自旋轨道耦合强度,需要修改源文件vasp_source_code_path/src/relativistic.F 并重新编译。如果要将强度减少一半,可以在relativistic.F的129行L●S项乘以0.5d0,如下:
DO I=0,1 DO J=0,1 DO M =1,2*LL+1 DO MP=1,2*LL+1 DLLMM(LMP+MP-1,LM+M-1,J+2*I+1)=DLLMM(LMP+MP-1,LM+M-1,J+2*I+1)+ & 0.5d0*SUM*LS(M,MP,I+2*J+1,LL) !!!line 129 relativistic.F file END DO END DO END DO END DO
注意,需要使用重新编译的vasp_ncl来运行计算。另外,最好比较0d0*L●S,1d0*L●S和vasp_std,vasp_ncl计算结果。 最后,请注意这是假的体系,只是用来辅助说明你研究体系中的物理规律。
No Comments
Leave a comment Cancel