24小时热门版块排行榜    

查看: 3355  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

含草

金虫 (小有名气)

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

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

含草

金虫 (小有名气)

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