应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp 计算磁矩结果比别人的要大
91/1返回列表
查看: 1604  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

林子6682

新虫 (初入文坛)

[交流] 【求助】vasp 计算磁矩结果比别人的要大

最近看了两篇关于Co2FeAl0.5Si0.5的文献（J. Phys: Condens. Matter 19(2007) 326216; J. Phys. D: Appl. Phys. 40(2007)1582-1586)，觉得比较有兴趣，便做了一下Co2FeAl0.5Si0.5这种四元合金的计算。结果我发现我得到磁矩为5.75，而这两篇文献的值均为5.5，经过多次计算我得到值都是5.75，我想请教大家究竟是怎么回事。有知道的请指点一下，谢谢！
以下是我的incar文件：
Startparameter for this Run:
   SYSTEM = L21 Co2FeAlSi-static-calculation
   ENCUT  = 500
   NWRITE = 2    write-flag & timer
   PREC   = Accurate
   ISTART = 0   job   : 0-new  1-cont  2-samecut
   ICHARG = 2    charge: 1-file 2-atom 10-const
   ISPIN  = 2    spin polarized calculation?
   MAGMOM = 8*3.0 4*4.0 2*0.0 2*0.0
   INIWAV = 1    electr: 0-lowe 1-rand  2-diag
   NELM   = 60
   NELMIN = 6; NELMDL= -5     # of ELM steps
   EDIFF  = 1E-04   stopping-criterion for ELM
   VOSKOWN = 1
   NPAR   = 1
   LORBIT = 11
   RWIGS  = 1.16 1.17 1.18 1.11
   LDAU   = .TRUE.
   LDAUTYPE = 4
   LDAUL = 2 2 1 1
   LDAUU  = 4.453 4.185 0 0
   LDAUJ  = 0.655 0.625 0 0

Electronic Relaxation 2
   IALGO  = 38    algorithm
   LDIAG  = T    sub-space diagonalisation
   LREAL  = Auto     real-space projection
   ROPT   = 2E-4 2E-4 2E-4 2E-4
   ALGO   = Fast
DOS related values:
   ISMEAR = -5
   SIGMA = 0.1  broadening in eV -4-tet -1-fermi 0
   NEDOS = 1000

Ionic Relaxation
   EDIFFG = 1E-03   stopping-criterion for IOM
   NSW    = 0   number of steps for IOM
   ISYM   = 2    0-nonsym 1-usesym 2-use
回复此楼

» 猜你喜欢
1楼 2009-10-29 14:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1,VIP+0):感谢参与 12-18 09:39
这个要看你用的参数和文献的是否一致
赞一下(11人) 回复此楼
相逢是缘！
2楼2009-10-29 15:27:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
不好意思，我也不知道为什么，我想请教下，就是这里面关于Co和Fe的U和J值你是怎么得到的,谢谢
赞一下(1人) 回复此楼
3楼2009-10-29 17:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+2,VIP+0):谢谢jghe专家解答 10-29 21:48
看你用的是什么泛函，什么赝势，还有 ISMEAR SIGMA 这个对磁矩影响很大
你得总磁矩从OUTCAR 还是 OSZICAR里得到的？
赞一下(7人) 回复此楼
4楼2009-10-29 20:27:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
引用回帖:
Originally posted by 含草 at 2009-10-29 17:14:
不好意思，我也不知道为什么，我想请教下，就是这里面关于Co和Fe的U和J值你是怎么得到的,谢谢

一般看文献，别人用多少就用多少
赞一下 回复此楼
5楼2009-10-29 20:27:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemix-chen

银虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1,VIP+0):谢谢参与 10-29 21:48
个人觉得可能是你的赝势选取的问题，你可以换一下赝势试试。。。
赞一下(7人) 回复此楼
科研需要激情
6楼2009-10-29 21:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):谢谢回帖交流 12-22 17:27
5.5和5.75差别很大吗？(5.75-5.5)/5.5=4.5%
赞一下(11人) 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
7楼2009-10-29 22:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

林子6682

新虫 (初入文坛)
各位大哥，我的U和J的值是从文献中找到的。这篇文献中算的J和U的值是用wien2k算的。我计算的时候选的是GGA＋U，和这两篇文献的一致。我采用这两篇文献中的U和J的值的时候，能够得到5.5的磁矩，但少数自旋的带隙（该体系为半金属）只有0.34eV，比别人的带隙（约1ev）要小很多，各位有很什么好建议么？
赞一下 回复此楼
8楼2009-12-17 22:49:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wo8165827

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by 林子6682 at 2009-12-17 22:49:
各位大哥，我的U和J的值是从文献中找到的。这篇文献中算的J和U的值是用wien2k算的。我计算的时候选的是GGA＋U，和这两篇文献的一致。我采用这两篇文献中的U和J的值的时候，能够得到5.5的磁矩，但少数自旋的带隙（ ...

我计算出来的带隙也是没有文献报道的宽，你的带隙宽度是怎么计算出来的。谢谢
赞一下(1人) 回复此楼
计算的魅力
9楼2009-12-22 16:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 林子6682 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 9/450 2026-08-02 16:32 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:32 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:31 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 16:27 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:21 by jHcrHMhbKCFV
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[基金申请] 面上提前没消息，有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
关闭
请填处理意见
关闭