Pony of Shadows

Queen of Truth & King of Practice


DFT-QE折腾记录

问题

作业说明: 请同学们自己寻找计算资源、自学密度泛函理论计算软件,计算单层MoS2相关性质,包括:参数测试、晶格常数测试、几何结构优化、计算态密度、能带和声子谱。(声子谱选做) 提示:对石墨烯晶格常数测试时,不需要考虑石墨烯厚度的变化;而单层MoS2是3原子层厚度,因此变xy平面内晶格常数时,MoS2厚度会变化。 请详细说明步骤,图文结合,上传pdf报告。

Install QE on Archlinux

# blas-openblas-0.3.28-1 and blas-3.12.0-5 are in conflict
sudo pacman -S blas-openblas
# Install QE from AUR
paru -S quantum-espresso

现在系统里多了一堆可执行程序:

  1. pw.x(Plane-Wave Self-Consistent Field)

    • 用于执行密度泛函理论 (DFT) 的自洽计算。它使用平面波基组和赝势方法来求解电子结构,计算材料的基态性质(如能带结构、总能量等)。
  2. ph.x(Phonon Calculations)

    • 用于计算声子(即晶体的振动模式)和相关的晶格动力学性质。适用于研究材料的热力学、声子频率、红外和拉曼散射以及超导电性等。
  3. dos.x(Density of States)

    • 用于计算态密度 (DOS)。通过态密度的计算,可以分析电子结构的分布,从而了解材料的电子性质。
  4. bands.x(Band Structure)

    • 用于计算能带结构。可以绘制材料的能带图,从而分析其导电性和带隙特性。
  5. pp.x(Post-Processing)

    • 提供了一些后处理工具,可以对密度、波函数等物理量进行处理和可视化。支持各种物理量的可视化,包括电荷密度、电子密度差分、库仑势等。
  6. neb.x(Nudged Elastic Band)

    • 用于过渡态搜索和最小能量路径 (MEP) 计算,适用于计算反应路径或扩散路径中的激活能。多用于研究分子和固体材料中的反应和迁移机制。
  7. cp.x(Car-Parrinello Molecular Dynamics)

    • 用于分子动力学模拟,基于Car-Parrinello方法。适用于在时间尺度上研究系统的动力学演变,广泛应用于液体、固体材料和表面化学等领域。
  8. td.x(Time-Dependent DFT)

    • 用于时间依赖性密度泛函理论 (TDDFT) 计算,可以计算材料的光学响应和激发态性质,适合研究光学材料和半导体的激发态动力学。
  9. q2r.x(Convert Q-Point to Real Space)

    • 将在 q 空间计算的声子频率和动力学矩阵转换为实空间,用于计算晶格的振动模式,特别适用于长波声子计算。
  10. matdyn.x(Phonon Dispersion and Thermal Properties)

    • 计算声子频散关系和材料的热力学性质(如比热、自由能等),主要依赖于从ph.x和q2r.x获得的数据。
  11. projwfc.x(Projected Wave Functions)

    • 用于投影态密度的计算,可以计算投影在特定原子或轨道上的态密度 (PDOS),对分析电子结构的本征特性很有帮助。

How to use QE

Find tutorials in quantum-espresso:tutorisals

使用Quantum ESPRESSO (QE) 进行计算时,一般会遵循以下工作流程,特别是针对电子结构计算的常见任务。具体流程可能会根据实际的研究需求有所不同,但通常包括以下主要步骤:

1. 准备输入文件

  • 主要的输入文件通常是.in文件,用于设置计算参数和结构信息。不同的计算类型(如自洽场计算、能带计算等)有特定的输入格式。
  • 典型参数包括晶胞结构、赝势类型、平面波截断能量、k点网格、电子自洽收敛条件等。

2. 赝势选择

  • Quantum ESPRESSO依赖于赝势来简化核和电子之间的相互作用。选择合适的赝势(例如PBE、USPP、PAW等)至关重要,通常赝势文件以.UPF格式存在。
  • 可以从Quantum ESPRESSO官方数据库或者其他赝势库中下载所需元素的赝势。

