QE软件支持哪些计算类型？

Question

QE（Quantum ESPRESSO）软件支持哪些常见的第一性原理计算类型？能否详细说明其在电子结构计算、几何优化、分子动力学模拟、能带结构与态密度计算，以及自洽场和非自洽场计算中的具体应用？此外，是否支持如范德华修正、自旋轨道耦合、杂化泛函等高级功能？这些计算类型在实际材料模拟中如何选择与配置？

Answer 1

一、Quantum ESPRESSO 支持的常见第一性原理计算类型

Quantum ESPRESSO（QE）是一款基于密度泛函理论（DFT）的开源第一性原理计算软件包，广泛应用于材料科学、凝聚态物理和化学等领域。其核心模块包括 pwx（pw.x）、phx（ph.x）、ppx（pp.x）等，支持多种类型的电子结构与材料模拟任务。

1. 电子结构计算

• 电子结构计算是 QE 的基础功能，通过求解 Kohn-Sham 方程获得体系的基态电子密度与总能。
• 使用 pwx 模块执行自洽场（SCF）迭代，收敛后可输出电荷密度、费米能级、能级分布等信息。
• 适用于绝缘体、半导体、金属等多种材料体系。
• 输入参数中需设置 calculation = 'scf'，并配置截断能（ecutwfc）、k点网格（K_POINTS）等关键参数。

2. 几何优化（结构弛豫）

几何优化用于寻找原子在势能面上的局部极小值，即稳定构型。

优化方法	对应关键词	适用场景
BFGS	ion_dynamics = 'bfgs'	大多数固体与分子系统
共轭梯度	ion_dynamics = 'cg'	复杂势能面或初态远离平衡时
分子动力学驱动	ion_dynamics = 'damp'	轻元素或软模系统

3. 分子动力学模拟（AIMD）

QE 支持基于 Born-Oppenheimer 近似的从头算分子动力学（ABINIT-MD），通过实时更新电子结构来驱动离子运动。

&IONS
ion_dynamics = 'verlet',
dt = 5.0,
/
&CONTROL
calculation = 'md',
/

可用于研究热力学性质、相变过程、扩散行为等动态演化问题。

4. 能带结构与态密度计算

在完成 SCF 计算后，可通过非自洽场（NSCF）计算获取高精度能带路径上的本征值。

1. 步骤一：执行 SCF 计算，生成收敛的电荷密度。
2. 步骤二：设置 calculation = 'nscf'，定义精细 k-path（如 Γ→X→M→Γ）。
3. 步骤三：使用 bands.x 工具提取能带数据。
4. 步骤四：调用 pp.x 计算态密度（DOS）与投影态密度（PDOS）。

5. 自洽场（SCF）与非自洽场（NSCF）计算

SCF 是所有计算的基础，目标是使电荷密度在迭代中收敛；NSCF 则用于固定电荷密度下进行单次对角化，常用于光学响应、能带绘制等后处理任务。

graph TD A[开始] --> B[初始化波函数] B --> C{是否收敛?} C -- 否 --> D[更新电荷密度] D --> E[求解Kohn-Sham方程] E --> C C -- 是 --> F[输出收敛电荷密度]

6. 高级功能支持情况

QE 不仅支持标准 LDA/GGA 泛函，还集成了多项高级物理效应建模能力：

• 范德华修正：通过 DFT-D2/D3 方法（Grimme 修正）或 vdW-DF 非局域泛函实现，适用于层状材料、吸附体系。
• 自旋轨道耦合（SOC）：需启用非共线磁性（noncolin = .true.）与旋轨耦合项，在重元素体系中尤为重要。
• 杂化泛函：通过 EXX 模块支持 HSE06、PBE0 等混合泛函，显著提升带隙预测精度，但计算成本较高。
• GW 与 Bethe-Salpeter：通过 epsilon.x, sigma.x, bse.x 实现准粒子修正与激子效应分析。

7. 实际材料模拟中的选择与配置策略

在实际应用中，应根据材料类型与研究目标合理配置计算流程：

目标	推荐方法	注意事项
金属稳定性分析	GGA+smearing	使用 degauss > 0.01 提高收敛性
二维材料带隙	HSE06 或 GW	避免 PBE 对带隙低估
表面吸附能	PBE-D3	必须包含范德华力
拓扑材料	PBE+SOC	开启非共线自旋
相变路径搜索	NEB via path.x	需多节点并行加速
光学吸收谱	BSE	依赖 GW 准粒子能带
铁电极化	Berry Phase (pwcond.x)	需超胞与细致k网
热导率预测	Phonon + BTE	结合 ph.x 与 q2r.x
缺陷形成能	Supercell + Charged Correction	注意尺寸效应与静电校正
催化活性位点	DFT+U 或 Hybrid	过渡金属d电子强关联