用VASP进行Partial Charge分析实例

vasp介绍(转帖)V ASP (V AMP)主页：http://cms.mpi.univie.ac.at/vasp/简介：V ASP是使用赝势和平面波基组，进行从头量子力学分子动力学计算的软件包，它基于CASTEP 1989版开发。

V AMP/V ASP中的方法基于有限温度下的局域密度近似（用自由能作为变量）以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。

这些技术可以避免原始的Car-Parrinello方法存在的一切问题，而后者是基于电子、离子运动方程同时积分的方法。

离子和电子的相互作用超缓V anderbilt赝势(US-PP)或投影扩充波(PA W)方法描述。

两种技术都可以相当程度地减少过渡金属或第一行元素的每个原子所必需的平面波数量。

力与张量可以用V AMP/V ASP很容易地计算，用于把原子衰减到其瞬时基态中。

功能：FeaturesSelf-consistent density functional method with plane wave basisAll-electron projector-augmented-wave (PA W) potentials covering all atoms of the periodic table Local density approximation (LDA) and generalized gradient approximation (GGA)Spin restricted and spin-polarizedSemi-relativistic and full spin-orbit relativisticNon-collinear magnetismLDA (GGA)+U for correlated systemsApplicable to bulk systems, surfaces, interfaces, and molecules (in supercell geometry)Total energies, forces, and the full stress tensorConcurrent relaxation of lattice parameters and atomic positionsAb initio molecular dynamicsGeneration of Monkhorst-Pack special k-pointsIntegration in k-space with smearing or tetrahedron method (with Blöchl correction)Band structure (spin restricted and spin polarized)Site, spin and partial-wave projected densities of statesPartial densities of states at Fermi levelCharge densities Spin densitiesUltra-soft pseudopotentials are offered for comparison平台：UNIX/LINUX相关软件：1. p4vasp主页：http://cms.mpi.univie.ac.at/odubay/p4vasp_site/news.php说明：处理V ASP xml格式输出文件的免费图形环境工具。

VASP计算实例目录一、氢气分子H2键长的计算 (3)1.基本文件 (3)2.赝势类型的选择 (3)3.截断能ENCUT参数的选择 (4)4.KPOINTS参数选择 (5)5.对晶格常数进行优化 (6)二、Si晶体晶格常数计算 (8)1.赝势类型选择 (8)2.截断能（ENCUT）参数的选定 (9)3.KPOINTS参数的选定 (11)4.SIGMA参数的选定 (12)5.晶格常数计算结果 (13)三、Si元素单原子能量计算 (14)1.由内聚能倒推单原子能量 (14)2.基本文件 (15)3.单原子能量计算 (15)四、Si的VASP力学常数计算 (16)1.计算所需文件 (16)2.计算与数据处理 (17)3.计算所用到的公式： (18)五、SI晶体的电子结构 (19)1.采用VASP计算能带的步骤 (19)2.电荷分布计算结果 (20)能带计算和结果 (21)3.态密度计算和结果 (21)六、Si晶体介电函数和光学性质的计算 (22)1.计算步骤 (22)2.用到的文件 (23)3.计算结果 (26)七、VASP的声子谱计算 (29)1.计算步骤 (29)2.基本文件 (30)3.声子谱、声子态密度计算和结果 (33)4.热学性质计算和结果 (34)八、化合物co2键长计算 (35)1.计算步骤 (35)2.基本文件 (35)一、氢气分子H2键长的计算1.基本文件准备基本文件INCAR、POTCAR、POSCAR、KPOINT以及脚本文件encut、k、optimize2.赝势类型的选择输入文件如下其中参数要靠经验初选INCAR:System = F2ISTART = 0ICHARG = 2NELMDL = 5ISMEAR = 0SIGMA = 0.1PREC = AccurateKPOINTS：Automatic meshM1 1 10 0 0POSCAR：O115.0 0.00 0.000.00 14.0 0.000.00 0.00 13.01D0.00 0.00 0.00分别选用五个贋势文件进行计算。

用VASP进行Partial Charge分析实例VASP Version : 4.6在这篇文章中，我将首先介绍Partial Charge的概念，以及如何用VASP具体的计算Partial Charge。

首先，所谓的Partial Charge是针对与Total Charge 来说的，指的是某个能量范围、某个K点或者某个特定的态所对应的电荷密度。

在文献中最常见的是价带顶部，导带底部，表面态或者局域态所对应的Partial Charge。

通过分析这些态所对应的Partial Charge，可以得到体系的一些性质，比如局域态具体的是局域在哪个原子上等。

我将通过具体的例子说明如何用VASP进行Partial Charge Analysis。

进行Partial Charge Analysis的第一步是进行自洽的计算，得到体系的电子结构。

这一步的计算采用通常的INCAR和KPOINTS文件。

在自洽计算结束后，我们需要保存WAVECAR文件。

（通过在INCAR文件中设置LWAVE=TRUE实现）在这个例子中，假设我们需要计算一个硅纳米线的导带和价带的Partial Charge。