3. 结构优化(可选)

  • 使用pw.x进行结构优化,以确保几何结构处于最低能量态。结构优化可以包括原子位置、晶胞参数和体积优化。
  • 优化计算可分为两步:
    1. 离子弛豫:固定晶胞参数,优化内部原子位置。
    2. 晶格优化:优化晶胞参数以及内部原子结构。

4. 自洽场 (SCF) 计算

  • 使用pw.x进行自洽场计算(SCF)。在SCF计算中,求解体系的电子密度,获得材料的基态能量。
  • SCF计算结果为后续步骤提供电子密度和波函数数据。

5. 非自洽场 (NSCF) 计算(可选)

  • 为了更精确地计算能带结构和态密度,通常在SCF之后进行非自洽场 (NSCF) 计算。NSCF计算使用更密集的k点网格,能够在电子密度不变的情况下精确计算特定路径上的电子能级。

6. 态密度 (DOS) 计算

  • 使用dos.x工具对SCF或NSCF计算的结果进行态密度 (DOS) 分析。
  • 态密度可以帮助理解电子结构分布和材料的导电性。

7. 能带结构 (Band Structure) 计算

  • 使用bands.x计算能带结构。通常在对称点路径上计算能带,绘制材料的能带图,从而直观展示电子结构。
  • 计算完成后,可以用后处理工具(如gnuplot、Matplotlib等)绘制能带图。

8. 声子计算和动力学性质(可选)

  • 若研究材料的振动特性或热力学性质,可以使用ph.x计算声子频率和晶格动力学矩阵。
  • 然后通过q2r.xmatdyn.x等后处理工具计算声子色散关系和热力学性质(如比热)。

9. 后处理与数据分析

  • Quantum ESPRESSO提供了pp.xprojwfc.x等工具,可以用于数据分析和可视化,例如电荷密度、库仑势、投影态密度 (PDOS) 等。
  • 将数据导入外部软件(如VESTA、VASPkit、Matplotlib等)进行可视化分析。

10. 分子动力学模拟(可选)

  • 如果要研究系统的动态行为,可以使用cp.x工具进行分子动力学 (MD) 模拟。分子动力学模拟适合模拟高温条件下的材料行为,如熔化、扩散、界面稳定性等。

典型工作流程总结

  1. 准备输入文件 -> 2. 赝势选择 -> 3. 结构优化(如需) -> 4. SCF计算 -> 5. NSCF计算(如需) -> 6. DOS和能带结构计算 -> 7. 声子和动力学计算(如需) -> 8. 后处理和数据分析 -> 9. 分子动力学模拟(如需)

整个流程高度模块化,可以根据研究目标选择不同的步骤进行。

Slove previous problem about MoS2

A.准备输入文件和下载赝势文件

不同于VASP,Quantum ESPRESSO的所有类型的参数设定放在一个输入文件里,用标题隔开为不同的section。 major sections比如:

  • &CONTROL
  • &SYSTEM
  • &ELECTRONS

optional sections比如:

  • &IONS
  • &CELL 上面列出的sections并不是所有的,完整列表见此页面。这个官方doc页面包含了输入文件的每个section的配置参数,是重要的参考。

每个section都要用/结束。另外,如果你的输入文件用到了上面5个sections,它们得按照上文列出的顺序依次排列。

在输入文件写完了所有的sections后,还可以加上一些additional information about your molecule or the crystal system. 它们的标题不以&开头,结尾也不用/。例如:

  • ATOMIC_SPECIES
  • ATOMIC_POSITIONS
  • K_POINTS
  • CELL_PARAMETERS

下面给出一份作为参考的输入文件imput-test.in (不保证它是合理的)

&CONTROL
    calculation = 'scf'             ! 执行 自洽场 (Self-Consistent Field, SCF) 计算
    restart_mode = 'from_scratch',  ! 从头开始计算,不依赖以前的计算结果, 在初次计算或更改了基础设置(如结构、k 点、截断能量等)后,推荐使用 from_scratch,以避免先前结果的影响; 如果设置为 'restart',Quantum ESPRESSO 会尝试从之前的中间结果(如电子密度、波函数)中恢复计算。
    prefix = 'MoS2-para_test',      ! 计算输出文件的前缀。
    outdir = './output',            ! 输出文件的目录。
    pseudo_dir = './pseudo',        ! 赝势文件的目录。
    nstep = 100,                    ! 最大计算步数。
/

