1、從VASP在線手冊中整理的非線性磁矩和自旋軌道結合的計算說明如下所示。非線性磁矩計算：1 )計算基于非磁性基態(tài)的WAVECAR和CHGCAR文件。除了INCAR之外ISPIN=2ICHARG=1或11！ 導入WAVECAR和CHGCAR文件LNONCOLLINEAR=.TRUEMAGMOM=注意：在非線性磁矩的計算中，請在x、y、z方向上施加磁矩MAGMOM=1 0 00 1 0！ 表示第一個原子在x方向，第二個原子在y方向上有磁矩指定MAGMOM的值的前提是ICHARG=1(沒有wave car和chgccar文件)，或ICHARG=1或11 (有wave car和chgccar文件)，但是

2、上次的計算是非磁性的(ISPIN=1)。磁各向異性能(自旋軌道耦合)的計算：注意： lsorbit=.true.lnoncollinear=.true .選項自動打開，自旋軌跡計算僅適用于PAW人工因素，而不適用于超軟人工因素。自旋軌道耦合效應定義了初始磁矩的方向，這意味著能量取決于磁矩的方向，即存在磁各向異性能量(MAE )。 如下所示。LSORBIT=.TRUEsax is=s _ XSS _ yss _ z (quantisationaxisforsspin )默認值： sax is=(0，0，1 )，也就是說，x方向有正的無限小的磁矩，z方向有磁矩。有兩種方法可以使初始磁矩的方向與選定

3、方向平行MAGMOM=x y z！ 本地管理在x，y，zSAXIS=0 0 1！ quan tisation axis并行到z或者MAGMOM=0 0 total_magnetic_moment！ locallmationmentcmomentparalleltossaxis (請注意按原子指定)SAXIS=x y z！ quantisationaxisparalleltovector (x，y，z)(01等)兩種方法原則上應該相等，但實際上第二種方法更準確。 第二種方法可以讀取現有的線性或非磁性計算(WAVECAR )文件，并繼續(xù)另一個旋轉方向的計算(SAXIS值不變，MAGMOM不變)。

4、讀取非線性磁矩計算的WAVECAR時，將自旋方向指定為平行于SAXIS。計算磁各向異性的推薦步驟如下1 )首先計算線性磁矩，生成WAVECAR和CHGCAR文件(注意加入LMAXMIX )。添加到INCAR :LSORBIT=.TRUEICHARG=11！ 非自一致運行，讀chg car啊！ 或者，雖然ICHARG=1被最優(yōu)化為易磁化軸，但在這種情況下，應該提高EDIFF的精度LMAXMIX=4！ fordelementsincreasestmixto4，f: LMAXMIX=6啊！ 你需要的是lmaxmixalreadyinthecollinearcalculationSAXIS=x y z

5、！ magnetic字段的方向好像是01nbands=2* numberofbandsofcollinearrun！ 正則表達式nbandsoutcarISYM=0！ switchoffsymmetry (isym=0) whenspinorbitcouplingiseselectedGGA_COMPAT=.FALSE . itimprovesthenumericalprecisionofgafornoncollinearcalculationsLORBMOM=.TRUE . 計算軌道磁矩。如果繼續(xù)計算，VASP讀取WAVECAR和CHGCAR使自旋量化方向(磁場方向)與SAXIS方向平行。最

6、后，可以比較各方向磁矩的能量差異。注意：第二步原則上使用自協商計算(ICHARG=1)，但旋轉直到與SAXIS平行的磁場成為基態(tài)。 例如，如果平行于易磁化軸，這個過程很慢，能量的變化很小，所以如果收斂標準不嚴格，則自談判計算停止。注意：根據SAXIS的方向，VASP的輸入輸出的磁矩和類似自旋量包含在INCAR中MAGMOM行、OUTCAR和PROCAR文件的總磁矩和局部磁矩、WAVECAR的類似自旋軌道和CHGCAR的磁密度。mag mom-tag:http:/CMS.MPI.univie.AC.at/vasp/mag mom _ tag.html # incar-mag momlnoncol