硅纳米线的结构如下：第二步是画出能带结构，以决定你需要画哪条能带的那个K点的态所对应的Partial Charge。

关于具体如何用VASP画能带，请参见用VASP4.6计算晶体硅能带实例一文。

我们得到硅纳米线的能带结构如下：画能带时有些小技巧。

你可以用一些支持列模块的编辑器，如UltraEdit，将OUTCAR里的各个K点所对应的本征值粘贴到Origin中。

这一步完成后，在Origin 中做一个矩阵转置，然后将K点坐标贴到第一列，并将其设为X坐标。

如此画出来的基本上就是能带图了。

在Origin中可以通过设置纵轴范围来更加清楚的区分费米能级附近的各条能带。

如上的硅纳米线所对应的能带结构图如下：决定画哪条能带，或者那些感兴趣的K点之后，有如下几种方法计算不同的Partial Charge。

CHGCAR £leThis £le contains the lattice vectors, atomic coordinates, the total charge density multiplied by the volume (r) ¤V cell on the £ne FFT-grid (NG(X,Y,Z)F), and the PAW one-center occupancies. CHGCAR can be used to restart VASP from an existing charge density, for visualisation the CHG £le should be used, since the PAW-one centre occupancies are dif£cult to parse. It is possible to avoid that the CHGCAR is written by settingLCHARG = .FALSE.in the INCAR £le (see section 7.45). In VASP, the density is written using the following commands in Fortran:WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)The x index is the fastest index, and the z index the slowest index. The £le can be read format-free, because at least in new versions, it is guaranteed that spaces separate each number. Please do not forget to divide by the volume before visualizing the £le!For spinpolarized calculations, two sets of data can be found in the CHGCAR £le. The £rst set contains the total charge density (spin up plus spin down), the second one the magnetization density (spin up minus spin down). For non collinear calculations the CHGCAR £le contains the total charge density and the magnetisation density in the x, y and z direction in this order.For dynamic simulation (IBRION=0), the charge density on the £le is the predicted charge density for the next step: i.e. it is compatible with CONTCAR, but incompatible with the last positions in the OUTCAR £le. This allows the CHGCAR and the CONTCAR £le to be used consistently for a molecular dynamics continuation job. For static calculations and relaxations (IBRION=-1,1,2) the written charge density is the selfconsistent charge density for the last step and might be used e.g. for accurate band-structure calculations (see section 10.3).Mind: Since the charge density written to the £le CHGCAR is not the selfconsistent chargedensity for the positions on the CONTCAR £le, do not perform a bandstructure calculation (ICHARG=11) directly after a dynamic simulation (IBRION=0) (see section 10.3).Band decomposed charge density (parameters) 分波电荷密度计算VASP can calculate the partial (band decomposed) charge density部分（分波）电荷密度 according to parameters specified in the file INCAR. It must be noted, that the densities calculated by VASP (including the band decomposed charge density) are always symmetrized对称的 using the space group空间群 and point group symmetry determined by VASP.根据设置INCAR中的参数，VASP可以计算分波电荷密度，算分波电荷密度通常考虑空间群和点群对称性。

一、含义Partial charge density：也称为Band decomposed charge densityPartial charge density计算，就是计算特定的某个(或某些)k点和本征值(这些k点和本征值是相互对应的)所对应的本征波函数的平方(也就是电荷密度)。

二、计算过程和参数设置计算Partial电荷密度的方法：1.结构驰豫2.计算静态，输出CHGCAR、W A VECAR3.计算能带，无须输出CHGCAR、W A VECAR4.做能带图，选点。

5.计算Partial电荷密度，在INCAR中设置需要计算的能带编号和K点编号，不需要KPOINTS文件（也可以放一个算静态或者算能带的KPOINTS文件，都不影响结果），读入静态的CHGCAR和WA VECAR。

6.计算结束后，用PARCHG.****.****文件删掉后面数字，改成PARCHG作图。

参数设置：以上第一步和以往的计算一样，参数就不介绍了。

第4步：选点1.需要一个小程序（band_plot来处理能带计算输出的OUTCAR和EIGENV AL），或者一个可执行脚本（band_plot放到计算目录的bin/目录下），两者可以相互转化。

此处，用后面这种，将band_plot放到./bin目录下：2.然后到能带计算完成的目录下，输入命令：band_plot3.然后会得到一个文件数据文件：band_plot.dat4.将这个文件下载，用origin作图5.调整图片坐标6.我们根据需要在图中选择一点，然后找出其坐标，为了操作方便需要放大图片View—Zoom In:为了选点方便，采用点模式用左侧快捷键来得到需要点的坐标，如图【X=04714，Y=0.41041】然后，我们到数据表格中，搜索0.41041这个数，（因为X=04714是第一列，一目了然）Edit-Find，出现对话框，输入需要找的数：从表中，我们得知，第一列，对应X轴，也对应K点，我们算了31个K点，我们找的那一点，正好是第31个K点。