一文讲清“DFT计算中的非自洽计算”！

非自洽计算是指在已有的自洽计算结果基础上，通过调整某些参数（如k点、电子密度或势函数）进行迭代计算的方法。与自洽计算不同，非自洽计算不需要重新从头开始计算体系的基态性质，而是利用初始自洽计算的结果作为输入，仅针对特定的物理量（如能带结构、态密度或光学性质）进行进一步的迭代计算。这种方法可以显著减少计算资源的消耗，同时保持较高的精度。

非自洽计算通常用于以下场景：

1. 能带结构计算：通过调整k点路径，计算体系的能带结构。

2. 态密度（DOS）计算：在已有的电子密度基础上，计算态密度分布。

3. 光学性质计算：基于已有的电子结构数据，计算材料的光学吸收谱。

4. 电子态投影分析：对特定原子或轨道的电子态进行详细分析。

非自洽计算通常包括以下步骤：

1. 自洽计算：首先进行自洽计算，获得体系的基态电子密度和波函数。这一步骤需要设置合适的参数，如k点密度、收敛条件等。

2. 参数调整：根据需要计算的物理量，调整输入文件中的参数。例如，对于能带结构计算，需要指定k点路径；对于态密度计算，则需要设置k点密度和投影原子。

3. 非自洽迭代：基于自洽计算的结果，进行非自洽迭代。这一步骤通常不需要重新求解薛定谔方程，而是通过修改初始条件（如电子密度或势函数）进行快速迭代。

4. 结果输出与分析：将计算结果保存并分析，例如生成能带图、态密度图或光学吸收谱。

1. 能带结构计算

在研究半导体材料时，非自洽计算常用于绘制能带结构图。例如，在VASP软件中，可以通过设置ICHAR=11参数并指定k点路径来实现非自洽计算。这种方法可以快速生成能带结构图，而无需重新进行完整的自洽计算。

2. 态密度（DOS）计算

在态密度计算中，非自洽计算可以利用已有的电子密度数据，通过调整k点密度来快速生成态密度图。例如，在Quantum ESPRESSO中，通过设置ICHAR=12参数并指定k点密度，可以高效地完成态密度计算。

3. 光学性质计算

非自洽计算还可以用于光学性质的预测。例如，在研究金属-半导体界面时，可以通过非自洽计算生成光学吸收谱，从而评估界面处的光吸收特性。

4. 电子态投影分析

在某些研究中，需要对特定原子或轨道的电子态进行详细分析。此时，可以通过非自洽计算提取投影态密度（PDOS），从而更深入地理解电子分布特性。

1. 优势

● 高效性：非自洽计算避免了重新进行完整的自洽迭代，大大减少了计算时间。

● 灵活性：可以根据需要快速调整参数，适用于多种物理量的计算。

● 适用范围广：从能带结构到态密度，再到光学性质，非自洽计算均可提供有效的解决方案。

2. 局限性

● 精度问题：虽然非自洽计算在大多数情况下精度较高，但在某些极端条件下可能会出现误差。例如，在处理强关联效应时，非自洽结果可能无法完全反映真实物理现象。

● 适用范围有限：对于某些复杂的物理问题（如多体相互作用），非自洽计算可能无法提供足够的信息。

非自洽计算是DFT计算中的重要方法之一，它通过利用已有的自洽计算结果，快速高效地完成特定物理量的迭代计算。这种方法在能带结构、态密度、光学性质等领域具有广泛应用，并显著提高了计算效率。然而，在某些复杂问题中，非自洽计算可能无法完全反映真实物理现象，因此需要结合其他方法进行综合分析。