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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于稀土元素的4f轨道
161/2返回列表12下一页
查看: 3977  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

含草

金虫 (小有名气)

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

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

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

第一原理资料汇编

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-09 23:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencemaji

至尊木虫 (文坛精英)
★
含草(金币+2):谢谢参与
含草(金币+1): 2010-09-10 20:24:41
你能把你的参数设置贴上来吗？
赞一下(1人) 回复此楼
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否？
赞一下(1人) 回复此楼
阿什顿联
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电子
赞一下(2人) 回复此楼
4楼2010-09-10 08:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellsaking

金虫 (正式写手)
★
含草(金币+2):谢谢参与
可能是POTCAR有问题吧？
赞一下(1人) 回复此楼
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作为核处理了
赞一下(2人) 回复此楼
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电子，你的参数再多也没有用啊！
赞一下(2人) 回复此楼
大家一起学才好。
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的回帖
相关版块跳转 我要订阅楼主 含草 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +13 koalala 2026-08-24 15/750 2026-08-24 13:04 by lugia2
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