&SYSTEM
    ibrav = 4,                ! 晶格类型为6角晶格
    a = 3.16000e+00,          ! 晶格常数 a 的值,单位是 Bohr (约等于 1.67 Å)
    c = 2.00000e+01,          ! 晶格常数 c 的值,选择较大的值来确保不同层间的隔离,避免层间相互作用。
    nat = 3,                  ! 原子数,共 3 个原子(1 个 Mo 和 2 个 S)。
    ntyp = 2,                 ! 元素种类数,Mo 和 S。
    ecutwfc = 60,             ! 平面波基组的截断能量(单位 Ry),决定计算精度和速度。
    ecutrho = 200,            ! ecutrho 是密度截断能量(单位为 Rydberg, Ry),用于平面波基组中电子密度的截断。它决定了描述电子密度的平面波展开的截断范围。ecutrho 通常是 ecutwfc 的 4 倍或 8 倍(取决于所用赝势的类型)
    occupations = 'smearing', ! 固定占据数,适合半导体材料。
    smearing = 'gauss',       ! 使用高斯展宽,用于金属计算,这里是冗余的。
    degauss = 0.01,           ! degauss 是电子态展宽(smearing)参数,单位是电子伏特 (eV)。它定义了 smearing 的宽度,在计算费米面附近电子占据时,对电子态引入的模糊程度。
    vdw_corr = 'DFT-D',       ! 包含范德瓦尔斯修正,以更好地描述层间相互作用。
    lspinorb = .true.,        ! 启用自旋轨道耦合
    nspin = 4,                ! 启用非共线自旋
    noncolin = .true.         ! 启用非共线自旋
/

&ELECTRONS
    conv_thr = 1.0d-9         ! 自洽场迭代的收敛阈值,设为较小值保证高精度。
/

ATOMIC_SPECIES                ! 定义原子种类和使用的赝势文件。
    Mo   95.94000   Mo.rel-pbe-spn-rrkjus_psl.1.0.0.UPF
    S    32.06600   S.rel-pbe-n-rrkjus_psl.1.0.0.UPF

ATOMIC_POSITIONS (angstrom)
    S   0.002804304   0.001613133   8.429890591
    S   0.002804304   0.001613133   11.570109409
    Mo  1.582007393   0.913365735   10.000000000

K_POINTS automatic
    6 6 1  0 0 0

下面讲解这份输入文件中需要注意的点:

参考:FQE-ICP

  1. &CONTROL的calculation参数指定了这个输入文件用于进行何种计算,包括但不限于下面的可选择参数
  • scf - Single point energy calculation without changing atomic coordinates or the cell parameters, i.e., no structure relaxation.
  • nscf - Non-self-consistency calculation to extract the desired properties such as the DOS in case you need a denser electronic momentum mesh (k-mesh), based on the previous SCF calculation.
  • bands - Calculation of electronic bands along a path of high-symmetry k-points based on the previous SCF calculation.
  • relax - Optimization of the atomic coordinates to minimize the forces while keeping the cell fixed. This calculation requires section &ions. All options for a single SCF calculation apply, plus a few others.
  • vc-relax - Optimization of both atomic coordinates and the cell parameters. Structure relaxation is important to minimize forces acting on the atoms as well as to reduce the stress in the cell. This calculation requires section &cell.
  1. 一些重要的与收敛有关的参数

