24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3516  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 含草 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 食品工程专硕一志愿中海洋309求调剂 +4 小张zxy张 2026-03-26 8/400 2026-03-28 19:25 by lbsjt
[考研] 0856求调剂 +11 zhn03 2026-03-25 12/600 2026-03-28 13:32 by 唐沐儿
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 311求调剂 +3 希望上岸阿小杨 2026-03-23 3/150 2026-03-28 07:57 by 热情沙漠
[考研] 085404求调剂,总分309,本科经历较为丰富 +4 来财aa 2026-03-25 4/200 2026-03-28 07:41 by 棒棒球手
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 考研调剂 +10 呼呼?~+123456 2026-03-24 10/500 2026-03-27 11:46 by wangjy2002
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
信息提示
请填处理意见