VASP中计算DOS时的设置问题

VASP计算DOS时的参数设置

2012/3/25

LORBIT = 110.4s-PDOS0.30.20.10.00.4-10-50510150.30.20.10.02.8-102.42.01.61.20.80.40.0-10-5051015p-PDOSd-PDOSISMEAR = -5ISMEAR = -1ISMEAR = 2-5051015E-EF (eV)

图1. 相同LORBIT值不同ISMEAR值时Pt(111)表层原子的PDOS

ISMEAR = 20.3s-PDOS0.20.10.00.20-10-50510150.150.100.050.002.8-102.42.01.61.20.80.40.0-10-5051015p-PDOSd-PDOSLORBIT = 1 (RWIGS) LORBIT = 11-5051015E-EF (eV)

图2. 同一ISMEAR值不同LORBIT值时Pt(111)表层原子的PDOS

以上计算均基于k点为Gamma_6x6x1时优化的构型（p(2x2) AL4 VL6模型），随后再做一次静态计算，此时取k点为Gamma_12x12x1计算DOS。

由图1可以看出，当ISMEAR取值为-1时（Fermi smearing with SIGMA=0.1）时DOS分布最为光滑，取值为2时（Methfessel-Paxton scheme of order 2 with SIGMA = 0.2）其次，最精确的是取值为-5时（tetrahedron method with Bl?chl correction）。此外，MP方案中的能带宽度偏窄，且当k点去的较小时，会出现负的态密度（如图3所示）。

ISMEAR = 2; SIGMA = 0.2; LORBIT = 10.4s-PDOS0.30.20.10.00.4-10-5051015p-PDOS0.30.20.10.0-102.4-5051015d-PDOS1.81.20.60.0-10-505Gamma_12x12x1Gamma_6x6x11015E-EF (eV)

图3. 取不同k点时Pt(111)表层原子的PDOS