[求助]vasp如何計(jì)算磁各向異性能?
請教一個(gè)vasp問題:
vasp計(jì)算磁各向異性能的時(shí)候,MAGMOM和SAXIS應(yīng)該如何設(shè)置呢?
比如一維Co單原子鏈 :
------Co-------Co-------Co-------Co------- ----->Z
鏈的軸沿向Z方向,每個(gè)元胞里一個(gè)Co原子,單個(gè)Co原子磁矩約2ub
現(xiàn)在計(jì)算磁矩沿z方向和x方向兩種情況的磁各向異性能,MAGMOM和SAXIS這兩項(xiàng)應(yīng)該怎么設(shè)置呢?
MAGMOM的三個(gè)分量根SAXIS有什么關(guān)系,它們倆跟POSCAR里晶胞的基矢又有什么關(guān)系?
我算出來的磁矩總是沿著z方向 (MAGMON只有Z分量有值,x和y分量都為零).
[ Last edited by wuchenwf on 2009-6-21 at 20:26 ]
返回小木蟲查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號(hào)
控制磁矩方向 要用SAXIS
你算的結(jié)果就是因?yàn)闆]你沒控制住方向
不過偶還是很久前做過這個(gè)測試，所以提供不了詳細(xì)的幫助^_^
自己多看看vasp的相關(guān)計(jì)算的說明書，還是有一些用處
我把幾種可能性都嘗試了一下:
----------------------------------------------------------------------------------------------------------
SAXIS MAGMOM Energy
0 0 1 0 0 2 -43.344340 eV
0 0 1 2 0 0 -43.344340 eV
1 0 0 0 0 2 -43.289052 eV
1 0 0 2 0 0 -43.289052 eV
----------------------------------------------------------------------------------------------------------
也就是說最后的能量只跟SAXIS有關(guān)系, 而跟MAGMON的設(shè)置沒有關(guān)系, 而且最后的OUTCAR里面磁矩只有z方向有分量, 而x,y方向都沒有分量(上面四種情況都是這樣的, 即使我初始設(shè)置了MAGMOM磁矩只有x方向有磁矩), 為什么呢
，
控制磁矩方向 要用SAXIS，這個(gè)我上面說了，在你的計(jì)算情況里面也是對應(yīng)的。
至于磁矩方向，你說的情況有點(diǎn)奇怪。我這里上傳一個(gè)圖片，你看看，是vasp workshop第4部分上的，你可以去官方網(wǎng)下載看看
這里是通過magmom控制磁矩，不過實(shí)際上計(jì)算的的時(shí)候，控制磁矩方向 要用SAXIS
關(guān)于磁矩分量，如果你是一個(gè)原子的話，那么z和x等向的話就沒事情了。如果是一個(gè)dimer，那么就有區(qū)別。你自己仔細(xì)看一下。這個(gè)計(jì)算我很久前做過的一個(gè)測試，也不是很熟悉呢
上面那個(gè)圖是noncollinear的情況,noncollinear計(jì)算不需要設(shè)置SAXIS的,和磁各向異性能計(jì)算好像不一樣的,手冊,workshop都看了,還是不得其解,我猜測是不是MAGMOM是相對于SAXIS坐標(biāo)的,就是磁矩方向總是沿著SAXIS方向?
十分感謝csfn,今天看了一下手冊,MAGMOM確實(shí)是相對于SAXIS坐標(biāo)的
不過具體之間的關(guān)系還沒搞清楚,
thank you!