密度泛函理论

电子结构理论的经典方法，特别是Hartree-Fock方法和后Hartree-Fock方法，是基于复杂的多电子波函数的。密度泛函理论的主要目标就是用电子密度取代波函数做为研究的基本量。因为多电子波函数有 3N 个变量（N泛 为电子数，每个电子包含三个空间变量），而电子密度仅是三个变量的函数，无论在概念上还是实际上都更方便处理。

虽然密度泛函理论的概念起源于Thomas-Fermi模型，但直到Hohenberg-Kohn定理提出之后才有了坚实的理论依据。Hohenberg-Kohn第一定理指出体系的基态能量仅仅是电子密度的泛函。

最初的HK理论只适用于没有磁场存在的基态，虽然现在已经被推广了。最初的Hohenberg-K

ohn定理仅仅指出了一一对应关系的存在，但是没有提供任何这种精确的对应关系。正是在这些精确的对应关系中存在着近似（这个理论可以被推广到时间相关领域，从而用来计算激发态的性质[6]）。

密度泛函理论最普遍的应用是通过Kohn-Sham方法实现的。 在Kohn-Sham DFT的框架中，最难处理的多体问题（由于处在一个外部静电势中的电子相互作用而产生的）被简化成了一个没有相互作用的电子在有效势场中运动的问题。这个有效势场包括了外部势场以及电子间库仑相互作用的影响，例如，交换和相关作用。处理交换相关作用是KS DFT中的难点。目前并没有精确求解交换相关能 EXC 的方法。最简单的近似求解方法为局域密度近似(LDA)。LDA近似使用均匀电子气来计算体系的交换能（均匀电子气的交换能是可以精确求解的），而相关能部分则采用对自由电子气进行拟合的方法来处理。

自1970年以来，密度泛函理论在固体物理学的计算中得到广泛的应用。在多数情况下，与其他解决量子力学多体问题的方法相比，采用局域密度近似的密度泛函理论给出了非常令人满意的结果，同时固态计算相比实验的费用要少。尽管如此，人们普遍认为量子化学计算不能给出足够精确的结果，直到二十世纪九十年代，理论中所采用的近似被重新提炼成更好的交

换相关作用模型。密度泛函理论是目前多种领域中电子结构计算的领先方法。 尽管密度泛函理论得到了改进，但是用它来恰当的描述分子间相互作用，特别是范德瓦尔斯力，或者计算半导体的能隙还是有一定困难的。

[编辑] 早期模型: Thomas-Fermi 模型

密度泛函理论可以上溯到由Thomas和Fermi 在1920年代发展的Thomas-Fermi模型。他们将一个原子的动能表示成电子密度的泛函，并加上原子核-电子和电子-电子相互作用（两种作用都可以通过电子密度来表达）的经典表达来计算原子的能量。

Thomas-Fermi模型是很重要的第一步，但是由于没有考虑Hartree-Fock理论指出的原子交换能，Thomas-Fermi方程的精度受到限制。1928年保罗·狄拉克在该模型基础上增加了一个交换能泛函项。

然而，在大多数应用中Thomas-Fermi-Dirac理论表现得非常不够准确。其中最大的误差来自动能的表示，然后是交换能中的误差，以及对电子相关作用的完全忽略。

[编辑] 导出过程和表达式

在通常的多体问题电子结构的计算中，原子核可以看作静止不动的（波恩-奥本海默近似），这样电子可看作在原子核产生的静电势 中运动。电子的定态可由满足多体薛定谔方程的波函数 描述：

其中 为电子数目， 为电子间的相互作用势。算符 和 称为普适算符，它们在所有系统中都相同，而算符则依赖于系统，为非普适的。可以看出，单粒子问题和比较复杂的多粒子问题的区别在于交换作用项 。目前有很多成熟的方法来解多体薛定谔方程，例如：物理学里使用的图形微扰理论和量子化学里使用的基于斯莱特行列式中波函数系统展开的组态相互作用（CI）方法。然而，这些方法的问题在于较大的计算量，很难用于大规模复杂系统的计算。

对能量泛函 求变分极值可以用不定算子的拉格朗日方法，这由科恩和沈吕九在1965年完成 [2]。这里我们使用如下结论：上面方程中的泛函可以写成一个无相互作用的体系的密度泛函

可以得到一系列的电子轨域 ，并由此求得原来的多体体系的电子密度