第一步优化结构,EDIFF=1E-6或者78... 第二步扩胞 phonopy -d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 mv POSC...
运行phonopy -d --dim="2 2 2",进行扩胞(利用phonopy软件,没有的话进行安装) phonopy -d --dim="2 2 2" (="没有空格") 备份原胞文件,复制扩胞文件 cp POSCAR POSCAR-unitcell cp POSCAR-001 POSCAR 确保POSCAR是扩胞的晶格,提交vasp计算,建议用提交作业脚本。 mpirun -np 8 vasp_std #不建议用 ...
不扩胞会有虚频一般
用0.1就行了。还有就是看人家都到底是扩的多大的包。你这不是应助贴，我没法给你评分 ...
请问题主解决了吗？我算的声子谱出现了一样的问题，求教下您知道什么原因吗
超胞3*3*2, 250个原子 k点4*4*4 IBZKPT 36个 这是我的Incar SYSTEM =Sr (SrTiO3) PREC = Accurate ENCUT = 600 SYMPREC=1E-2 ISYM=0 #关闭对称性# IBRION = 8 EDIFF = 1.0e-08 IALGO = 38 ISMEAR = 0; SIGMA = 0.1 #对于金属ISMEAR=1,SIGMA=0.2 LREAL = Auto ADDGRID = .TRUE. LWAVE...
应该没问题啊...I had being in trouble for many weeks for this problem. Here is my experience/...
NSW=0 #注意这里和 DFPT 方法不同 IBRION=-1 #还有这里 EDIFFG=-1E-7 ISMEAR=0 SIGMA=0.01 ②扩胞: ~/thirdorder_vasp.py sow a b c -d 其中~是家目录的意思,这里把Thirdorder目录里面的thirdorder_vasp.py放到家目录下面,比较方便调用,当然也可以放到自己习惯的路径。
1. 原胞体系有16个原子，如果计算GAMMA点振动频率需要建立超胞吗？多大的超胞比较好？2. 从OUTCAR的...