第一步优化结构,EDIFF=1E-6或者78... 第二步扩胞phonopy -d --dim="4 4 1"这里说明一下,大致就是要让扩胞后大于10埃量级,比如我算的是体系如果是3.8埃,乘4大于10埃,这里会生成很多POSCAR,鉴于我们这个是密度泛函微扰理论,所以那些POSCAR-001之类的没用,详见另一种方法。 第三步准备文件计算 mv POSCAR...
不扩胞会有虚频一般
运行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 #不建议用 ...
而DFPT原则上不需要扩胞，但因为vasp原因，计算DFPT还是需要扩胞，请问我了解的地方有错误的地方吗？
楼主，DFPT计算是LREAL=A, 计算的结果对吗？ 我看到很多例子为了追求计算精度都是LREAL=.F ...
计算的是SiP-monolayer扩胞扩了521,虚频减小了但是剩下的0.几怎么也消不掉,如图img]img]但是OUTCAR里的f/i不是很大358f/i=0.000916THz0.0057572PiTHz0.030562cm-10.003789meV359f/i=0.001395THz0.0087672PiTHz0.046543cm-10.005771meV360f/i=0.002415THz0.0151742PiTHz0.0
前面是DFPT方法,后面是有限位移法。 〇 需要准备的软件:vasp+phonopy ①准备高精度优化好的结构POSCAR,以及POTCAR ②使用命令phonopy -d --dim='5 5 1',意思是扩胞成5*5*1的超胞。 ③之后会生成一堆文件,其中有POSCAR-00*和SPOSCAR,分别对应不同的方法密度泛函微扰法(DFPT)以及有限位移法。
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放到家目录下面,比较方便调用,当然也可以放到自己习惯的路径。
应该没问题啊...I had being in trouble for many weeks for this problem. Here is my experience/...
用0.1就行了。还有就是看人家都到底是扩的多大的包。你这不是应助贴，我没法给你评分 ...