使用密度泛函微扰方法计算材料力常数时,由于使用超胞结构和高精度的设置等等,计算时会产生诸多报错,其中exceeds NMAX_DEG则是一个较为难解决的问题。 笔者以往写过一则SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG报错提示的解决办法,即弃用密度泛函微扰,改用有限位移方法。DFPT计算报错internal error in SETUP_DEG_...
第一步优化结构,EDIFF=1E-6或者78... 第二步扩胞phonopy-d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 mv POSCAR ...
34:20 VASP的BoltzTraP安装 09:52 VASP视频教程-重新编译以进行dfpt计算 02:50 VASP视频教程-DFPT计算 09:06 VASP视频教程-boltz (BoltzTraP) 12:01 VASP视频教程-DFPT计算 奕星模拟 164 0 VASP视频教程-输入输出提交计算任务 奕星模拟 1217 1 VASP视频教程-电荷局域分析 奕星模拟 1937 0 ...
本文将简单介绍在DFPT计算声子谱的基础上,如何通过Phonopy后处理得到自由能、热容和熵等热学性质相关数据并绘图。 操作步骤 (1) 首先对结构扩胞,采用DFPT法进行声子谱计算,得到的声子谱没有虚频; (2) 获得FORCE_CONSTANTS文件基础上,准备mesh.conf文件,文件内容如下; (3) 输入如下命令,即可获得thermal_properties....
DFPT的两种主要形式是由Baroni和Gonze提出的,虽然这两种形式可以被证明是等价的,但在实现上存在差异,这可能导致一种方法优于另一种方法。Baroni形式的核心是获得一系列方程,这些方程可以用Green函数方法自洽地求解;Gonze形式是基于Kohn-Sham能量泛函的扰动展开,导致了类似于零阶问题的偶数级展开的变分问题。一些要点...
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...
请问一下各位使用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...
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/...
请问一下各位使用DFPT计算介电线性响应对于外场的那部分电子步不收敛,发生震荡。请问有什么好的解决办法...