在进行密度泛函理论(DFT)计算时,为了获得精确的结果,必须为特定材料和计算任务选择合适的参数,并确保这些参数的收敛性。以下是对这些参数的详细理解和如何通过能量最小化优化这些参数的解释:

  • 平面波动能量截断(ecutwfc

    • ecutwfc 表示平面波基组的截断能量,也就是表示系统中电子波函数的基组复杂程度。
    • 较高的截断能量可以增加基组的数量,使得电子密度分布和总能量的计算更加精确。但更高的 ecutwfc 会增加计算的成本。
    • 优化 ecutwfc 的过程通常通过观察总能量的变化:逐步增加 ecutwfc 值,直到计算的总能量变化小于指定阈值,这样可以确保选取的 ecutwfc 在保证计算精度的同时也不会导致计算时间过长。
  • 布里渊区采样的 k 点网格(k-points

    • k 点网格用于布里渊区的采样,决定了电子态的离散化精度。
    • 增加 k 点密度可以提高计算的精度,但也会增加计算的复杂度。
    • 优化 k 点网格通常是通过测试不同的 k 点密度(如 5×5×5、7×7×7 等)观察总能量的变化,直到 k 点密度对总能量的影响不显著为止。
  • 自洽能量收敛阈值(conv_thr

    • conv_thr 是自洽迭代过程中的能量收敛标准,用于控制 SCF 迭代的精度。
    • 较小的 conv_thr 可以提高计算精度,但可能导致更多的 SCF 迭代,增加计算时间。
    • 优化 conv_thr 的过程通常通过减少其数值并观察能量的收敛性,确保在适当的计算精度下也不会增加过多的计算成本。
  1. 关于赝势文件

这次我是在QE官网下载了一些Ready-to-use PP tables。

寻找适合的赝势文件(pseudopotential files)是成功进行 Quantum ESPRESSO 计算的重要步骤。以下是如何找到和选择赝势文件的详细步骤:

  • 使用 Quantum ESPRESSO 官方赝势库 可以通过以下步骤获取赝势文件:

    • 从 Quantum ESPRESSO 项目的 QE Pseudopotentials Page 下载。
    • 或前往quantum-simulation站: PSlibrary.
    • 在库中,可以根据元素、交换-相关函数(如 PBE、LDA)、赝势类型(如 Norm-Conserving、Ultrasoft 等)筛选适合的赝势文件。
  • 选择赝势类型 赝势的类型决定了计算的精度和效率,常见的赝势类型包括:

    • Norm-Conserving Pseudopotentials (NCPP):严格遵守规范守恒,但需要更高的截断能量(ecutwfcecutrho)。
    • Ultrasoft Pseudopotentials (USPP):允许更低的截断能量,适合大系统。
    • Projector Augmented-Wave (PAW) Pseudopotentials:近似完全电子结构计算,适合高精度计算。

选择赝势类型时,需结合计算任务和系统需求,通常优先选用已被广泛验证的赝势。

  • 注意交换-相关泛函 每个赝势文件都对应于特定的交换-相关泛函(XC Functional),例如 LDA、PBE、PBEsol 等。在 Quantum ESPRESSO 的 DFT 计算中,确保选择的赝势文件与输入文件中指定的 XC 泛函一致。例如,如果使用 PBE 泛函,赝势文件也需要使用 PBE 泛函生成。

  • 确认赝势文件的版本和适用范围 下载赝势文件时,注意查看文件名中的版本信息和适用范围。例如,Mo_ONCV_PBE-1.0.upf 表示这个赝势是使用 PBE 泛函生成的 Mo 元素赝势,版本是 1.0。

  • 验证赝势文件的适用性 在某些计算(如高精度能带结构、态密度、声子谱)中,需要验证所选赝势的适用性。可以通过以下方式验证:

    • 测试收敛性:检查 ecutwfcecutrho 的值,观察计算结果(如总能量)随截断能量变化的收敛情况。
    • 对比实验值或已发表的结果:可以与已知的实验数据或文献中的理论计算结果进行对比,以验证赝势的可靠性。
  • 赝势文件的格式 Quantum ESPRESSO 支持多种赝势文件格式,常见的格式包括:

    • .UPF 格式:这是 Quantum ESPRESSO 推荐的赝势格式,包含了赝势的各类信息。
    • .PSF 格式:在早期的 Quantum ESPRESSO 版本中也被使用,但现在已不常用。
  • 存放赝势文件的位置 下载的赝势文件需要放置在 pseudo_dir 指定的目录中。例如,如果输入文件中指定 pseudo_dir = './pseudos',那么赝势文件应存放在 pseudos 目录下,以便 Quantum ESPRESSO 正确调用。

  • 常用赝势库和资源 除了 Quantum ESPRESSO 官方库,还有一些第三方赝势库可以参考,包括:

    • SSSP (Standard Solid State Pseudopotentials):为量子材料设计的高质量赝势库,适用于 PBE 和 PBEsol 泛函。
    • PseudoDojo:提供广泛的赝势文件,覆盖了多个函数和赝势类型。
    • Materials Project 和 Open Quantum Materials Database (OQMD):这些数据库中包含使用特定赝势的计算数据和资源。

通过这些资源和步骤,可以找到并验证适合具体材料系统和计算任务的赝势文件。

B. 赝势文件的选择

网站上中每个元素(比如MoS)都列出了一堆赝势文件,该怎么选择?

在选择单层 MoS₂ 的钼 (Mo) 元素赝势时,以下几个因素可以帮助你在这些选项中做出合适的选择:

1. 赝势类型

  • PAW (Projector Augmented-Wave)USPP (Ultrasoft Pseudopotentials) 是两种常见的赝势类型:
    • PAW 赝势更接近全电子计算,适合高精度计算,但可能需要较高的计算资源。
    • USPP 赝势允许更低的截断能量,适合更大规模的计算,通常效率更高,但在精度上可能略低于 PAW。
    • 如果你的计算需要高精度(例如对电子结构和态密度有高要求),推荐使用 PAW。否则,如果计算效率是主要考虑,可以选择 USPP

2. 交换-相关泛函类型

  • PBEPBESOL 是两种不同的广义梯度近似 (GGA) 泛函:
    • PBE (Perdew-Burke-Ernzerhof) 是通用 GGA 泛函,适合大多数固态材料。
    • PBESOL 是针对固态系统优化的 GGA 泛函,通常在晶格常数和结构优化方面精度更高。
    • 如果计算的目的是研究结构或晶格常数,PBESOL 会是更好的选择;如果研究电子结构或能带,PBE 更常用。

3. 相对论处理

  • 有些赝势是 Scalar Relativistic (SR),而有些是 Full Relativistic (FR)
    • Full Relativistic 赝势考虑了自旋-轨道耦合(SOC),对于重元素(如 Mo)有时会显著影响计算结果。
    • 如果你的计算涉及自旋相关的性质(如磁性或自旋极化态)或你需要高精度的能带结构,可以选择 Full Relativistic 赝势。
    • 对于一般的结构优化和基于非磁性的电子结构计算,Scalar Relativistic 已经足够。

4. 赝势版本

  • 版本 1.0.0 和 0.2 的区别在于生成的代码版本不同(v6.3 和 v5.0.2)。通常较新版本的赝势经过改进和更严格的测试,推荐使用较新的版本(1.0.0)以确保兼容性和精度。

总结

  • PAW 赝势适合高精度需求; USPP 赝势允许更低截断能量,通常效率更高。
  • PBE比较通用; PBESOL 泛函适合结构优化。
  • Full Relativistic 适合需要考虑自旋-轨道耦合的情况;相反的则是 Scalar Relativistic

选择建议(以Mo为例)

综合这些因素,针对单层 MoS₂ 计算的常见选择方案如下:

B. 尝试计算一次

这样用前面的输入文件计算一次

pw.x < input_1para-test.in > ./output/MoS2-test.out

运行结束后,可在文件./output/MoS2-test.out看到计算过程中所有的总览信息

运行结束后,可在./output目录下看到运行结果相关文件:

./output/
├── MoS2-test.out
├── MoS2-test.save
│   ├── charge-density.dat
│   ├── data-file-schema.xml
│   ├── Mo.rel-pbe-spn-rrkjus_psl.1.0.0.UPF
│   ├── S.rel-pbe-n-rrkjus_psl.1.0.0.UPF
│   ├── wfc10.dat
│   ├── wfc11.dat
│   ├── wfc12.dat
│   ├── wfc13.dat
│   ├── wfc1.dat
│   ├── wfc2.dat
│   ├── wfc3.dat
│   ├── wfc4.dat
│   ├── wfc5.dat
│   ├── wfc6.dat
│   ├── wfc7.dat
│   ├── wfc8.dat
│   └── wfc9.dat
└── MoS2-test.xml

以上的输出文件和目录结构显示了 Quantum ESPRESSO 执行计算时生成的文件。下面是每个文件和目录的含义及需要关注的部分:

1. MoS2-test.out (输出文件)

这是 Quantum ESPRESSO 的标准输出文件,包含了计算过程中所有的日志信息。您需要关注以下几部分:

  • 计算状态:是否计算成功,是否达到收敛。
  • 自洽场收敛:检查 SCF(自洽场)计算的收敛情况,确保电子密度和能量收敛到设定的阈值(通常是 conv_thr 设置的阈值)。
  • 总能量:检查计算的总能量,特别是在优化结构时是否有稳定的总能量变化。
  • 可能的错误和警告:注意是否有错误信息,尤其是与赝势、原子位置或计算参数相关的问题。

2. MoS2.save/ (保存计算的目录)

这是 Quantum ESPRESSO 为了存储计算结果和中间数据而生成的文件夹。里面存储了用于后续计算的文件,包括波函数、密度数据等。

重要文件和数据:

  • charge-density.dat

    • 存储了电荷密度分布的数据。您可以用它来分析材料的电子密度分布。
  • Mo.rel-pbe-spn-rrkjus_psl.1.0.0.UPFS.rel-pbe-n-rrkjus_psl.1.0.0.UPF

    • 这两个文件是您使用的赝势文件(Pseudopotentials)。它们描述了 Mo 和 S 元素的相对论性效应和电子结构。
  • wfcX.dat(例如 wfc1.dat, wfc2.dat, 等)

    • 这些是计算过程中生成的波函数数据文件。每个文件对应不同的波函数数据,可以用于后续的分析(例如态密度计算)。如果您需要使用这些波函数进行进一步分析(如电子态密度、能带结构等),这些文件会非常重要。
  • MoS2.xml

    • 这个文件包含了与计算相关的结构和参数信息,它是一个 XML 文件,适合用于后处理分析和可视化工具导入。可以通过一些可视化工具(例如 XCrysDenVESTA)加载这个文件来可视化计算结果。

3. 需要特别关注的内容

  • test.out 中的错误信息:如果计算没有成功完成,您需要在 test.out 中查找错误或警告消息,解决这些问题后重新进行计算。

  • charge-density.dat:如果您需要分析电荷分布或进行更详细的电荷密度分析,可以关注该文件。

  • wfcX.dat 文件:这些文件包含了波函数数据,是后续分析中可能需要的内容。如果您的计算结果需要进一步分析,例如生成能带结构或态密度图,这些文件是必要的。

  • MoS2.xml:如果您使用可视化工具来查看结果,这个 XML 文件非常重要,它包含了您计算的基本信息和结构。

总结:

  • test.out:检查计算是否成功、是否收敛,关注错误和警告信息。
  • MoS2.save/:保存计算结果,包括电荷密度、赝势文件和波函数数据。
  • charge-density.datwfcX.dat:关注这些文件,用于电荷密度分析和波函数数据分析。
  • MoS2.xml:用于后处理和可视化。

确保这些文件没有问题,尤其是 test.out 中是否有错误信息,它是您计算是否成功的重要指标。如果有错误或警告,需要及时检查并修复。

C. 如何阅读计算日志MoS2.out

我们刚才的计算得到了一个几百行的日志文件MoS2.out

我试图解读它的内容

  • 首先输出程序版本信息
  • 然后开始读取输入文件
  • 读取赝势文件
  • 计算抽样的k点及权重
  • 下面是自洽场 (SCF) 计算的过程
    • 展示自洽场计算的参数信息(G-向量数量,FFT 维度,)和动态内存估计
    • 负核电荷检查,以及显示初始电荷的值
    • 循环的自洽场 (SCF) 迭代过程
      • 每次迭代中,Quantum ESPRESSO 会更新电子密度并逐步调整系统能量,直到收敛条件达到设定的阈值。
      • 每次迭代会展示一些参数,包括total energy, estimated scf accuracy(逐步下降到设置的精度),total cpu time spent up等。其中negative rho表示负电子密度:尽管存在负值(理想情况下为 0),但值非常小,通常不会引发问题。
      • 迭代中能量逐步下降,显示出 SCF 计算逐渐稳定,这表明计算无明显数值不稳定性。从最后一次迭代来看,能量已经收敛至高精度要求。
  • End of self-consistent calculation 提示自洽场(SCF)计算结束
  • 显示k点信息:平面波数量 (PWs), 每个k点下列出了多个能带的能量值(单位为电子伏特 eV),按从低到高的顺序排列。每个数值表示特定能带在该 k 点上的能量。这次计算中每个能量值出现两次,表示该能级具有自旋简并;可以关注能量值靠近零的能带,因为费米能级(通常设为 0 eV)附近的能带决定了材料的导电性和半导体特性。
  • k点信息的最后给出了Fermi energy
  • 最后!开头的文段总结了 Quantum ESPRESSO 计算的最终结果和运行性能统计。以下是各项内容的具体含义 这一部分内容总结了 Quantum ESPRESSO 计算的最终结果和运行性能统计。以下是各项内容的具体含义:

1. 费米能级

  • 费米能级 (the Fermi energy is 0.1290 eV):这是电子的最高占据能级,标志着能带结构中电子填充的上限。费米能级的位置对于确定材料的导电性至关重要。

2. 总能量和收敛精度

  • 总能量 (total energy = -181.67467619 Ry):这是自洽场计算的最终收敛总能量,是系统的稳态能量。
  • 估计的 SCF 精度 (estimated scf accuracy < 3.7E-10 Ry):收敛精度非常高,满足了设定的收敛要求,说明计算的结果具有较高的精度。
  • 展宽贡献 (smearing contrib. (-TS)):由于采用了展宽技术(如高斯展宽),-TS 表示展宽的负熵贡献,通常在金属系统中用于处理费米能级附近的电子态分布。在此为 -0.00000000 Ry,几乎没有影响。
  • 内能 (internal energy E=F+TS):E 是自由能 F 加上展宽项 TS 后的内能,与总能量一致。

3. 能量分项

  • Quantum ESPRESSO 将总能量分解成不同物理贡献,以帮助理解系统的能量组成:
    • 单电子贡献 (one-electron contribution):描述电子的单粒子贡献。
    • Hartree 贡献 (hartree contribution):描述电子间的库伦相互作用。
    • 交换-相关贡献 (xc contribution):由交换-相关泛函(例如 PBE)描述电子的量子力学关联效应。
    • Ewald 贡献 (ewald contribution):描述了离子之间的静电相互作用。
    • 色散校正 (Dispersion Correction):用于校正长程范德瓦尔斯力,对层状材料和分子相互作用的描述非常重要。

4. 计算收敛信息

  • 收敛 (convergence has been achieved in 9 iterations):表明计算在 9 次迭代后达到自洽,符合预期的收敛标准。

5. 输出文件写入

  • Writing all to output data dir ./output/MoS2-test.save/:将计算结果(包括 XML 数据文件、电子密度、赝势和波函数)写入指定的输出目录中,为后续的分析提供了基础数据。

6. 性能统计

计算的 CPU 和墙上时间(实际运行时间)统计,为了解不同部分的计算效率提供了信息。

  • 主模块统计

    • init_runelectrons 等模块描述了初始化和电子结构计算的总体耗时。
  • 子例程统计

    • c_bandssum_band 等例程详细分解了电子结构计算中的时间消耗情况,例如波函数初始化、能带求和、密度混合、FFT 变换等。
    • fftfftw 是快速傅里叶变换(FFT)的相关模块,显示了在波函数和密度计算中 FFT 的调用次数和时间消耗。
  • 并行计算统计 (PWSCF):总 CPU 时间为 4 分 42 秒,总墙上时间为 1 分 22.9 秒,表明使用了并行计算来加速。

7. 终止时间

  • 最后记录了计算终止时间为 20: 8:22 9Nov2024

总结

这部分内容包括计算的最终收敛信息、能量分解、性能统计和文件写入信息:

  • 费米能级、总能量、收敛精度是关键的物理结果。
  • 能量的不同组成项为电子、离子相互作用的具体贡献提供了信息。
  • 性能统计显示了计算的时间消耗情况。

D. 选定参数后进行结构优化

E. 计算态密度,能带,声子谱的命令

  1. 首先要先进行自洽场计算(SCF)得到电子密度,然后进行非自洽场(NSCF)计算得到波函数信息。
pw.x < scf.in > scf.out
pw.x < nscf.in > nscf.out
  1. 在 NSCF 计算完成后, 可进行DOS计算得到态密度:
dos.x < dos.in > dos.out
  1. 在 NSCF 计算完成后,可计算能带:
bands.x < bands_post.in > bands_post.out
  1. 在在完成 SCF 计算后,使用 ph.x 程序进行声子计算
ph.x < ph.in > ph.out
# 声子计算的后处理
q2r.x < q2r.in > q2r.out
matdyn.x < matdyn.in > matdyn.out