24小时热门版块排行榜

返回列表

当前主题已经存档。

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

yuyaolun

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2422.7
散金: 26
帖子: 57
在线: 13.2小时
虫号: 566220
注册: 2008-05-30
性别: GG
专业: 电化学

[交流] 【求助】怎么我算的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图，有很多尖峰，我不知道是哪一步算的不正确，希望牛人们耐着性子帮我讲解一下。。。谢谢！！

回复此楼