DFT计算如何考虑温度? 四种间接引入热力学效应的计算应用

密度泛函理论（DFT）在理论上无法直接设置温度参数，本文对其原因进行了探讨，并介绍了当前主流的几种间接考虑温度效应的方法，包括从头算分子动力学（AIMD）、谐波频率分析、声子谱计算、统计热力学修正与自由能计算。

通过这些方法，可以在不同程度上引入热振动、热容、熵等热力学贡献，使计算结果更接近实验温度下的真实情况。

除此之外，文中还通过列举多个实际应用案例，说明这些方法在催化、材料热稳定性、相变、吸附热力学和电池材料等领域具有重要价值，强调了在DFT计算中引入温度效应对于提升模拟准确性和实验一致性的重要性。

引 言

密度泛函理论（DFT）是一种严格针对基态（通常等效于0K静态晶体结构）的电子结构方法。在DFT的计算框架中，电子体系处于最低能量状态，原子核的运动被视为外在参数并且通常被固定（即采用Born–Oppenheimer近似）。因此，在一般的DFT计算中没有显式的温度参数。

简单地说，DFT求解的是电子基态能量和对应的势能面，而非含热运动的自由能面。换言之，DFT计算给出了“静态晶格”的能量。即使在0K下，真实体系的原子核仍然存在量子零点振动，但对于含重元素的材料这种零点能很小，可忽略不计。然而，DFT本身并不包含原子热振动和热熵效应，只提供静态的总能和势能面。

因此，无法在标准DFT输入中直接设置温度参数。虽然在DFT计算中可以采用费米–狄拉克占据数来模拟电子的有限温度占据（通常称作电子温度或电子弛豫宽展），但这只影响电子输运性质和计算收敛，并不能模拟原子的热运动，也不等同于实际的热力学温度。真正的“温度”对应于原子核的动能分布或晶格振动激发，这需要通过分子动力学或统计热力学方法来考虑。

此外，1980年代Mermin将DFT推广到了有限温度，理论上可以最小化自由能 F=E-TS得到有限温电子体系的自洽场方程。然而，实际应用中此类有限温度泛函很少使用，主要只用于极高温度的等离子体模拟。在常规材料计算中，人们通常还是使用零温的Kohn–Sham方程，并通过其他手段间接考虑温度效应。

ProfM在Matter Modeling论坛中指出：对于晶格温度，需要同时求解电子方程和核的动力学方程；若将原子固定，那么DFT计算既不包含晶格温度也不包含零点振动。因此，可以理解为：DFT本质上求解的是给定原子排列下的电子基态能量，不包含原子的热运动自由度。这一点也是DFT在绝大多数情况下能很好描述室温性质的原因，因为势能面通常对温度变化不敏感；但对于那些存在结构相变或轻原子振动显著的体系，仅靠静态DFT计算是不够的。

间接考虑温度效应的方法

尽管DFT计算中没有直接设置温度的变量，实际研究中可以通过多种方法间接引入温度效应，以获得接近实验条件的热力学量或动力学行为。

常用的方法包括：从头算分子动力学（AIMD）、谐波频率（简正模）分析、声子谱计算、统计热力学修正以及自由能计算等。以下分别介绍这些方法及其原理：

从头算分子动力学（AIMD）：AIMD将DFT计算嵌入分子动力学模拟，以模拟给定温度下原子的热运动。在AIMD中，通过赋予原子初速度并使用恒温（如Nosé–Hoover）或恒能（NVE）系综演化体系，轨迹自动采样了不同构型。

在每一步，电子结构由DFT计算决定，原子在势能面上演化。这样可以直接获得温度下的平均结构、扩散系数、动力学过程等信息。

例如，郑朝阳等综述了AIMD在含能材料中的应用，指出可以通过设置不同温度分析晶体结构、热力学、力学等性质的温度效应。

DOI: 10.11858/gywlxb.2015.02.001

Novko等人通过对Ag(111)表面吸附氮原子的AIMD模拟，比较了80K和700K下吸附概率和动力学，发现温度对吸附路径和能量耗散有细微影响。

总之，AIMD可以在模拟中人为设定体系温度，从而直接考察温度效应，但计算开销很大，且一般只能采样较小体系和有限时间。

谐波频率（简正模）分析：对于分子或有限体系，可以在DFT优化的基态构型上进行小扰动，求解哈密顿量的二阶导数矩阵来获得振动频率。

每个振动模式类似于量子谐振子，可以计算其零点能和热力学贡献。在统计热力学框架下，振动模式的配分函数

从而得到振动自由能

熵 Sv等。简单地说，谐波频率分析可以将DFT静态能量修正为包含振动贡献的热力学量。Tyberius等人提到：在谐波近似下解决振动问题可以给出基态的零点能，这对轻原子尤为重要。

Galimberti和Sauer等利用DFT计算得到甲烷、乙烷等分子在沸石中的振动模式，将这些频率投影到振动密度状态中，计算出配分函数并获得吸附的吉布斯自由能，精确地反映了温度效应。

