(1)需要提前准备文件:compile_converter.exe、converter.f90、qgroup.f90、phonon_sym.f90、diel_sym.f90、fort.224、fort.114、fort.115、fort.201 (2)检查compile_converter.exe和converter.f90路径;检查fort.114 4、ph计算: cd ../../tmp2/ cp -r scf/ ph1 cp -r scf/ ph2 cd ../work2/full ph...
1. 计算phonon需要电子自洽计算的信息,所以首先要在当前目录提交一个scf的任务。其中ph.x的输入文件中要加上recover=.true. 输出文件包括output_0文件夹,P0.scf.out和P0.dyn0 #!/bin/bash #PBS -N titanite-ph #PBS -l nodes=node06:ppn=16 #PBS -l walltime=3000:00:00 #PBS -e error #PBS -o...
(6)finally, put the phono.dat into origin, matlab, MS excel and etc. Reference: phonon manual.
DFPT 今天主要来聊一聊QE是怎么通过DFPT理论计算声子的。 在DFPT之前,人们计算声子主要有两种方法:直接超胞方法(frozen phonon)和介电方法。 超胞方法的好处在于非常直接,而且可以处理线性效应和非线性效应。而且只需要知道价态的信息就可以了(介电方法需要知道导带的信息),所以自冾计算可以使用效率比较高的自冾迭代...
最后用phonopy处理,发现自由能随温度变化差别很大。 产生这个不同的原因是什么? 有没有可能算出来本来就会一样? 附注,我的计算首先都是把结构优化了好几遍才做phonon计算的。另外对称性用的默认的ISYM=2。 这个不同是不是也跟对称性很有关系? 恳请各位大神出现指点一下,谢谢...
你参照phonopy的网站即可，有两种模式，一种是有限差分，一种是线性响应，你google搜索phonopy，之后到phonopy网站，找VASP+phonopy，然后按说明做即可！ 如果只算gamma点的，vasp中把IBRION改成6或者8就行，你搜一下IBRION这个选项future_wl(站内联系TA)我们用vasp+phonon算的，不过phonon是收费的。
首先，你这个frequency band 和phonon dos对不上，如果band上有虚频，会在dos里表现出来。band的结果，...
In particular, all the pronounced peaks of the phonon DOS curve are located at positive frequencies ranging from 0.36 to 8.34 THz, confirming the dynamical stability of the titled crystal. Further, when we consider the influ- ence of the hydrostatic pressure in the present calculations, we find...
C. Phonon spectrumD. Thermodynamics propertiesOrthorhombic CaSnO3.Density-functional theory (DFT) within the LDA and GGA approaches.Electronic band structure.Vibrational normal modes.Raman and Infrared spectra.doi:10.1016/j.jpcs.2014.09.016Moreira, E....
装好phonopy之后，你参照phonopy的网站即可，有两种模式，一种是有限差分，一种是线性响应，你google搜索phonopy，之后到phonopy网站，找VASP+phonopy，然后按说明做即可！ 如果只算gamma点的，vasp中把IBRION改成6或者8就行，你搜一下IBRION这个选项future_wl(站内联系TA)我们用vasp+phonon算的，不过phonon...