| 查看: 1659 | 回复: 13 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yuyaolun木虫 (小有名气)
|
[交流]
【求助】怎么我算的Pt(100)的DOS尖峰这么多？
|
||
|
我想算Pt(100)界面的d-band center，应该是通过pdos来计算吧。 结构是3*3的supercell，四层一共36个原子。 我是这么做的，一，结构优化 INCAR PREC = Accurate ENCUT = 400 LREAL = AUTO ISMEAR = 1 NELMIN = 8 EDIFF = 1E-4 EDIFFG = -0.001 NSW = 100 IBRION = 1 ALGO = Fast NSIM = 4 LORBIT = 10 LWAVE = .FALSE. LCHARG = .FALSE. 这里有第一个问题： LORBIT具体应该怎么设定？在做结构优化的时候，我不需要DOSCAR和PROCAR，但是manual P66里LORBIT默认值是0，它又需要我输入RWIGS值，我不清楚应该怎么设定。当然这里LORBIT的设定并不影响结构优化的结果，可是下面的计算就出问题了。 二、用第一步的CONTCAR，静态自洽得到CHGCAR INCAR: PREC = Accurate ENCUT = 400 LREAL = AUTO ISMEAR = 1 NELMIN = 5 EDIFF = 1E-4 EDIFFG = -0.01 NSW = 0 IBRION = -1 ISTART = 0 ALGO = Fast NSIM = 4 LORBIT = 11 LWAVE = .FALSE. LCHARG = .TURE. KPOINTS: 6*6*6 这个有第二个问题： 按照manual P126的DOS calculation里做静态计算主要应该设定NSW=0和IBRION=-1，这里应该没有用到第一步计算的WAVECAR吧……可是P43 DOSCAR里为什么出现了这么一句话：Mind: For relaxations, the DOSCAR is usually useless. If you want to get an accurate DOS for the final configuration, first copy CONTCAR to POSCAR and continue with one static (ISTART=1; NSW=0) calculation. ISTART=1？那是说我应该算第一步的时候计算WAVECAR吗？是不是manual写错了？ 三、从第二步的CHGCAR计算DOS INCAR: PREC = Accurate ENCUT = 400 LREAL = AUTO ISMEAR = -5 NELMIN = 5 EDIFF = 1E-4 EDIFFG = -0.01 NSW = 100 IBRION = 1 ALGO = Fast NSIM = 4 ICHARG = 11 LORBIT = 11 LWAVE = .FALSE. LCHARG = .FALSE. 这是第三个问题： 我是应该设定LORBIT=11，不读RWIGS，还是应该设定LORBIT=0再给出RWIGS值？我用的PAW_PBE POTCAR。Manual在P126里写着Mind that VASP can calculate partial DOS. Partial DOS are very powerful for the analysis of the electronic DOS (see section 6.32---RWIGS). 可是在P65的RWIGS和LORBIT里写的很含糊。前面说算PDOS貌似应该给出RWIGS值，但是说LORBIT的时候又说对于PAW方法用不着给RWIGS？ 另外，Manual里写着如果我给了RWIGS值，那么应该是VASP writes a line Volume of Typ 1: 98.5 % to the OUTCAR file. You should use a RWIGS value which yields a volume of approximately 100%. 问题是我查了OUTCAR，怎么没有这句话呢？我怎么确定给的RWIGS是不是正确的？ 从POTCAR里可是有两个RWIGS值啊！ RWIGS = 2.750; RWIGS = 1.455 wigner-seitz radius (au A) 我应该用哪一个？ 还有，我需不需要设定NEDOS值？这个值有什么用？我看到版里有不少人用5000或者10000，而默认值只有300。 最后是我算的Pt DOS图，有很多尖峰，我不知道是哪一步算的不正确，希望牛人们耐着性子帮我讲解一下。。。谢谢！！ |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
yuyaolun
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2422.7
- 散金: 26
- 帖子: 57
- 在线: 13.2小时
- 虫号: 566220
- 注册: 2008-05-30
- 性别: GG
- 专业: 电化学
6楼2010-01-06 16:29:03
yuyaolun
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2422.7
- 散金: 26
- 帖子: 57
- 在线: 13.2小时
- 虫号: 566220
- 注册: 2008-05-30
- 性别: GG
- 专业: 电化学
2楼2010-01-06 16:23:06
yuyaolun
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2422.7
- 散金: 26
- 帖子: 57
- 在线: 13.2小时
- 虫号: 566220
- 注册: 2008-05-30
- 性别: GG
- 专业: 电化学
3楼2010-01-06 16:24:14
yuyaolun
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2422.7
- 散金: 26
- 帖子: 57
- 在线: 13.2小时
- 虫号: 566220
- 注册: 2008-05-30
- 性别: GG
- 专业: 电化学
4楼2010-01-06 16:26:36
20