24小时热门版块排行榜    

查看: 3302  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 含草 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见