| 查看: 1698 | 回复: 4 | |||
[交流]
****求助！***！！！关于金属表面的功函数计算和数据处理的问题 已有2人参与
|
|
我现在想做一组功函数对一组金属的得失电子能力进行比较 但是在进行功函数计算的时候遇到问题如下 1. 在vasp5.2里。计算的真空能级到底是用LVTOT = ture 还是LVHAR = ture 也就是说真空能级的计算是否应该包含exchange的部分 2. 我用这两个关键词算出来的，用 vtotav 处理后得到的都是 一个E-3数量级的非常小的真空能级的数值，很明显是不对的。我不知道是不是计算中出错 我的INCAR是这样的 1 Startparameter for this Run: 2 NPAR = 4 3 LPLANE = .TRUE. 4 ENCUT = 380 eV 5 6 NWRITE = 2; LPETIM=F write-flag & timer 7 ISTART = 0 job : 0-new 1-cont 2-samecut 8 ICHARG = 2 charge: 1-file, 2-atom, 10-const 9 ISPIN = 2 10 11 Electronic Relaxation 1 12 NELM = 60; 13 EDIFFG = -0.03 stopping-criterion for ELM 14 EDIFF = 1E-04 15 16 Ionic Relaxation 17 NSW = 300 number of steps for IOM 18 IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG 19 20 DOS related values: 21 ISMEAR = 0 22 SIGMA = .05 broadening in eV -4-tet -1-fermi 0-gaus 23 ALGO = Normal 24 25 # IDIPOL = 3 26 # LDIPOL = .TRUE. 27 # MAGMOM = 72*0.0 -6.0 -6.0 -6.0 6.0 6.0 6.0 6.0 28 29 LWAVE = F write WAVECAR 30 # LCHARG = T write CHGCAR 31 LCHARG = .FALSE. 32 LREAL=.FALSE. 33 LVTOT=.TRUE. 我的POSCAR是如下 Zn 2 1.00000000000000 3 2.2477611784190601 -1.2977455214342299 0.0000000000000000 4 -0.0000000000000020 2.5954910428684599 0.0000000000000000 5 0.0000000000000000 0.0000000000000000 36.5084650340021994 6 Zn 7 5 8 Direct 9 0.6666666666666643 0.3333333333333357 0.3515395121485313 10 0.6666666666666643 0.3333333333333357 0.6484604878514687 11 0.3333333333333357 0.6666666666666643 0.4253677314687111 12 0.3333333333333357 0.6666666666666643 0.5746322685312889 13 0.6666666666666643 0.3333333333333357 0.5000000000000000 我不知道错误在哪里 3 请问功函数的图怎么画出来的，我用vtotav 处理后只给了一个数值，请问如何保存数据或者画出功函数的那个变化图Sample Text [ Last edited by airyleay on 2013-4-2 at 16:46 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
27博士
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 纳米管功函数计算
已经有27人回复
- 求助 Origin 自定义函数输入问题
已经有4人回复
- 求助关于协作通信中节点数据功率的问题
已经有4人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 关于Origin 中输入新函数并拟合数据的问题，重谢200个
已经有8人回复
- 【求助】关于函数返回值的问题!
已经有3人回复
- 【求助】求助，如何计算液体100个原子的双体分布函数
已经有4人回复
- 【求助】Ce金属表面氧化膜的处理(晶体生长)
已经有3人回复
- 【求助】金属锌的费米能级是多少？
已经有11人回复
- 【求助】半导体功函数
已经有20人回复
2楼2013-04-02 16:45:36
3楼2013-04-02 17:04:02
fdd096030079
新虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 481.4
- 散金: 20
- 红花: 5
- 帖子: 112
- 在线: 63小时
- 虫号: 1973034
- 注册: 2012-09-04
- 专业: 凝聚态物性 II ：电子结构
4楼2016-12-03 10:42:20
5楼2018-12-25 10:55:11