7、linear :3358 CMS.MPI.univie.AC.at/vasp/lnoncollinear _ tag.htmllsorbit-tag http:/CMS.MPI.univie.AC.at/vasp/vasorbit _ tag.html問題：在第一次線性計算中得到WAVECAR和chgccar文件的必須是靜態(tài)計算的WAVECAR和chgccar文件嗎？ 能動態(tài)優(yōu)化嗎？ 需要在靜態(tài)計算中使用NUPDOWN鎖定磁矩嗎？進行非線性磁矩或自旋軌道耦合計算時，是否需要重新優(yōu)化結構？我目前的做法是在加入LMAXMIX=4結構優(yōu)化之后，使用LMAXMIX=4靜態(tài)計算(ICHARG=2，LWA

8、VE=.TRUE .LCHARG=.TRUE.)進行高收斂標準的靜態(tài)soc自協商計算，以考慮soc的影響，從而考慮空間spin-orbit請求幫助！最近用spin-orbit計算系統(tǒng)磁總是錯誤的，揭示了計算中的幾個問題和疑問，尋求蝦的答案與使用spin-orbit的spin-polarizedINCAR文件有什么區(qū)別？ 除了需要追加的情況LSORBIT=.TRUESAXIS=X Y ZNBANDS是共用線的兩倍，MAGMON也在xyz方向以外沒有必要追加。GGA_COMPAT=.FALSE。ISYM=0怎么了？ (這是在VASP手冊中看到的。)2 .自旋量子軸SAXIS=X Y Z應該如何設定

9、？ 應該與系統(tǒng)的網格結構有關，軸的設定應該基于實空間網格來設定，還是基于逆空間網格來設定，請舉例說明一下根據VASP手冊所提供的方法，在計算spin-orbit之前，應該在計算公共線的WAVECAR和CHGCAR之后開始spin-orbit的計算。但是，如果我在計算時引進以前計算過的WAVECAR和CHGCAR的話，計算會錯誤，是否應該計算出不知道是什么原因的共同線的WAVECAR和CHGCAR，來計算spin-orbit呢？計算spin-orbit不需要超細胞嗎？ (如果不摻雜或摻雜)如何分析使用超細胞的計算結果呢？對假的勢頭(例如GGA-PBE等)有要求，對算法(例如ALGO=Normal

10、，ALGO=Damped等)有要求，對計算精度等有要求嗎得出答案1 .在1. ISYM=0為on時相加，在用GGA交換泛函時用GGA_COMPAT=.FALSE .改善了精度。2. SAXIS和MAGMOM有兩個共同設定，如手動介紹的那樣，默認為01方向，是相對的實空間正交坐標系。 一個是SAXIX為01方向，此時MAGMOM為任意方向，另一個是MAGMOM為0z或0-z方向，SAXIX為任意方向，另一個可以讀取以前存在的共用線或非共用線的電荷密度和波函數，精度很高。3 .導入前計算出的WAVECAR和CHGCAR的錯誤原因是，KPOINTS將SO設為on時對稱性發(fā)生了變化。 使用SO生成的I

11、BZKPT，復制到KPOINTS，重新計算WAVECAR和CHGCAR，讀出即可。4 .不需要超細胞！5.SO只支持PAW假貨勢，精度要求高，能量精度EDIFF一般為10-7eV2)SOC版本：cp化妝文件. MPI化妝文件. SOC在makefile.SOC中更改CPP=$ (CPP _ )-dmpi-dhost=Linux IFC-difc- dcache _ size=5000-dpgf 90-davoid alloc-dngzhalf - dmpi _ block=-duse _ collect ve-dscalapack -DRPROMU_DGEMV -DRACCMU_DGEMV中刪除-DNGZhalf然后make-f化妝文件. SOC得到了vasp，得到了mvvasvaspvasp.MPI.SOC.nebMAE (磁各向異性能) -非共線磁矩計算系統(tǒng)=Fe/gral real=自動ALGO=FastIALGO=48ISYM=0ISTART=1ICHARG=11ENCUT=500NPAR=2ISMEAR=0； SIGMA=0.2GGA=91； VOSKOWN=1GGA_COMPAT