人们从化学实践出发，对分子的电子结构一直非常关注。多年来，人们总是想探讨分子中原子的性质，特别是分子中原子的净电荷，提出了各种计算方法，并进行了不断地改进，然而至今都还不能令人满意。特别是原子净电荷划分、计算还远不够完善合理，使得净电荷有时与实验的预期相差太远，许多净电荷异常到不可思议的地步，致使人们对量化计算出的电荷失去信心。

在量化图形中，电子密度差Δρ尤其重要！它应是传统的原子净电荷概念的扩大和空间化，它直接从MO而来，真实地、科学抽象地表现了电荷在分子空间分布的净变化，不再只将净电荷人为地局限、归结、划分在某个原子上。

我们如果仅看分子轨道图、电子密度图，象上面的HF分子的σ键电子密度图，它是据HF分子第3个MO即σ键作的电子密度图，它们那么浑圆、憨厚、质朴，很难看出什么。它是一片云，难察其意。而一旦作出了电子密度差，成键原子前后变化、分子之间络合、氢健、吸附等，电子增减的净变化，被鲜明地、细微地突出了，可以发现许多新东西、意想不到的东西，从而引人入思。

如从上面的HF分子的σ键电子密度差图Δρ，并参照MO系数来看,这个σ键可以分析为三个内含：

第一、原子以基组

14

9

VASP中电子态密度计算的流程

主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算

以Al－FCC为例子

第一步 结构优化

输入文件（INCAR, POTCAR, POSCAR, KPOINT）    

INCAR文件

System=Al

ISTART=0

ISMEAR=1

SIGMA=0.2

ISPIN=2

GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01

IBRION=2

NSW=50

ISIF=2 (OR 3)

NPAR=10

POTCAR 文件直接在势库中拷贝

POSCAR文件

Al

4.05

1.0  0.0  0.0

0.0  1.0  0.0

0.0  0.0  1.0

4

Direct

0.0        0.0  0.0

0.5        0.5  0.0

0.5        0.0  0.5

0.0  0.5  0.5

KPOINT 文件

Automatic generation