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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP计算体系性能没有f能带电子
251/1返回列表
查看: 6865  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

langren24

新虫 (小有名气)

[交流] 【求助】VASP计算体系性能没有f能带电子

求助：用VASP计算SmCo5磁学性质，始终没能出现f能带电子，这是为什么呢？
         求高手指教？

这是我的INCAR文件：文件代码特别是磁学设置有什么问题吗？求高手指教
还有就是MAGMOM的设置，磁矩的大小，哪些原子是自旋向上，哪些原子是自选向下啊？
# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
SYSTEM = Sm Co5  (P6/mmm) ~ 947 minimized: C:/Documents and Settings/Administrator/Lb/SmCo5-107255(2001).cif (VASP)

    PREC = Accurate
   ENCUT = 500.000
  IBRION = -1
     NSW = 0
    NELM = 60
  NELMIN = 2
   EDIFF = 1.0e-06
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
    ALGO = Normal (blocked Davidson)
LSORBIT = .TRUE.
   SAXIS = 1 0 0
MAGMOM = 2*0 2 2*0 2 2*0 2 2*0 2 2*0 2 2*0 -2
    LDAU = .TRUE.
LDAUTYPE = 1
   LDAUL = 2 3
   LDAUU = 4.2 7.1
   LDAUJ = 0.8 1.1
LMAXMIX = 6
  INIWAV = 1
  ISTART = 0
  ICHARG = 0
   LWAVE = .TRUE.
  LCHARG = .TRUE.
ADDGRID = .FALSE.
   LVTOT = .TRUE.
    LELF = .TRUE.
  ISMEAR = 1
   SIGMA = 0.2
   LREAL = .FALSE.
RWIGS = 1.16 1.62

[ Last edited by langren24 on 2011-4-7 at 10:20 ]
回复此楼

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

first principle

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-04-07 10:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六198

金虫 (正式写手)
★
zzy870720z(金币+1): 鼓励讨论交流 2011-04-08 08:33:56
langren24(金币+1): 2011-04-08 15:08:47
langren24(金币+1): 2011-04-08 15:35:12
也许是赝势的问题，呵呵，我猜的！
赞一下(1人) 回复此楼
2楼2011-04-07 18:12:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★
youzhizhe(金币+1): 谢谢提示。 2011-04-07 21:43:50
站内搜索，前面有这个问题的讨论。
赞一下(1人) 回复此楼
3楼2011-04-07 21:20:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
★
bingmou(金币+1): 感谢交流！ 2011-04-18 17:02:23
关于VASP计算体系性能只有s p d没有出现f能带电子 ，站内搜索并没有相关的帖子！
赞一下(1人) 回复此楼
4楼2011-04-08 15:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanglei200

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
langren24(金币+1): 2011-04-19 15:47:09
langren24(金币+1): 2011-04-19 15:50:49
感觉你的INCAR设置有问题啊!NSW=0,INRION=-1，是静态的计算的计算啊，怎么后面又写NELM，NELMIN呢？？
赞一下(1人) 回复此楼
5楼2011-04-18 16:54:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanglei200

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by langren24 at 2011-04-08 15:16:08:
关于VASP计算体系性能只有s p d没有出现f能带电子 ，站内搜索并没有相关的帖子！

sorry，我看错了。这点没问题。
赞一下(1人) 回复此楼
6楼2011-04-18 16:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
ellsaking(金币+1): 谢谢交流 2011-04-19 10:21:39
langren24(金币+1): 2011-04-19 15:47:13
INCAR上面貌似没有什么问题 可能是你选择的赝势的价电子数太低了
赞一下(2人) 回复此楼
7楼2011-04-19 09:31:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by liuhuimeiwoo at 2011-04-19 09:31:41:
INCAR上面貌似没有什么问题 可能是你选择的赝势的价电子数太低了

