第一性原理||VASP计算声子谱(以SiO2-HP为例)。那么本期将为大家介绍另一种方法,除了直接法还可以应用密度泛函微扰理论(DFPT)来进行计算声子谱。 DFPT是直接计算得到force constants文件,在经过后处理得到对应的声子谱。今天,今天我们简单描述通过VASP 结合DFPT怎么计算声子谱,以NaCl为例,详细参考官网(VASP-DFPT & ...
使用密度泛函微扰方法计算材料力常数时,由于使用超胞结构和高精度的设置等等,计算时会产生诸多报错,其中exceeds NMAX_DEG则是一个较为难解决的问题。 笔者以往写过一则SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG报错提示的解决办法,即弃用密度泛函微扰,改用有限位移方法。DFPT计算报错internal error in SETUP_DEG_...
DFPT的两种主要形式是由Baroni和Gonze提出的，虽然这两种形式可以被证明是等价的，但在实现上存在差异，这可能导致一种方法优于另一种方法。Baroni形式的核心是获得一系列方程，这些方程可以用Green函数方法自洽地求解；Gonze形式是基于Kohn-Sham能量泛函的扰动展开，导致了类似于零阶问题的偶数级展开的变分问题。一些要点...
2、声子谱计算方法 声子谱的计算主要有两种方法,一种是直接法,另一种是微扰密度泛函方法(DFPT)。 直接法,或称frozen-phonon方法,是通过在优化后的平衡结构中引入原子位移,计算作用在原子上的Hellmann-Feynman力,进而由动力学矩阵算出声子色散曲线。直接法的缺陷在于它要求声子波矢与原胞边界supersize正交,或者原胞足够...
考虑了微扰理论在DFT中的应用,并用这种形式推导出可以计算晶体材料中声子和电场响应的方程。 DFPT的两种主要形式是由Baroni和Gonze提出的,虽然这两种形式可以被证明是等价的,但在实现上存在差异,这可能导致一种方法优于另一种方法。Baroni形式的核心是获得一系列方程,这些方程可以用Green函数方法自洽地求解;Gonze形式是...
往期推文中我们已经详细介绍了运用直接法来计算声子谱,具体可点击链接进行查看第一性原理||VASP计算声子谱(以SiO2-HP为例)。那么本期将为大家介绍另一种方法,除了直接法还可以应用密度泛函微扰理论(DFPT)来进行计算声子谱。
第一步:精确优化结构 EDIFF =1E-8 EDIFFG=1E-5 增大Ecut 第二步:拿第一步优化出来的结构命名为POSCAR-unitcell计算声子谱 INCAR设置:IBRION = 8 NSW=1 第三步:扩胞 phonopy -d --dim="2 2 2" -c POSCAR-unitcell mv SPOSCAR POSCAR 以上一步生成的 SPOSCAR 为 POSCAR,准备 KPOINTS、POTCAR 和 INC...
DFPT的两种主要形式是由Baroni和Gonze提出的,虽然这两种形式可以被证明是等价的,但在实现上存在差异,这可能导致一种方法优于另一种方法。Baroni形式的核心是获得一系列方程,这些方程可以用Green函数方法自洽地求解;Gonze形式是基于Kohn-Sham能量泛函的扰动展开,导致了类似于零阶问题的偶数级展开的变分问题。
第一步:精确优化结构 EDIFF =1E-8 EDIFFG=1E-5 增大Ecut 第二步:拿第一步优化出来的结构命名为POSCAR-unitcell计算声子谱 INCAR设置:IBRION = 8 NSW=1 第三步:扩胞 phonopy -d --dim="2 2 2" -c POSCAR-unitcell mv SPOSCAR POSCAR 以上一步生成的 SPOSCAR 为 POSCAR,准备 KPOINTS、POTCAR 和 INC...
VASPDFPT 方法计算声子谱 第一步:精确优化结构 EDIFF =1E-8 EDIFFG=1E-5 增大 Ecut 第二步:拿第一步优化出来的结构命名为 POSCAR-unitcell 计算 声子谱 INCAR 设置:IBRION = 8 NSW=1 第三步:扩胞 phonopy -d --dim="2 2 2" -c POSCAR-unitcell mv SPOSCAR POSCAR 以上一步生成的 SPOSCAR 为 PO...