这一类方法通常假设谐波近似成立，忽略高阶非线性耦合，但对于许多固、液、气体系的热力学量已经非常实用。

DOI: 10.1021/acs.jctc.1c00519

声子谱计算（晶格动力学）：对于周期性固体，可以利用密度泛函微扰理论或有限位移法计算声子谱。求解晶格动力学本征值问题可得到各声子模式频率 $\omega_{q\nu}$，并据此构造声子分布。

声子频率可用来计算晶体的振动内能、振动自由能、热容、热膨胀系数等热力学性质。

Togo和Tanaka指出：首先用DFT获得不同体积下的静态能量 E(V)，然后计算各体积对应的声子频率得到声子自由能Fph(T,V)；再通过自由能最小化

得到吉布斯自由能和热力学平衡体积。

他们利用铝（Al）的例子演示了这种计算：通过对不同晶格常数下求声子谱并拟合状态方程，获得了随温度变化的体积、热容等，并与实验良好吻合。

需要注意的是，这通常采用谐波或准谐波近似，即假设频率随温度很小变化，因此在接近熔点时可能失效。不过，声子谱计算是评价材料相变温度、热稳定性、界面热阻等性质的重要工具。

DOI: 10.1016/j.scriptamat.2015.07.021

统计热力学修正和自由能计算：基于以上振动计算结果，可以对DFT静态能量进行热力学修正。具体地，将DFT给出的静态能作为内能基准，再加上零点能、振动热容和熵等贡献，得到焓或吉布斯自由能。

例如，在气相分子体系中，常用理想气体近似结合分子振动和转动贡献求导热焓H(T)和熵S(T)。在固体中，可用声子自由能得到Helmholtz自由能F(T,V)，进而通过G=F+PV得到吉布斯自由能。简单公式如在谐波近似下：

其中E0为DFT静态能。通过这些统计修正，可从DFT能量推导出温度相关的物性。Galimberti等的方法也属于这一类：他们通过DFT-MD生成振动密度态，在该基础上构造配分函数，计算吸附的自由能和熵，成功重现了不同分子在实验条件下的吸附热力学变化。

综合而言，上述方法可以分别或组合使用来引入温度效应。AIMD可直接模拟热动力学过程；频率分析和声子计算可得到热力学修正；统计修正和自由能计算则将这些贡献整合到最终的热力学量中。

实际研究中，往往首先用DFT计算基态能，然后加上振动自由能等修正来获得接近真实条件的吉布斯自由能或反应自由能。

代表性研究与应用案例

表面吸附与化学反应：Novko等在Ag(111)表面模拟了氮原子的吸附过程，并通过对比80K和700K下的AIMD轨迹，发现表面温度对吸附几率、吸附后能量损失等影响有限。此类研究说明了利用恒温AIMD可以直接考察温度对表面化学动力学的作用。

DOI:10.1103/PhysRevB.96.085437

此外，Galimberti等通过DFT-MD计算得到的振动频率精确评估了甲烷、乙烷在沸石H-CHA中的吸附自由能，揭示了温度下吸附热力学的重要性。

DOI: 10.1021/acs.jctc.1c00519

晶体结构相变和热力学性质：基于声子计算的研究也很多。例如，Togo和Tanaka的工作利用声子自由能计算了Al的热膨胀系数和比热，其计算结果与实验值高度吻合。

研究者常用声子谱分析钙钛矿、金属及半导体材料的相变温度，通过自由能曲线判断高温相的稳定性。

在钛金属中，不同相的声子计算揭示了BCC相在低温下不稳定（存在虚频），只有在高温才通过振动振幅扭曲到稳定结构，这解释了实验中1100K以上才出现α→β相变的现象。

DOI: 10.1016/j.scriptamat.2015.07.021

化学热力学与电池材料：在电化学领域，研究人员利用DFT结合振动贡献计算锂离子电池材料的开路电压温度依赖。例如，通过计算锂插层与脱层反应的自由能，可以获得温度对电池电势的影响。

以上实例说明，AIMD、谐波/声子计算及统计修正等方法能够将温度引入DFT计算，在催化、能源材料、固体物理等研究领域发挥重要作用。

研究者通过这些方法可以更准确地预测材料的高温结构、热稳定性、相变行为以及化学反应的自由能变化，从而指导实验设计和工艺优化。

结 论

综上所述，标准DFT计算不包含显式温度参数，因为它求解的是给定原子配置下的基态能量。

要考虑真实条件下的温度效应，研究者需要借助额外方法：AIMD模拟可直接采样热运动；频率分析和声子谱计算则可获得振动自由能和熵；统计热力学修正和自由能计算将这些贡献整合，从而得出体系的高温热力学性质。

这些方法在实践中常常结合使用，以便获得接近实验的预测结果。各类材料和反应体系中已有大量成功案例表明，通过上述技术可以有效补偿DFT的静态近似，使计算结果与实验温度下的现象相一致。

在报道或使用DFT结果时，关注温度效应并进行必要的修正，是提高模拟可靠性和与实验对比意义的关键。