使用密度泛函微扰方法计算材料力常数时,由于使用超胞结构和高精度的设置等等,计算时会产生诸多报错,其中exceeds NMAX_DEG则是一个较为难解决的问题。 笔者以往写过一则SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG报错提示的解决办法,即弃用密度泛函微扰,改用有限位移方法。DFPT计算报错internal error in SETUP_DEG_...
前面是DFPT方法,后面是有限位移法。 〇 需要准备的软件:vasp+phonopy ①准备高精度优化好的结构POSCAR,以及POTCAR ②使用命令phonopy -d --dim='5 5 1',意思是扩胞成5*5*1的超胞。 ③之后会生成一堆文件,其中有POSCAR-00*和SPOSCAR,分别对应不同的方法密度泛函微扰法(DFPT)以及有限位移法。 ④把原本的单...
本文将简单介绍在DFPT计算声子谱的基础上,如何通过Phonopy后处理得到自由能、热容和熵等热学性质相关数据并绘图。 操作步骤 (1) 首先对结构扩胞,采用DFPT法进行声子谱计算,得到的声子谱没有虚频; (2) 获得FORCE_CONSTANTS文件基础上,准备mesh.conf文件,文件内容如下; (3) 输入如下命令,即可获得thermal_properties....
DFPT的两种主要形式是由Baroni和Gonze提出的，虽然这两种形式可以被证明是等价的，但在实现上存在差异，这可能导致一种方法优于另一种方法。Baroni形式的核心是获得一系列方程，这些方程可以用Green函数方法自洽地求解；Gonze形式是基于Kohn-Sham能量泛函的扰动展开，导致了类似于零阶问题的偶数级展开的变分问题。一些要点...
密度泛函微扰理论(DFPT)是一种特别强大和灵活的理论技术,它允许在密度泛函框架内计算这些性质,从而有助于理解这些过程背后的微观量子力学机制,同时也为理论发展提供了严格的试验场。系统对外部扰动的响应可使用DFT计算,并添加一些扰动势;然而,由于此类方法涉及通过在外部扰动强度变化时进行的一系列单点能量计算来获得系统...
本文介绍VASP结合Phonopy计算声子谱和声子态密度的过程，包括相关软件安装、计算流程及示例。首先，声子谱的计算主要有两种方法：有限位移法和微扰密度泛函方法（DFPT）。有限位移法在优化后的平衡结构中引入原子位移，计算作用在原子上的力，进而得到声子色散曲线。然而，该方法在复杂系统中使用受限，尤其是对于...
VASP+DFPT+Phonopy 计算声子过程 以Ammonia为例,基本翻译自http://atztogo.github.io/phonopy/vasp.html#how-to-run。 1. 结构优化. 输入文件如下: INCAR !System settings SYSTEM = Ammonia PREC = accurate low, normal, accurate ISTART = 0 job : 0-new 1-cont 2-samecut...
一、声子态密度与投影声子态密度的获取 1. **声子态密度**：使用mesh.conf文件准备计算环境，先通过直接法计算得到FORCE_SETS文件，随后运用DFPT计算FORCE_CONSTANTS文件。执行相关命令后，可以得到total_dos.dat文件，通过作图直观展示声子态密度。2. **投影声子态密度**：同样通过mesh.conf文件，准备计算...
DFPT法需要使用SPOSCAR进行计算(单个任务)。可参考王宁博士在B站的视频[7]。笔者在这里贴出自己的代码仅供参考。 计算前需要将原高精度优化的POSCAR 重命名为POSCAR-unitcell SPOSCAR命名为POSCAR。 #INCAR for DFPTGlobal Parameters ISTART=0(Read existing wavefunction;ifthere)ICHARG=2(Non-self-consistent: GGA/...
VASP+DFPT+Phonopy 计算声子过程 以Ammonia为例,基本翻译自http://atztogo.github.io/phonopy/vasp.html#how-to-run。 1. 结构优化. 输入文件如下: INCAR !System settings SYSTEM = Ammonia PREC = accurate low, normal, accurate ISTART = 0 job : 0-new 1-cont 2-samecut...