动能密度

本文整理自 Aurélien Lherbier 于 2009 年撰写的两份 Abinit 技术报告 \cite{lherbier2009implementation,lherbier2009test}：先回顾平面波电子结构计算的基础，以及波函数、电子密度和动能密度在晶体对称性下的处理方式；再在此基础上构造电子局域函数（ELF），并用孤立 H 原子和 Li 原子检验它在 Abinit 模守恒赝势框架中的实现。

$$ ELF(\mathbf{r}) = \frac{1}{1+\left(\dfrac{D(\mathbf{r})}{D^0(\mathbf{r})}\right)^2} $$