修改第一处NMAX_DEG的赋值,默认数值为48,我这里修改为480。 然后返回安装目录,重新提交make命令进行编译。 重编译过程只与刚刚修改过的subrot_cluster.F文件相关,不会影响vasp的正常使用。 使用重编译后的vasp再次进行DFPT计算即可正常完成任务。
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...
VASP 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,准...
第二步扩胞phonopy -d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 ...
请问一下各位使用DFPT计算介电线性响应对于外场的那部分电子步不收敛,发生震荡。请问有什么好的解决办法...
然后再通过后处理得到Hessian矩阵。IBRION=5其实相当于上面的过程全部在VASP内部进行,IBRION=6用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...
我的计算采用VASP+DFPT+phonopy(IBRION=8然后用phonopy处理)超胞用的6*6*1,K点为9*9*1 这是我...
你参照phonopy的网站即可,有两种模式,一种是有限差分,一种是线性响应,你google搜索phonopy,之后到phonopy网站,找VASP+phonopy,然后按说明做即可! 如果只算gamma点的,vasp中把IBRION改成6或者8就行,你搜一下IBRION这个选项future_wl(站内联系TA)我们用vasp+phonon算的,不过phonon是收费的。
VASP 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,准...