24小时热门版块排行榜    

查看: 2763  |  回复: 8

weiwei128

金虫 (小有名气)

[交流] 【求助】vasp磁矩计算 已有4人参与

我最近用vasp计算体系的态密度和磁矩时遇到了一些问题, 希望有大侠们指点,谢谢.
   1. 体系优化之后, 态密度和磁矩的计算可以同时进行吗? 即设置
          Istart=0,
          ISpin=2,
          magtom= 52*2.0,
          Icharge=2,
          IBrion=-1,
          NSW=0,
          Lorbit=11
       同时删除了以前产生的文件. 我是看以前的一些handson例子有这样算的,不知道行不行?
  2. 如果态密度和磁矩不能同时算的话, 那应该怎么办? 采用自洽来算态密度?然后磁矩?....

希望大家指点指点,感激不尽!
回复此楼

» 收录本帖的淘帖专辑推荐

vasp计算资料

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)

weiwei128(金币+1):谢谢 2010-05-12 16:31:52
在INCAR中设置RWIGS后,磁矩就会在OUCAR中输出,在体系优化的时候就可以得到磁矩,不需要等到算态密度的时候
如果是在并行机上,还需要加入NPAR=1
2楼2010-05-12 15:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

★ ★
zzy870720z(金币+2):感谢交流 2010-05-12 16:27:36
weiwei128(金币+5):非常感谢 2010-05-12 16:32:56
引用回帖:
Originally posted by weiwei128 at 2010-05-12 15:36:24:
我最近用vasp计算体系的态密度和磁矩时遇到了一些问题, 希望有大侠们指点,谢谢.
   1. 体系优化之后, 态密度和磁矩的计算可以同时进行吗? 即设置
          Istart=0,
          ISpin=2,
          magtom ...

计算态密度时就会产生磁矩,也就是一般的静态计算得到体系的基态以后,态密度,磁矩(或者称自旋密度)都已经产生了,如何得到和你采用什么赝势和关联有关。
如果你是采用PAW方法,那么直接在做静态计算时在INCAR里面添加LORBIT=11或者12就可以在静态计算结束时往OUTCAR里面输出磁矩,在最后的magnetization那一段就是。
如果采用的是USPP方法,好像应该在INCAR里面设置RWIGS什么的,我没有用过,具体的不太清楚如何设置。
3楼2010-05-12 16:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by quby7503 at 2010-05-12 15:44:04:
在INCAR中设置RWIGS后,磁矩就会在OUCAR中输出,在体系优化的时候就可以得到磁矩,不需要等到算态密度的时候
如果是在并行机上,还需要加入NPAR=1

NPAR值的设置好像是关于如何分配并行任务的,和体系的最终计算结果没有什么关系。
4楼2010-05-12 16:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiwei128

金虫 (小有名气)

引用回帖:
Originally posted by quby7503 at 2010-05-12 15:44:04:
在INCAR中设置RWIGS后,磁矩就会在OUCAR中输出,在体系优化的时候就可以得到磁矩,不需要等到算态密度的时候
如果是在并行机上,还需要加入NPAR=1

NPAR是并行时分配任务时的节点设置吧? 应该不会影响磁矩和态密度结果的
5楼2010-05-12 16:38:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiwei128

金虫 (小有名气)

引用回帖:
Originally posted by nkleof at 2010-05-12 16:25:18:

计算态密度时就会产生磁矩,也就是一般的静态计算得到体系的基态以后,态密度,磁矩(或者称自旋密度)都已经产生了,如何得到和你采用什么赝势和关联有关。
如果你是采用PAW方法,那么直接在做静态计算时在IN ...

我用的是PAW-GGA,那这样的话我的设置应该是没问题吧?
6楼2010-05-12 16:39:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)

weiwei128(金币+1):谢谢回复 我们用的是NPAR=4, 也设置了Lorbit=11 2010-05-12 17:20:00
VASP4.5   NPAR不等于1    no site projected DOS
VASP4.6   NPAR不等于1, LORBIT=0-5    no site projected DOS
7楼2010-05-12 16:59:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
8楼2010-05-12 18:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

★ ★
weiwei128(金币+1):是的 应该是大写 这里没注意 谢谢 2010-05-12 19:19:06
zzy870720z(金币+2):感谢交流 2010-05-13 08:27:03
引用回帖:
Originally posted by weiwei128 at 2010-05-12 16:39:22:


我用的是PAW-GGA,那这样的话我的设置应该是没问题吧?

你的INCAR里面的tag都很奇怪啊,我记得都应该是大写的,而且有的标签你写错了,比如
ICHARG     不应该写成ICHARGE
另外如果是做自洽计算的话,ISMEAR是个重要的设置,根据体系不同要设定不同的值,金属一般选1或者2,非金属选5.
9楼2010-05-12 18:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weiwei128 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见