如何用VASP计算单个原子的能量

in VASP all energies are referred to the the reference state for whichthe potential was generated!to determined the energy of the grounstate of the atom place the atomin a larger non cubic box to break initial symmetry (ie11Ax10Ax9A)use the G point onlyINCAR:ISPIN = 2 ! spin polarizedISMEAR = 0 ; SIGMA = 02 ! for small sigma conv for TM is diffMAGMOM = 2 ! initial magnetic momentone should use the energy value energy without entropy of the OUTCAR fileEcoh=E_total-nEatom大约翻译如下:在VASP中，所有的能量是相对于产生赝势的组态来说的，所以并不能将计算得到的能量当成真正的原子的基态能量。为了得到原子的基态能量，将原子放在一个非正方体的晶格内，比如说11Ax10Ax9A的晶格（这样做是为了消除简并）。KPOINTS只用G撒点大概的INCAR文件可以如下：ISPIN = 2 ! spin polarizedISMEAR = 0 ; SIGMA = 02 ! for small sigma conv for TM is diffMAGMOM = 2 ! initial magnetic moment最后计算结合能的公式应该为

使用的软件：VASP, VESTA, SshClient

Charge density difference of system AB: ∆ρ = ρAB − ρA − ρB

结构优化后 的自洽计算

A, B 各自保留优化后的 POSCAR 部分，编辑删除 POSCAR 文件

POTCAR 与 POSCAR 信息相符即可，KPOINTS 与自洽计算时相同，提交作业，计算完成后，获得三个 CHGCAR 文件，使用 SshClient 导出。

使用的软件：VASP, Origin, SshClient, vaspkit051

在完成结构优化和静态计算后，拷贝 scf 文件夹为 optic

编辑 optic 文件夹下 INCAR

提交作业，计算完成后，可查看 OUTCAR 文件，有计算得出介电常数 (Dielectric constant) 的实部和虚部。

使用 vaspkit 程序对计算结果进行处理。

输入命令 opt 产生 REALIN 和 IMAGIN 文件，依次输入命令 kit 和 51 ，屏幕显示如下

输出文件 ABSORBdat，REFRACTIVEdat，REFLECTIVITYdat，EXTINCTIONdat 和 ENERGYLOSSSPECTRUMdat，依次为 absorption coefficient, refractive coefficient, reflectivity coefficient, extinction coefficient and energy-loss function。导出使用 Origin 作图即可。

附： 相关光学性质计算公式

vasp中的与pwscf有区别吗？feierxide(站内联系TA)我的意思是d原来简并的5个轨道劈裂后的每个轨道对应的态密度怎么画出……很多文献都给了这样的结果coogert(站内联系TA)== 不清楚。 帮顶下。。gdflly08(站内联系TA):tiger05:bingmou(站内联系TA)Originally posted by feierxide at 2011-03-30 16:55:16:我的意思是d原来简并的5个轨道劈裂后的每个轨道对应的态密度怎么画出……很多文献都给了这样的结果 在VASP中就是Pdos这个方法来获得的，具体方法是通过局部的原子轨道和体系的波函数进行积分得到，所以选坐标很重要，因为原子轨道依赖于坐标系jugengfans(站内联系TA)按照xy yz z2 xz x2顺序 排列 如果有自旋的话每一列轨道变成紧挨着的两列。如果分析劈裂，请将晶体场的x,y,z轴与你坐标系的x,y,z轴对齐。具体方法是在vasp的poscar文件里面旋转基矢huazhorg(站内联系TA)每个轨道都有一列单独的数据，在gnuplot中，用plot 'filename' u 1:2 w l, 'filename' u 1:3 w l,'filename' u 1:4 w l，……这样就可以把所有的轨道用不同的颜色画在一张图里了。fengshiquan(站内联系TA)你的意思是不是想把你图上的五个d轨道的劈裂态密度分开做，这可以在origin中做到goldenfisher(站内联系TA)用projwfcx这个程序来算，pdos就出来了franch(站内联系TA)问题是这些轨道有重叠的情况下，，怎么区分画出来？feierxide(站内联系TA)输出文件中，每个轨道都有一列单独的数据，分别画就行了:tiger07:franch(站内联系TA)Originally posted by feierxide at 2011-03-31 19:49:12:ms的确没有，用pwscf，好像wien2k也可以:) 看了一下例8，计算之后有个文件：nipdos_atm#1(Ni)_wfc#1(s)看了一下例8，计算之后有个文件：nipdos_atm#1(Ni)_wfc#1(s)每个轨道都有一列单独的数据，在gnuplot中，用plot 'filename' u 1:2 w l, 'filename' u 1:3 w l,'filename' u 1:4 w l，……这样就可以把所有的轨道用不同的颜色画在一张图里了。

VASP是采用的Berry phase方法来计算半导体和绝缘体材料的Born有效电荷。这种方法的介绍可以参考VASP手册上提到的文献。
如VASP手册上的介绍，在采用Berry phase方法来计算某个原子某个方向的Born有效电荷时，有两大步：
1) 计算离子未移动时(即平衡状态时)的极化(polarization); 2)手动移动所要计算的原子的位置，计算由原子偏离平衡态时的极化。每一大步中有包括了四个小步骤: i)自恰计算体系以得到收敛的电荷密度(CHG, CHGCAR文件)和波函数(WAVECAR); ii)计算G1方向上的电子极化; iii)计算G2方向上的电子极化; iv)计算G3方向上的电子极化。
最后从每一小步中的OUTCAR文件中找出极化各个部分的贡献: ev (电子项), bp (Berry phase项), ion (离子项)，进行简单的数值加减计算，按差分方法计算得到Born有效电荷。 对每类原子的各个方向按同样的步骤进行计算处理，在计算时可根据晶体的对称性，只需要计算Born有效电荷张量的某些分量即可。

---