24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

airyleay

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6.5
帖子: 22
在线: 11.9小时
虫号: 1591021
注册: 2012-01-30
性别: MM
专业: 理论和计算化学

[交流] **求助！*！！！关于金属表面的功函数计算和数据处理的问题已有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 ]

回复此楼