第一步优化结构,EDIFF=1E-6或者78... 第二步扩胞phonopy-d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 mv POSCA
整理自天玑算的王二狗:天玑算-科研服务丨热电教程(4):VASP+Phonopy计算声子谱和二阶力常数矩阵_哔哩哔哩_bilibili 前面是DFPT方法,后面是有限位移法。 〇 需要准备的软件:vasp+phonopy ①准备高精度优化好的结构 POSCAR,以及 POTCAR ②使用命令 phonopy -d --dim='5 5 1',意思是扩胞成5*5*1的超胞。 ③之...
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...
(1) 首先对结构扩胞,采用DFPT法进行声子谱计算,得到的声子谱没有虚频; (2) 获得FORCE_CONSTANTS文件基础上,准备mesh.conf文件,文件内容如下; (3) 输入如下命令,即可获得thermal_properties.yaml文件; (4) 为了绘图方便,接着输入如下命令,将thermal_properties.yaml文件转换成thermal.dat文件; (5) 用Notepad++打...
phonopy-bandplot --gnuplot>PBAND.dat 来将数据存放到PBAND.dat文件中,进行后续处理。 TIPS:关于虚频,使用上述方法得到的声子谱常常在gamma点附近有虚频,解决方法我知道的一是尽可能的再去优化原胞结构,二是尝试使用更大的超包进行计算,当然那样运算量会更大。 新手一个,仅供参考。
利用VASP计算碳纳米管的能带、声子谱及声子态密度，首先需要了解前两节课所讲解的直接法与密度泛函微扰理论（DFPT）计算声子谱的基本步骤。通过第一性原理的学习，我们已经掌握了以SiO2-HP和NaCl为例的声子谱计算方法。要深入探索声子相关的其他性质，例如声子态密度、自由能、热容、熵、焓、声子速度与均方...
声子谱的计算主要有两种方法,一种是直接法,另一种是微扰密度泛函方法(DFPT)。 直接法,或称frozen-phonon方法,是通过在优化后的平衡结构中引入原子位移,计算作用在原子上的Hellmann-Feynman力,进而由动力学矩阵算出声子色散曲线。直接法的缺陷在于它要求声子波矢与原胞边界supersize正交,或者原胞足够大使得Hellmann-Feynm...
除此之外，即便对DFPT来说VASP也没办法处理声子BZ采样q点为Gamma以上的情况，得通过超晶胞的声子谱进行...
计算完成后，需准备一个band.conf文件，其中包含原子名称、扩胞倍数、高对称点坐标及对应名称。在DFPT方法中，FORCE_CONSTANTS参数应设为“READ”，遵循既定的格式编写。执行三条命令即可获得二阶力常数矩阵FORCE_CONSTANTS和声子谱数据551.dat。对于有限位移法，步骤从准备文件开始，直接接续到上述方法。需...
第一性原理||VASP计算声子谱(以SiO2-HP为例)。那么本期将为大家介绍另一种方法,除了直接法还可以应用密度泛函微扰理论(DFPT)来进行计算声子谱。 DFPT是直接计算得到force constants文件,在经过后处理得到对应的声子谱。今天,今天我们简单描述通过VASP 结合DFPT怎么计算声子谱,以NaCl为例,详细参考官网(VASP-DFPT & ...