24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4510  |  回复: 16

ghzhong

金虫 (小有名气)

[交流] 【讨论】非共线磁性参数设置 已有4人参与

请教虫友:

我想做一个体系的非共线磁性(磁性原子是zig-zag链状排列),看看磁结构怎么样,比如:每个原子在xyz方向上的投影磁矩(含方向+-)、次各向异性等。这时我在INCAR中只设置非共线参数开关能得到这些磁性信息吗?

在这种计算中,INCAR中自旋轨道耦合参数开关不加可以吗?我总是认为,考虑自旋轨道耦合后,加的磁场又把原子磁矩方向固定到磁场方向了,怎么还能让原子进行非共线磁场收敛呢。

请教有经验的VASP虫友,谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wo8165827

金虫 (小有名气)


ghzhong(金币+2):谢谢参与
ghzhong(金币+1): 2010-05-21 02:14:27
1:计算非共性磁性计算,首先要重新编译下VASP。将Makefile里的-DNGXhalf 和 -DNGZhalf去掉。

2: add the following lines to INCAR
LNONCOLLINEAR=.True.
LSORBIT=.True.
SAXIS = # please give the spin quantization axis here, like 0 0 1 for the z-axis)
MAGMOM= # please give a triplet of numbers for each atom here, and please have a look at the manual (chapter non-collinear calculations and spin-orbit tag) on how the direction of the magnetic moments has to be defined with respect to the spin-quantization axis)
LORBMOM=.True.
ISYM= -1

» 本帖已获得的红花(最新10朵)

计算的魅力
5楼2010-03-15 15:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ghzhong

金虫 (小有名气)

顶上来,别沉了
2楼2009-12-06 10:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlshi

金虫 (正式写手)

★ ★
ghzhong(金币+2):谢谢参与
qasd(金币+1,VIP+0):xiexie~ 12-7 19:55
ghzhong(金币+1): 2010-05-21 02:13:41
非共线计算默认 自旋轨道耦合参数开 ,可以根据说明书上的计算一下,看看能不能得到相应的结果
3楼2009-12-06 15:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chgzhong

木虫 (正式写手)

★ ★
ghzhong(金币+2):谢谢参与
qasd(金币+1,VIP+0):xiexie~ 12-7 19:55
ghzhong(金币+1): 2010-05-21 02:13:50
非共性计算可直接按照说明书上的介绍设置参数,加上自旋轨道耦合后计算的结果与各向异性场相关!
4楼2009-12-06 15:29:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟


ghzhong(金币+2):谢谢参与
ghzhong(金币+1): 2010-05-21 02:14:33
比较复杂,现在共线的还没算好,帮支持下
6楼2010-03-15 15:19:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chimegreen

木虫 (正式写手)

Dr.

★ ★ ★
ghzhong(金币+2):谢谢参与
ice_rain(金币+2):谢谢分享 2010-03-15 22:21
ghzhong(金币+1): 2010-05-21 02:14:55
请教一个vasp问题:
vasp计算磁各向异性能的时候,MAGMOM和SAXIS应该如何设置呢?

比如一维Co单原子链 :
------Co-------Co-------Co-------Co-------  ----->Z
链的轴沿向Z方向,每个元胞里一个Co原子,单个Co原子磁矩约2ub

现在计算磁矩沿z方向和x方向两种情况的磁各向异性能,MAGMOM和SAXIS这两项应该怎么设置呢?

MAGMOM的三个分量根SAXIS有什么关系,它们俩跟POSCAR里晶胞的基矢又有什么关系?

我算出来的磁矩总是沿着z方向 (MAGMON只有Z分量有值,x和y分量都为零).  
相关回复:

作者: csfn   发布日期: 2008-06-08

控制磁矩方向 要用SAXIS

你算的结果就是因为没你没控制住方向

不过偶还是很久前做过这个测试,所以提供不了详细的帮助^_^
自己多看看vasp的相关计算的说明书,还是有一些用处

作者: cationly   发布日期: 2008-06-09

QUOTE:
Originally posted by csfn at 2008-6-8 21:20:
控制磁矩方向 要用SAXIS

你算的结果就是因为没你没控制住方向

不过偶还是很久前做过这个测试,所以提供不了详细的帮助^_^
自己多看看vasp的相关计算的说明书,还是有一些用处
我把几种可能性都尝试了一下:
----------------------------------------------------------------------------------------------------------
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有关系, 而跟MAGMON的设置没有关系, 而且最后的OUTCAR里面磁矩只有z方向有分量,
而x,y方向都没有分量(上面四种情况都是这样的, 即使我初始设置了MAGMOM磁矩只有x方向有磁矩), 为什么呢?

作者: csfn   发布日期: 2008-06-09

控制磁矩方向 要用SAXIS,这个我上面说了,在你的计算情况里面也是对应的。

至于磁矩方向,你说的情况有点奇怪。我这里上传一个图片,你看看,是vasp workshop第4部分上的,你可以去官方网下载看看
未命名.JPG

作者: csfn   发布日期: 2008-06-09

这里是通过magmom控制磁矩,不过实际上计算的的时候,控制磁矩方向 要用SAXIS

关于磁矩分量,如果你是一个原子的话,那么z和x等向的话就没事情了。如果是一个dimer,那么就有区别。你自己仔细看一下。这个计算我很久前做过的一个测试,也不是很熟悉呢

作者: cationly   发布日期: 2008-06-09

QUOTE:
Originally posted by csfn at 2008-6-9 09:16:
这里是通过magmom控制磁矩,不过实际上计算的的时候,控制磁矩方向 要用SAXIS

关于磁矩分量,如果你是一个原子的话,那么z和x等向的话就没事情了。如果是一个dimer,那么就有区别。你自己仔细看一下。这个计算 ...
上面那个图是noncollinear的情况,noncollinear计算不需要设置SAXIS的,和磁各向异性能计算好像不一样的,手册,workshop都看了,还是不得其解,我猜测是不是MAGMOM是相对于SAXIS坐标的,就是磁矩方向总是沿着SAXIS方向?

作者: cationly   发布日期: 2008-06-09

十分感谢csfn,今天看了一下手册,MAGMOM确实是相对于SAXIS坐标的
不过具体之间的关系还没搞清楚,
thank you!
快乐科研,快乐生活。
7楼2010-03-15 15:29:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chimegreen

木虫 (正式写手)

Dr.

ghzhong(金币+1): 2010-05-21 02:15:12
一个论坛上很久以前的帖子,希望有所帮助
快乐科研,快乐生活。
8楼2010-03-15 15:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuling985

银虫 (小有名气)


ghzhong(金币+2):谢谢参与
ghzhong(金币+1): 2010-05-21 02:15:16
虽然注册时间短,爱上小木虫了
9楼2010-04-22 18:27:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)


ghzhong(金币+2):谢谢参与
ghzhong(金币+1): 2010-05-21 02:15:23
支持VASP
大家一起来读代码,解剖它
尽早国产化
10楼2010-04-22 18:49:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghzhong 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +20 让我上岸吧阿西 2026-04-09 20/1000 2026-04-15 14:40 by 黑科技花岗岩
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26 慕绝cc 2026-04-09 30/1500 2026-04-14 18:50 by 蔡苏阳
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
信息提示
请填处理意见