24小时热门版块排行榜    

查看: 1737  |  回复: 4

leila3409

捐助贵宾 (正式写手)

[求助] VASP 孤立原子 求磁矩

我用VASP想算bcc Fe磁矩是没有问题的,2.32 miu_B,几乎接近于实验值,也与文献上其他结果相同。对于一个孤立的Fe原子,理论值应该是4.0 miu_B的。可是我无论如何得不到这个值,设置同bcc的INCAR,这样的设置有问题吗?求高手!!!
  SYSTEM =  NI
   ISTART =  0;
   ISPIN  =  2;
   EDIFF =  1E-06;
   ISMEAR =  1;
   ICHARG = 2;
   NELM = 1000;
   IBRION = 2;
   EDIFFG = 1E-05;
   ISIF = 3;
   ALGO  = fast;
   LREAL  =  F;
   NPAR = 1;
   LORBIT = 1;
   WOSKOWN = 1;
   NSIM = 4;
   PREC = high
   NSW=0
   ENCUT=350
   RWIGS =1.39115
回复此楼

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

VASP

» 猜你喜欢

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

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

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖


fzx2008(金币+1): 谢谢回帖交流! 2011-11-26 12:30:55
leila3409(金币+25): 2011-11-26 14:27:03
关键是你有没有孤立掉这个粒子,孤立要加很大的真空层
2楼2011-11-26 10:56:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)

【答案】应助回帖

leila3409(金币+25): 2011-11-26 14:27:10
你是怎么构建一个一个孤立的Fe原子的?
paper我要好多好的的paper
3楼2011-11-26 14:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)

引用回帖:
2楼: Originally posted by liliangfang at 2011-11-26 10:56:23:
关键是你有没有孤立掉这个粒子,孤立要加很大的真空层

使用20*20*20的超晶胞。孤立原子坐标设为(0,0,0)即可。发现问题了,超晶胞变得很大时,系统的总磁矩与每个原子的局域磁矩有很大差别,这时只要取总磁矩/原子数就行了。

我对超晶胞变大,OUTCAR里magnetization与out里的mag差别变大还是不太理解。请指教!
4楼2011-11-26 14:30:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)

引用回帖:
3楼: Originally posted by liuhuimeiwoo at 2011-11-26 14:23:00:
你是怎么构建一个一个孤立的Fe原子的?

使用20*20*20的超晶胞。孤立原子坐标设为(0,0,0)即可。发现问题了,超晶胞变得很大时,系统的总磁矩与每个原子的局域磁矩有很大差别,这时只要取总磁矩/原子数就行了。

我对超晶胞变大,OUTCAR里magnetization与out里的mag差别变大还是不太理解。请指教!
5楼2011-11-26 14:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leila3409 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:42 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:30 by cITaGg3p5edr
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-18 03:18 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-18 01:54 by cITaGg3p5edr
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-17 23:05 by cuPZTDXS3VOd
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 22:53 by cuPZTDXS3VOd
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 15:29 by 0ngsZVSCTS7z
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:27 by ObVzqQQh1rrE
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:16 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 02:14 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-17 01:16 by Tql5LQhh5rLK
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见