1 历史与发展
1.1 早期量子理论背景
密度泛函理论的形成，建立在早期量子力学与多体问题研究的基础之上。20世纪上半叶，人们已经认识到电子之间的相互作用会显著影响原子、分子和固体的性质，但直接求解多电子薛定谔方程通常极为困难。为降低复杂度，研究者先后提出了Hartree方法、Hartree-Fock近似等方案，尝试用较少的变量描述电子行为。
这些早期方法虽然推动了计算量子化学的发展，但仍然需要处理多电子波函数，计算量会随着电子数迅速增长。由此，如何用更简洁的物理量来表征体系，成为后续理论突破的重要方向。电子密度作为一种只依赖空间坐标的量，因其直观性和可操作性逐渐受到重视。
1.2 Hohenberg-Kohn定理的提出
1964年，Hohenberg和Kohn提出了两条奠基性定理，证明在基态条件下，电子密度可以唯一决定体系的全部性质。这个结论从根本上改变了多体量子问题的研究思路：与其直接追踪高维波函数，不如将注意力集中在三维空间中的电子密度上。
这一结果具有重要意义，因为它为“以密度代替波函数”提供了严格理论基础。自此之后，研究者能够把复杂体系的基态能量、结构和其他物理量，表达为电子密度的泛函，从而打开了密度泛函理论的研究道路。
1.3 Kohn-Sham框架的建立
尽管Hohenberg-Kohn定理证明了电子密度的重要性，但并未直接给出可实际计算的形式。1965年，Kohn和Sham提出了一套更具可操作性的框架：把真实的相互作用多电子体系，映射为一个非相互作用电子体系，通过引入有效势来重现真实密度。
这一方案极大增强了理论的实用性。Kohn-Sham框架使得DFT不再只是抽象的存在性证明，而成为可以用迭代数值方法求解的计算工具。后来，大量电子结构计算软件都以这一框架为核心。
1.4 现代DFT的发展与应用扩展
随着计算机性能提升和数值算法成熟，DFT逐步从基础理论工具发展为材料模拟与分子计算的主力方法。不同类型的交换-相关泛函不断出现，涵盖从局域近似到杂化泛函、从静态体系到时间依赖体系的多种场景。
如今，DFT已广泛用于分子结构优化、晶体能带分析、表面吸附、缺陷研究、催化机理探索以及功能材料设计等领域。它以相对适中的计算成本，提供了较高的预测能力，因此在理论化学与凝聚态物理中都占据核心地位。
2 基本原理
2.1 电子密度的核心地位
在密度泛函理论中，电子密度是最基础的变量。对任意空间位置而言，电子密度表示该点附近找到电子的概率分布特征。与多电子波函数相比，电子密度只依赖三个空间坐标，形式上更为简洁，也更适合描述大体系。
理论的基本思想是：体系的基态性质并不必然需要完整波函数才能得到，只要知道电子密度，就可以推导出体系的能量、结构和其他相关信息。这种表述方式使DFT成为一种“以密度为中心”的理论。
2.2 Hohenberg-Kohn第一定理
第一定理指出，在给定电子数的条件下，体系的基态电子密度唯一决定外势，进而唯一决定哈密顿量和全部基态性质。换言之，不同的外部势对应不同的基态密度，二者具有一一对应关系。
这一结论的关键在于“唯一性”。它说明电子密度不是附属量，而是足以编码体系全部信息的基本变量。也正因为如此，DFT可以把很多物理问题转化为对密度的研究。
2.3 Hohenberg-Kohn第二定理
第二定理建立了变分原理：对于任意满足条件的试探密度，其对应的能量泛函值都不低于真实基态能量，只有真实基态密度才能使能量取最小值。这个结论为实际求解提供了方向。
在计算中，人们通常从一个初始密度出发，逐步调整，直到能量达到稳定极小值。第二定理因此成为DFT数值迭代的理论依据，也是自洽场方法的重要基础。
2.4 变分原理与基态能量
DFT的核心计算目标之一，是通过变分方式寻找最低能量对应的密度。体系的总能量可以写成电子动能、外势能、电子-电子库仑能以及交换-相关能等部分，其中后者往往需要借助近似形式来处理。
变分原理保证了，只要泛函形式和数值过程合理，就能逐步逼近真实基态。由于这一机制，DFT特别适合用于结构优化和相对稳定性比较，例如判断不同构型谁更稳定。
2.5 交换-相关泛函的概念
交换-相关泛函是DFT中最关键也最难精确处理的部分。它汇集了电子交换效应、相关效应以及动能差异中无法直接由简单项描述的那一部分。可以说，理论精度在很大程度上取决于这一项的近似质量。
由于交换-相关能没有通用的精确表达式，实际计算必须采用各种近似泛函。不同泛函在精度、适用范围和计算成本上各有特点，这也是DFT研究持续发展的重要动力。
3 Kohn-Sham方程
3.1 约化为单粒子体系
Kohn-Sham方法将相互作用电子体系映射为一个辅助的单粒子问题。每个电子在有效势中独立运动，但所有单电子轨道共同决定总电子密度。这样做的好处在于，原本难以处理的多体问题被拆解为一组相对易解的方程。
虽然这个辅助体系并不等同于真实体系，但只要构造得当，它就能重现真实的基态电子密度。也正因如此，Kohn-Sham方案成为现代DFT最常用的实现形式。
3.2 有效势的构成
Kohn-Sham方程中的有效势通常由三部分构成：外势、Hartree势和交换-相关势。它们分别描述离子核或外场对电子的作用、电子之间的经典库仑排斥，以及量子交换相关效应。
这种分解方式具有清晰的物理意义，也便于数值计算。前两项可以较明确地写出，而最后一项则由具体泛函近似来给定。
3.2.1 外势
外势通常来自原子核对电子的吸引，也可以包含外加电场、边界条件或其他外部作用。对于分子和固体而言，外势决定了体系的基本几何框架，是电子分布形成的起点。
在实际计算中，外势往往由原子核位置和电荷数确定。若采用赝势方法，则外势会被适当平滑，以减轻核附近的数值困难。
3.2.2 Hartree势
Hartree势表示电子密度产生的平均库仑排斥，它对应经典电荷分布之间的相互作用。该项可由电子密度直接计算，因此形式相对简单。
不过，Hartree势只描述平均静电作用，不包含量子交换和真正的相关效应。正因如此，还必须加入交换-相关势来修正这部分不足。
3.2.3 交换-相关势
交换-相关势是交换-相关泛函对电子密度的泛函导数，体现电子之间的量子修正。它决定了Kohn-Sham方程中最难准确给出的部分。
不同近似下，交换-相关势的具体形式会明显变化。常见泛函的差异，往往最终反映在这一势函数的细节上。
3.3 自洽场迭代过程
Kohn-Sham方程通常需要通过自洽场迭代求解。计算先从一个初始密度出发，构造有效势并求解单粒子方程，得到新的轨道和电子密度；随后用新密度更新势，再次求解，直到前后密度变化足够小。
这一过程类似反复校正的循环。若初始猜测合理、混合算法稳定，迭代通常可以收敛到一个自洽解。对于复杂体系，收敛性往往是数值计算中的关键问题之一。
3.4 本征值与轨道解释
Kohn-Sham方程的解给出一组本征值和轨道。轨道在数学上是辅助对象，用于构造电子密度；而本征值则常被近似解释为单电子能级，尤其在讨论价带、导带或分子轨道时较为常见。
需要注意的是，这些本征值并不等同于真实多体激发能。它们主要服务于密度重建和电子结构分析，因此在物理解释上应保持适当谨慎。
4 交换-相关泛函
4.1 局域密度近似（LDA）
局域密度近似假设某一点的交换-相关能，只取决于该点的局域电子密度，并把体系视为局部类似均匀电子气。它是最早、最基础的近似之一，在许多简单体系中表现稳定。
LDA通常对结构和总能计算较为可靠，但在键长、弱相互作用和电子局域性较强的问题上可能存在偏差。尽管如此，它仍是理解后续泛函发展的重要起点。
4.2 广义梯度近似（GGA）
GGA在局域密度的基础上进一步引入密度梯度信息，使泛函不仅依赖密度本身，还考虑密度变化的空间分布。相较LDA，它通常能更好描述分子键合、表面性质和化学环境变化。
GGA已成为非常常用的近似之一。许多材料计算和化学计算都会以GGA作为默认选择，再根据体系特点进行修正。
4.3 杂化泛函
杂化泛函将部分Hartree-Fock交换与DFT交换-相关项结合起来，兼顾波函数方法与密度泛函方法的优点。由于引入了精确交换成分，杂化泛函在分子反应能、能隙和局域电子结构方面通常更准确。
其代价是计算量显著增加，尤其在周期体系中更为明显。因此，杂化泛函常用于对精度要求较高但体系规模仍可控制的场合。
4.4 双杂化与高阶近似
双杂化泛函进一步引入相关能的更高阶处理，通常结合二阶微扰思想，精度往往优于普通杂化泛函。此类方法更接近高精度量子化学计算，但成本也更高。
高阶近似还包括更复杂的梯度修正、元GGA等形式。它们旨在提高对不同化学环境的普适性，但往往伴随更复杂的参数选择和数值实现。
4.5 自相互作用误差与修正
在某些近似中，单个电子会错误地与自身产生库仑作用，这被称为自相互作用误差。该误差会影响电荷局域化、离子化能和分子电荷分布等结果。
针对这一问题，研究者提出了多种修正方案，包括专门的校正泛函和改进的交换项处理。虽然无法完全消除所有误差，但可以显著提升特定体系的表现。
4.6 经验泛函与半经验泛函
经验泛函和半经验泛函通常通过实验数据或高精度结果拟合参数，以改善对某类体系的预测效果。它们在某些应用中非常实用，尤其适合快速评估和工程化计算。
不过，这类方法的适用范围较依赖训练数据，外推能力有时不如更具普适性的泛函。因而在使用时，需要明确其优化目标与适用条件。
5 计算实现
5.1 基组与表示方法
DFT计算需要把连续空间中的电子波函数或密度表示为可数值处理的形式，这就引入了基组和离散化方案。不同表示方法会影响精度、效率以及适用体系类型。
常见策略包括平面波、原子轨道和实空间网格等。它们各有优势，实际选择通常取决于体系周期性、局域性以及程序实现方式。
5.1.1 平面波基组
平面波基组适合周期性体系，如晶体和表面模型。其数学形式规整，收敛性较好，并且便于利用傅里叶变换和并行计算。
不过，平面波方法通常需要配合赝势来处理核附近的快速变化，因此更适合晶体材料和周期边界条件下的模拟。
5.1.2 原子轨道基组
原子轨道基组以局域函数为基础，接近化学直觉，适合分子和有限体系。它在小中等规模体系中通常较高效，尤其便于描述局域键合和轨道特征。
这类基组的精度取决于函数选取和级数扩展，常见问题包括基组不完备误差和基组重叠带来的数值偏差。
5.1.3 实空间网格方法
实空间网格方法直接在空间离散点上表示电子密度和势函数，概念上较直观。它便于处理复杂边界和非周期体系，也适合某些大规模并行计算。
该方法的主要挑战在于网格密度与计算成本之间的平衡。网格过粗会影响精度，过细则会带来较大开销。
5.2 赝势与全电子方法
赝势方法用一个平滑的有效势替代原子核及内层电子的详细作用，从而减少计算难度。它特别适用于平面波计算，能够显著降低所需基组规模。
全电子方法则直接处理所有电子，不对核心区域做简化，因而在描述核附近性质时更为精确。两者之间的选择，通常取决于体系精度要求与计算资源。
5.3 数值收敛与参数选择
DFT计算结果对参数较为敏感，常见包括基组截断、k点采样、网格精度以及自洽收敛阈值等。若参数设置不足，可能导致能量、力和电子结构出现明显误差。
因此，收敛测试是标准计算流程的一部分。通常需要逐步提高参数，检查关键物理量是否稳定，以确保结论可靠。
5.4 并行计算与软件实现
现代DFT软件通常支持并行计算，以应对大体系和高精度需求。并行化可以分配不同k点、平面波、原子或空间区域的计算任务，从而提升效率。
软件实现层面还涉及矩阵对角化、密度混合、快速傅里叶变换和数据管理等模块。一个成熟的程序往往不仅要正确求解方程，还要兼顾稳定性、可扩展性和用户可操作性。
6 理论扩展
6.1 时间依赖密度泛函理论
时间依赖密度泛函理论用于研究体系在外场作用下的动态响应，常见于激发态光谱、电子跃迁和超快过程分析。它将电子密度扩展到时间变量，处理对象从静态基态转向时变体系。
该理论在描述光学性质和响应函数方面具有重要价值，但其精度仍依赖于时间依赖交换-相关近似的质量。
6.2 约束密度泛函理论
约束密度泛函理论通过对电子占据、局域电荷或其他条件施加约束，研究特定态或特定过程。它常用于分析电荷转移、激发过程和反应中间态。
这一方法能够在一定程度上增强对目标态的控制，因此常被用于辅助解释复杂体系中的电子重排。
6.3 局域密度近似的局限修正
针对LDA在非均匀体系中的局限，研究者提出多种修正思路，如引入梯度信息、考虑非局域效应或采用更复杂的交换-相关表达。此类改进旨在提升对分子和低维材料的适应性。
这些修正并非简单替代，而是针对具体不足进行补强。不同方法适合不同问题，没有单一方案能够覆盖所有情形。
6.4 DFT+U方法
DFT+U是在常规DFT基础上加入局域库仑修正项的方法，常用于处理过渡金属氧化物、稀土体系等局域d或f电子较强的材料。它能改善电子局域化和磁性描述，缓解部分自相互作用误差。
该方法的核心在于引入额外参数U。虽然参数需要谨慎选取，但在许多强局域体系中，DFT+U是兼顾效率与准确性的实用方案。
6.5 轨道自由度与自旋极化处理
在某些体系中，电子的自旋和轨道自由度会显著影响能量和结构。自旋极化DFT允许不同自旋通道具有不同电子密度，因此适合描述磁性材料和开壳层体系。
若进一步考虑轨道占据的不对称性，还可更细致地分析局域电子行为。这类处理提高了理论的表达能力，也增加了计算复杂度。
7 适用范围与局限
7.1 基态性质的描述能力
DFT最擅长描述体系的基态性质，包括几何结构、总能量、电荷分布和部分振动特征。对于许多实际问题，它能在合理成本下提供较可靠结果。
因此，DFT常被视为研究稳定构型与相对能量差的首选工具。特别是在材料筛选和分子优化中，其实用性非常突出。
7.2 激发态与带隙问题
虽然DFT可以间接提供能带结构信息，但标准DFT并不天然适合精确描述激发态。常规Kohn-Sham本征值往往低估带隙，尤其在半导体和绝缘体中更为常见。
针对激发态问题，通常需要采用时间依赖DFT、杂化泛函或更高阶多体方法进行修正。也就是说，标准DFT更适合基态，而非完整激发谱。
7.3 强关联体系中的挑战
在强关联体系中，电子之间的相互作用极强，单粒子近似往往难以捕捉真实物理。此时，普通DFT可能无法正确描述绝缘性、磁性和电子局域化。
这类问题通常需要引入更专门的方法，如DFT+U或更高层次的多体理论。强关联体系也是DFT研究中最具挑战性的领域之一。
7.4 长程相互作用与色散作用
普通局域或半局域泛函对长程色散作用的描述往往不足。范德华相互作用、π-π堆积和分层材料中的弱耦合，常常因此被低估。
为改善这一点，研究中常引入色散修正或非局域泛函。相关方法使DFT在软物质、吸附和分子晶体研究中更加可靠。
7.5 近似带来的误差来源
DFT的误差主要来自交换-相关泛函近似、基组不完备、赝势误差以及数值收敛不足等方面。不同来源的误差可能相互叠加，影响最终结论。
因此，在使用DFT时需要结合体系特点、对照实验数据并进行收敛检查。对结果的解释也应保留适度谨慎，尤其在边界案例中更是如此。
8 主要应用
8.1 分子结构与反应路径
DFT广泛用于分子几何优化、键长键角预测以及反应过渡态搜索。通过比较不同构型的能量，可以分析反应路径和能垒大小。
在有机反应、配位化学和小分子催化研究中，这类应用尤为常见。它为理解反应机理提供了较为直观的电子结构视角。
8.2 晶体结构与能带计算
在固体物理中，DFT常用于晶体结构稳定性、态密度和能带结构分析。它能够帮助判断金属、半导体和绝缘体的电子特征。
通过比较不同晶型或不同原子排列的总能，研究者还能评估材料的热力学稳定性。这使DFT成为材料设计中的基础工具。
8.3 表面、界面与吸附研究
DFT特别适合研究表面吸附、界面耦合和薄膜生长等问题。由于这些体系往往具有局域电子重排和不完整配位，电子结构分析非常重要。
在吸附研究中，DFT能够计算吸附能、键合方式和电荷转移情况，从而解释分子如何与表面相互作用。
8.4 缺陷、掺杂与催化
缺陷、空位、替位掺杂和杂质态会显著改变材料性质。DFT可以分析这些微观结构对电学、磁学和化学反应性的影响。
在催化领域，DFT常用于研究活性位点、吸附中间体和反应势垒，为筛选催化材料提供理论依据。
8.5 光电材料与功能材料设计
DFT在光电材料、压电材料、储能材料和二维材料等方向都有广泛应用。通过电子结构与能量分析，可以预测材料是否适合用于发光、导电、储能或传感。
它尤其适合做材料筛选和结构优化的前期工作。虽然并非所有性质都能完全精确预测，但其趋势判断能力很强。
9 相关理论与比较
9.1 与Hartree-Fock方法的比较
Hartree-Fock方法使用单行列式波函数，并显式考虑交换作用，但通常忽略电子相关。DFT则通过交换-相关泛函统一处理交换与相关，常在计算效率上更有优势。
在许多体系中，DFT对总能和结构的表现往往更实用，而Hartree-Fock在某些局域化和交换主导问题上有其独特价值。两者各有适用范围。
9.2 与从头算量子化学方法的关系
DFT通常被视为从头算方法的一种，因为它从基本量子力学出发，不依赖经验势来描述体系主体。与传统波函数型从头算方法相比，DFT往往计算代价更低。
不过，“从头算”并不意味着完全无近似。DFT的准确性仍依赖于所选泛函，因此在严格意义上它是基于第一性原理的近似理论。
9.3 与多体微扰理论的关系
多体微扰理论从相互作用电子体系出发，逐步加入修正项，适合处理某些弱关联或准粒子问题。DFT则以密度为核心变量，更注重基态性质的高效求解。
在实际研究中，两者有时可互补使用。DFT可提供可靠的初始电子结构，而多体微扰方法可进一步修正激发性质和准粒子能级。
9.4 与分子动力学的耦合
DFT常与分子动力学结合，形成第一性原理分子动力学，用于模拟有限温度下原子运动与结构演化。这样可以在不依赖经验势的情况下研究化学反应和材料变形。
这种耦合方法精度较高，但成本也较大，因此通常适用于中等规模体系和较短时间尺度的动力学过程。
9.5 与实验结果的对应关系
DFT计算结果通常需要与实验数据相互印证，如晶格常数、振动频率、吸附能和电子谱等。对比实验不仅能检验理论可靠性，也有助于识别泛函的适用边界。
在多数情况下，DFT能较好再现实验趋势，但对高精度定量预测仍需结合更先进方法。因而它更多被视为一种强大的理论辅助工具。
10 重要人物与代表性工作
10.1 Walter Kohn的贡献
Walter Kohn是密度泛函理论的核心奠基者之一。他不仅提出了关键理论框架，还推动了DFT从概念证明走向可计算方法。
Kohn的工作使电子密度成为研究多体系统的中心变量，其影响遍及凝聚态物理、化学与材料科学。他也因此获得了极高学术声誉。
10.2 Pierre Hohenberg的贡献
Pierre Hohenberg与Kohn共同提出了奠基性的定理，为DFT提供了严格理论基础。他的贡献主要体现在从数学上证明电子密度与基态性质之间的决定性关系。
这一贡献虽然发表于较早阶段，却对后续数十年的理论发展产生了深远影响。没有这一步，后来的Kohn-Sham形式也难以建立。
10.3 Kohn-Sham方案的影响
Kohn-Sham方案是DFT能广泛应用的关键。它把抽象的密度泛函理论转化为可执行的单粒子方程，使电子结构计算变得现实可行。
今天几乎所有主流DFT程序都建立在这一思想之上。它不仅是一个数学工具，也是一种极具影响力的物理建模方式。
10.4 后续泛函发展的代表研究
在Kohn-Sham框架之后，众多研究者持续改进交换-相关泛函，从LDA到GGA、元GGA、杂化泛函，再到各种色散修正和强关联方案，形成了丰富的理论谱系。
这些工作不断提升DFT的准确度和适用范围，也让它能够应对更复杂的化学与材料问题。泛函发展至今仍是该领域的重要研究方向。
11 术语与常见概念
11.1 电子密度
电子密度是描述单位体积内电子分布的物理量。它反映了电子在空间中的聚集与扩散情况，是DFT最核心的输入和输出之一。
11.2 基态
基态指体系能量最低的量子态。DFT主要关注基态，因为Hohenberg-Kohn定理和变分原理都是围绕基态建立的。
11.3 功能泛函
功能泛函是“以函数为自变量的函数”。在DFT中，它表示某个物理量不是直接依赖数值，而是依赖整个电子密度分布。
11.4 自洽场
自洽场是指体系的输入密度与输出密度最终相互一致的状态。DFT计算通常通过反复迭代寻找这种稳定解。
11.5 交换-相关能
交换-相关能是对电子交换和相关效应的综合描述。它是DFT中最难精确求出的部分，也是各种近似泛函的核心所在。
11.6 第一性原理
第一性原理是指从基本物理定律出发，不依赖经验参数或仅依赖极少经验输入进行计算。DFT常被归入这一类方法。
12 相关软件与工具
12.1 常见DFT计算程序
常见的DFT程序包括面向分子体系、周期体系和大规模并行体系的多种软件。它们通常提供结构优化、电子结构、振动分析和能量计算等功能。
不同程序在基组类型、泛函支持、并行效率和用户界面方面各有侧重。研究者通常根据具体任务选择合适工具。
12.2 输入输出与结果分析
DFT计算的输入通常包括初始结构、泛函选择、基组参数和收敛条件等；输出则包括总能、电子密度、能带、态密度和力等信息。
结果分析往往需要结合图形和数值表格，检查结构是否收敛、能量是否稳定以及电子分布是否合理。规范的后处理是保证结论可信的重要环节。
12.3 可视化与后处理工具
电子密度、轨道和能带信息通常需要借助可视化工具进行展示。通过等值面、二维切片和能带图，研究者可以更直观地理解体系电子结构。
后处理工具还可用于分析电荷转移、局域态和键合特征，使DFT结果更便于解释与传播。
12.4 常用数据库与材料平台
随着高通量计算发展，许多材料数据库和在线平台积累了大量DFT结果，涵盖结构、能量和物性信息。它们可用于材料筛选、数据对比和计算复现。
这些平台提升了DFT研究的共享性和效率，也推动了计算材料学的发展。通过与数据库结合，研究者能够更快地定位候选材料和验证趋势。