DF计算简例分析.ppt

DFT计算简例分析 1D氢链 1D氢链的DFT计算 1D氢链的DFT计算 Mokhorst-Pack mesh点：为简单起见只取1个k-point。 初始电荷密度：假设n(r)是均匀的， n(r)=Z/a。 确定基函数，构造H, S并解Kohn-Sham方程。 1D氢链的DFT计算 基函数 取 1D氢链的DFT计算 Hamiltonian: 1D氢链的DFT计算 计算 其中动能 1D氢链的DFT计算 为了使Hamiltonian更加简单，下面假定（物理上未必合适） 假定 中所有的项都由赝势的项 所支配。 取 写 假定Vpp(G)为实，使得Vpp(-G)= Vpp(G) 1D氢链的DFT计算 将2x2矩阵H（k=0）对角化，得到本征矢和本征值如下 1D氢链的DFT计算 计算新密度 因为每个原子只有1个电子，所以最低能带是占有带，电子密度是 * a V(r) a=格常数 L=体系的总长度 V(r)=周期性（吸引）赝势 G=倒格矢 Z=1,价电子数 a V(r) a=格常数 L=体系的总长度 V(r)=周期性（吸引）赝势 G=倒格矢 Z=1,价电子数 可见，T是一个对角的2x2矩阵。 因为Vpp(2?/a)0，所以最低能态是 这个n(out)与所用的n(in)是相当不同的，将n(in) 与 n(out)混合可以得到一个新的n(r)。计算新的Veff，并开始下一次迭代。一直到自洽。 *