选择赝势的价电子数太低了应该怎么解决呢？是换一种赝势，还是修改赝势的设置啊？求教了
赞一下 回复此楼
8楼2011-04-19 15:49:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)
换一种赝势 除非是你自己写一个 要不然你直接改的话 是没用的
赞一下 回复此楼
9楼2011-04-20 09:29:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
zhang668(金币+1): 多谢提醒 2011-04-21 17:12:37
还有一种可能性 我是突然想起来的 我记得我以前做加U计算的时候  如果加在了d电子上 f电子上随便你怎么写都没用了 一个元素不能加两个轨道的U  但是我不知道你这个两个元素上加了不同轨道的u会不会也有影响 还是很久以前做过的了
赞一下(2人) 回复此楼
10楼2011-04-20 17:18:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by liuhuimeiwoo at 2011-04-20 17:18:49:
还有一种可能性 我是突然想起来的 我记得我以前做加U计算的时候  如果加在了d电子上 f电子上随便你怎么写都没用了 一个元素不能加两个轨道的U  但是我不知道你这个两个元素上加了不同轨道的u会不会也有影响 还是很 ...

是啊，我用LSDA+U,一个元素U加在d电子上，一个元素U加在f电子上，难道这会导致不能出现f电子吗？还是那个d电子元素不用加U？
赞一下 回复此楼
11楼2011-04-21 12:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
gzqdyouxia(金币+2): 谢谢交流 2011-04-22 08:03:08
还是很久以前做的了 我也不确定 我当时是在同一个原子上加的 不知道不同原子上加有没有影响 你可以这样测试一下：你把那个d上的去掉 看会不会出来f电子的结果 要是还是没有 就不是这个问题了 当然结果不对 只是用来测试加在不同轨道上会不会有影响的
赞一下(2人) 回复此楼
12楼2011-04-21 16:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

howberg

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
gzqdyouxia(金币+1): 谢谢提示 2011-04-22 08:03:18
貌似是赝势的问题，有的赝势里根本就不包含f电子，去哪里找f电子的性质，？
你找找有f电子的赝势，不要管元素种类，看看算出来有没有f电子这一项，我好像记得我玩过这个，记不清楚了！
赞一下(2人) 回复此楼
13楼2011-04-21 22:15:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by howberg at 2011-04-21 22:15:51:
貌似是赝势的问题，有的赝势里根本就不包含f电子，去哪里找f电子的性质，？
你找找有f电子的赝势，不要管元素种类，看看算出来有没有f电子这一项，我好像记得我玩过这个，记不清楚了！

VASP里的赝势不是自带的吗，比如GGA-PBE，有f电子的赝势是什么意思啊?

我现在使用的只能选择GG-AMO5 GGA-PBE GGA-PBEsol GGA-rPBE这几种

[ Last edited by langren24 on 2011-4-22 at 16:25 ]
赞一下 回复此楼
14楼2011-04-22 15:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

howberg

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by langren24 at 2011-04-22 15:10:09:
VASP里的赝势不是自带的吗，比如GGA-PBE，有f电子的赝势是什么意思啊?

我现在使用的只能选择GG-AMO5 GGA-PBE GGA-PBEsol GGA-rPBE这几种

[ Last edited by langren24 on 2011-4-22 at 16:25 ]

pot      potpaw_GGA    pot_GGA  potpaw_PBE    potpaw
这是我们所拥有的赝势，有没有f电子，你Vi  一下那个赝势不就看到了?
   VRHFIN =U:  [Xe,5d,4f]
这是potpaw_GGA下的U元素的赝势里看到的，说明他有f电子，是在Xe核的基础上加了d电子和f电子而生成的赝势;
   VRHFIN =La : [core=Kr4d]
这是potpaw_GGA下的La元素的赝势里看到的，说明他没有f电子，只是在Kr核的基础上加了d电子而已.
赞一下(1人) 回复此楼
15楼2011-04-22 18:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
若问下楼主，怎么确认有没有出现f电子？看能带图吗？
赞一下(1人) 回复此楼
16楼2011-04-22 19:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by howberg at 2011-04-22 18:56:45:
pot      potpaw_GGA    pot_GGA  potpaw_PBE    potpaw
这是我们所拥有的赝势，有没有f电子，你Vi  一下那个赝势不就看到了?
   VRHFIN =U:  [Xe,5d,4f]
这是potpaw_GGA下的U元素的赝势里看到的，说明他有 ...

那么具体怎么找有f电子的赝势啊？应该怎么操作啊？
赞一下 回复此楼
17楼2011-04-25 14:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by howberg at 2011-04-22 18:56:45:
pot      potpaw_GGA    pot_GGA  potpaw_PBE    potpaw
这是我们所拥有的赝势，有没有f电子，你Vi  一下那个赝势不就看到了?
   VRHFIN =U:  [Xe,5d,4f]
