第一步优化结构,EDIFF=1E-6或者78... 第二步扩胞phonopy -d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 mv POSCAR...
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...
第一步:精确优化结构 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...
第一步:精确优化结构 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...
(2)vasp+phonopy计算动力学矩阵有有限位移法和DFPT两种方法,我了解到有限位移法要求声子波矢与原胞边界...
也有文献是用DFPT算的。算声子的INCAR已贴。单独的那个声子谱是提高精度之后的结果,即力收敛是EDIFFG=-...
你参照phonopy的网站即可,有两种模式,一种是有限差分,一种是线性响应,你google搜索phonopy,之后到phonopy网站,找VASP+phonopy,然后按说明做即可! 如果只算gamma点的,vasp中把IBRION改成6或者8就行,你搜一下IBRION这个选项future_wl(站内联系TA)我们用vasp+phonon算的,不过phonon是收费的。
我的计算采用VASP+DFPT+phonopy(IBRION=8然后用phonopy处理)超胞用的6*6*1,K点为9*9*1 这是我...
没事,LDA和PBE计算的频率有些差别很正常。另外你可以用IBRION=6算声子,其他和IBRION=8是一样的。好...
vasp是不能直接计算声子谱的,至少我不知道啦, 我们经费不多,所以是用的phonopy接口vasp计算的,当然也还在摸索中,不知道其他大神有没有好的方法,或者用pwscf计算声子刘仕晨(站内联系TA)集群装phonopy可能有点麻烦了,因为装之前要装很。