24小时热门版块排行榜    

查看: 3897  |  回复: 15

含草

金虫 (小有名气)

[交流] 【求助】关于稀土元素的4f轨道 已有3人参与

大家好,我在计算中采用了PAW GGA-PBE势, 来计算含有稀土元素的磁性性质,计算结果中发现根本没有考虑到4f 电子,稀土元素的离子磁矩基本上为0,而且在PDOS中也没有4f轨道的贡献,请问这是什么原因?谢谢
回复此楼

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

第一原理资料汇编

» 猜你喜欢

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

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

sciencemaji

至尊木虫 (文坛精英)


含草(金币+2):谢谢参与
含草(金币+1): 2010-09-10 20:24:41
你能把你的参数设置贴上来吗?
2楼2010-09-10 07:23:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

163.com

金虫 (著名写手)

实习版主


含草(金币+2):谢谢参与
含草(金币+1): 2010-09-10 20:24:44
含草(金币+1): 2010-09-10 20:24:50
VASP是可以算4f电子的,你用的是VASP否?
阿什顿联
3楼2010-09-10 08:41:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)

★ ★
含草(金币+2):谢谢参与
zzy870720z(金币+1):谢谢提示 2010-09-10 09:17:11
含草(金币+1): 2010-09-10 20:24:56
检查下POTCAR中是否包含4f电子
4楼2010-09-10 08:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellsaking

金虫 (正式写手)


含草(金币+2):谢谢参与
可能是POTCAR有问题吧?
ihaveadream
5楼2010-09-10 08:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

★ ★
含草(金币+2):谢谢参与
含草(金币+1): 2010-09-10 20:25:00
163.com(金币+1):谢谢回帖交流 2010-09-10 22:09:27
有可能选的赝势中4f作为核处理了
6楼2010-09-10 09:15:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)

ISTART = 0
ICHARG = 2

ENCUT=500
PREC=Accurate    medium,high,low
ISMEAR=-5
EMIN=-20
EMAX=20
NEDOS=500

ISPIN=2
MAGMOM=8*5  1*9  8*-7  7*7  32*0

ALGO =  Fast   # RMM-DIIS algorithm for electrons
LREAL=Auto    real-space projection
EDIFF = 1E-4 ; EDIFFG = -1E-2

#RWIGS=1.302 1.302 0.820
NPAR=1
LORBIT=12

AMIX=0.2
BMIX=0.0001  !most zero,but 0 will crash some versions
AMIX_MAG=0.8
BMIX_MAG=0.0001 !most zero,but 0 will crash some versions

LDAU=.TRUE.
LDAUTYPE=2
LDAUL=2 3 2 -1
LDAUU=4.000 6.000  4.000  0.000
LDAUJ=0.000  0.000  0.000  0.000
LDAUPRINT=2

LMAXMIX=6
LWAVE=.T.     write WAVECAR

上面是我的设置,我采用的是VASP,我这里面含有稀土Dy元素,但是好象只有Dy-3赝势,算出来的根本就没有f轨道的贡献。那我应该怎么办,不知道有没有出来新的势,就是考虑了4f电子。
7楼2010-09-10 20:30:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fantouwen

金虫 (小有名气)

将军

★ ★
含草(金币+2):谢谢参与
cenwanglai(金币+1):en 2010-09-11 11:48:18
LZ要报把你的POTCAR弄来给大家看看,如果POTCAR本来就没有f电子,你的参数再多也没有用啊!
大家一起学才好。
8楼2010-09-11 11:00:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)

PAW_PBE Dy_3 06Sep2000
9.00000000000000000
parameters from PSCTR are:
   VRHFIN =Dy : [core=Xe4]
   LEXCH  = PE
   EATOM  =   651.0671 eV,   47.8521 Ry

   TITEL  = PAW_PBE Dy_3 06Sep2000
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        1    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =    2.300    partial core radius
   POMASS =  162.500; ZVAL   =    9.000    mass and valenz
   RCORE  =    3.000    outmost cutoff radius
   RWIGS  =    3.000; RWIGS  =    1.588    wigner-seitz radius (au A)
   ENMAX  =  155.729; ENMIN  =  116.797 eV
   RCLOC  =    2.304    cutoff for local pot
   LCOR   =        T    correct aug charges
   LPAW   =        T    paw PP
   EAUG   =  276.179
   DEXC   =    -.043
   RMAX   =    3.695    core radius for proj-oper
   RAUG   =    1.300    factor for augmentation sphere
   RDEP   =    3.148    radius for radial grids
   QCUT   =   -3.383; QGAM   =    6.766    optimization parameters

   Description
     l     E      TYP  RCUT    TYP  RCUT
     0   .000     23  3.000
     0   .800     23  3.000
     1   .000     23  3.000
     1   .000     23  3.000
     2   .000     23  3.000
     2   .000     23  3.000
   Error from kinetic energy argument (eV)
   NDATA  =      100
   STEP   =   20.000   1.050
  63.0      60.2      58.8      55.9      54.4      51.4      48.4      46.8
  43.8      40.7      39.2      36.3      33.4      30.6      27.9      25.3
  22.9      20.5      18.4      15.4      13.6      11.9      9.70      8.39
  7.22      5.69      4.42      3.70      2.80      2.08      1.52      1.09
  .762      .525      .355      .237      .138      .936E-01  .597E-01  .460E-01
  .363E-01  .311E-01  .274E-01  .238E-01  .200E-01  .161E-01  .125E-01  .873E-02
  .649E-02  .465E-02  .364E-02  .321E-02  .291E-02  .266E-02  .229E-02  .192E-02
  .153E-02  .111E-02  .822E-03  .654E-03  .575E-03  .533E-03  .489E-03  .414E-03
  .337E-03  .253E-03  .194E-03  .163E-03  .149E-03  .138E-03  .119E-03  .936E-04
  .716E-04  .585E-04  .526E-04  .494E-04  .429E-04  .341E-04  .260E-04  .226E-04
  .211E-04  .189E-04  .153E-04  .119E-04  .104E-04  .973E-05  .846E-05  .659E-05
  .549E-05  .515E-05  .453E-05  .364E-05  .303E-05  .285E-05  .245E-05  .195E-05
  .174E-05  .161E-05  .129E-05  .110E-05
END of PSCTR-controll parameters
local part
62.3609576231044684
应该是没有带f电子,但是对于Dy暂时只有这一个赝势,不知道有没有出新的赝势。
9楼2010-09-12 09:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chimegreen

木虫 (正式写手)

Dr.


含草(金币+2):谢谢参与
含草(金币+1): 2010-09-13 09:25:50
potcar的问题
快乐科研,快乐生活。
10楼2010-09-12 21:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 含草 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +73 AI8RaGaPaSCR 2026-08-07 76/3800 2026-08-14 08:45 by IUa8yhHGR7hP
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 30/1500 2026-08-14 07:55 by 医学老男孩
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见