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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP计算体系性能没有f能带电子
251/1返回列表
查看: 6959  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
信息提示
关闭
请填处理意见
关闭