这是potpaw_GGA下的U元素的赝势里看到的，说明他有 ...

这是我的赝势文件POTCAR:
VRHFIN =Sm : [core=Xe4]是什么意思啊？是不是没有f电子啊？

POTCAR:
PAW_PBE Sm_3 07Sep2000
11.0000000000000000
parameters from PSCTR are:
   VRHFIN =Sm : [core=Xe4]
   LEXCH  = PE
   EATOM  =   973.0754 eV,   71.5190 Ry

   TITEL  = PAW_PBE Sm_3 07Sep2000
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        1    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =    2.400    partial core radius
   POMASS =  150.360; ZVAL   =   11.000    mass and valenz
   RCORE  =    2.820    outmost cutoff radius
   RWIGS  =    3.000; RWIGS  =    1.588    wigner-seitz radius (au A)
   ENMAX  =  182.035; ENMIN  =  136.526 eV
   RCLOC  =    2.413    cutoff for local pot
   LCOR   =        T    correct aug charges
   LPAW   =        T    paw PP
   EAUG   =  290.851
   DEXC   =    -.079
   RMAX   =    3.417    core radius for proj-oper
   RAUG   =    1.300    factor for augmentation sphere
   RDEP   =    2.877    radius for radial grids
   QCUT   =   -3.658; QGAM   =    7.316    optimization parameters
赞一下 回复此楼
18楼2011-04-25 22:41:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

howberg

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by langren24 at 2011-04-25 22:41:14:
这是我的赝势文件POTCAR:
VRHFIN =Sm : [core=Xe4]是什么意思啊？是不是没有f电子啊？

POTCAR:
PAW_PBE Sm_3 07Sep2000
11.0000000000000000
parameters from PSCTR are:
   VRHFIN =Sm : [core=Xe4 ...

他那个4是什么意思了，没说是s电子，也没说是f电子，但是对于Sm来讲Xe原子实再加上2个s电子和5个f电子才对啊，这个我没见过，不太懂，你在自己的机器上试试吧！
用有f电子的试试，要是计算出现了f电子，用同样的设置看看Sm有啥特点！
赞一下(1人) 回复此楼
19楼2011-04-26 10:15:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

howberg

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
对了，计算磁学性质用wien2k比较好，vasp计算磁学性质比不了wien2k,wien2k是全电子计算，磁学性质好的不得了！
赞一下(1人) 回复此楼
20楼2011-04-26 10:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saikun

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
对了，你这个potcar是没有F电子的。
价电子是
5s(2) 5p(6) 5d(1) 6s(2)

你应找个有F电子的potcar来用（但cutoff不敢想象! )
赞一下(1人) 回复此楼
21楼2011-04-26 15:36:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
请问怎么看s, p ,d ,f 电子的能级啊
赞一下(1人) 回复此楼
22楼2011-04-26 17:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by saikun at 2011-04-26 15:36:45:
对了，你这个potcar是没有F电子的。
价电子是
5s(2) 5p(6) 5d(1) 6s(2)

你应找个有F电子的potcar来用（但cutoff不敢想象! )

"你应找个有F电子的potcar来用"应该怎么找啊，还有我的potcar只有Sm_3  Sm两种啊？我用的是Sm_3
赞一下 回复此楼
23楼2011-04-26 19:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saikun

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by langren24 at 2011-04-26 19:57:10:
"你应找个有F电子的potcar来用"应该怎么找啊，还有我的potcar只有Sm_3  Sm两种啊？我用的是Sm_3

"Sm" 有F电子

5s(2) 6s(2) 5p(6) 5d(1) 4f(5)

16个价电子,但 EATOM=2218.8 eV!
赞一下(1人) 回复此楼
24楼2011-04-26 21:34:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langren24

新虫 (小有名气)
引用回帖:
Originally posted by saikun at 2011-04-26 21:34:04:
"Sm" 有F电子

5s(2) 6s(2) 5p(6) 5d(1) 4f(5)

16个价电子,但 EATOM=2218.8 eV!

你有Sm_3的赝势吗，应该如何解决呢？现在真的很急迫解决这个问题啊，否则后面的工作没法进行了！
赞一下 回复此楼
25楼2011-04-27 14:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 langren24 的主题更新
251/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 人气不行了 +9 fansofjerry 2026-08-21 9/450 2026-08-24 17:48 by zhaozhengxin
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