求助：在磁性物體表面放置吸附物之計算

shihchengli

求助各位老師：
我目前在計算PET單體在ZnFe2O4(110)之吸附能，ZnFe2O4根據文獻報導具有反磁性，我已經進行Spin-polarized DFT+U計算得到合理的晶體結構，計算之晶體常數為a=b=8.438 Å及c=8.426 Å，與文獻報導之實驗數據a=b=c=8.451 Å吻合。後續我使用該結構建構2×3×1 ZnFe2O4(110)表面，且添加了20 Å的真空層，該slab模型具有7層，我將下面三層鎖住並進行結構優化。
            
在取得優化完畢的slab model，我進一步在該表面上方放置PET單體，並進行結構優化。
然而我卻發現，在吸附物模型算第一步的SCF energy時，有部分原子的磁矩與原先沒放吸附物的slab model相差很多，想在這邊請教各位老師有沒有什麼方法可以解決這個問題？

INCAR設置如下(粗體字為相對於slab模型，在吸附物模型中額外添加的參數設定)：
SYSTEM = ZnFe2O4
ENCUT = 520.000000
EDIFF = 1.00e-05
EDIFFG = -0.02
GGA = PE
ISMEAR = 0
SIGMA = 0.05
PREC = Accurate
NELMIN = 4
NELM = 400
NSW = 200
IBRION = 2
ISIF = 2
LAECHG = .FALSE.
LWAVE = .TRUE.
LREAL = Auto
NPAR = 4
ISPIN = 2
AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions
AMIX MAG = 0.8
BMIX MAG = 0.0001 #almost zero, but 0 will crash some versions
LORBIT = 11
MAGMOM = 32*0 32*6.4 24*-6.4 112*0 24*0
IVDW = 11
LDAU      = .TRUE.
LDAUTYPE  = 2
LDAUL     = 2 2 -1
LDAUU     = 4.70 4.50 0.00 0.00 0.00
LDAUJ     = 0.00 0.00 0.00 0.00 0.00
LDAUPRINT = 2
LMAXMIX   = 6
LDIPOL = .TRUE.
IDIPOL = 3