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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP计算体系性能没有f能带电子
251/1返回列表
查看: 6005  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +12 生如夏花… 2026-03-11 14/700 2026-03-17 10:56 by lovewei0727
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +5 大火山小火山 2026-03-16 7/350 2026-03-17 10:26 by peike
[考研] 化学工程321分求调剂 +10 大米饭！ 2026-03-15 13/650 2026-03-17 06:57 by laoshidan
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
关闭
请填处理意见
关闭