关注:
1)只谈问题
问题:
1) PWSCF 自洽计算脚本
2)QE优化计算脚本
3)QE声子谱计算脚本
4)QE能带结构、态密度计算
5) QE电声耦合常数计算脚本
6)QE电声耦合常数后处理脚本:获取超导温度
7) PWSCF如何有效建立晶体模型及查看建立模型的正确性
具体问答
7) PWSCF如何有效建立晶体模型及查看建立模型的正确性
http://blog.sciencenet.cn/blog-567091-725337.html
如设置ibrav=0,则依靠CELL_PARAMETERS指定晶格类型及所属对称性
http://web.mit.edu/xcrysden_v1.4.1/www/XCRYSDEN/doc/pwscf.html
ibrav = 0
#### ibrav=0, celldm(1) = 1.8897268777743552, ####直接将cell parameters中的A单位转化为Bohr
CELL_PARAMETERS{ alat | bohr | angstrom }
v1(1) v1(2) v1(3)
v2(1) v2(2) v2(3)
v3(1) v3(2) v3(3)
最好不要直接利用VASP优化后的坐标系及原子坐标直接建PWscf结构模型;因为VASP优化可能是打破对称性优化的,因此,其给出的坐标系是不带对称性的。
如果直接利用该坐标系建立QE结构模型,则:
1) QE优化非常耗时;
2)优化得到的结果可能很诡异。
因此,正确的做法是:
将VASP优化后的结果,导入MS中施加对称性后,再通过vesta转化为VASP格式的坐标文件。此时,即可利用重新调整后的坐标系及原子坐标,进行QE建模。
&system
ibrav=0, celldm(1) = 1.8897268777743552, ####直接将cell parameters中的A单位转化为Bohr
nat=14,
ntyp=2,
ATOMIC_POSITIONS (crystal) (VASP中的原子坐标直接引用)
H 0.000000000 0.500000000 0.750000000
H 0.000000000 0.500000000 0.250000000
H 0.750000000 0.000000000 0.500000000
H 0.250000000 0.000000000 0.500000000
H 0.500000000 0.750000000 0.000000000
H 0.500000000 0.250000000 0.000000000
H 0.500000000 0.000000000 0.250000000
H 0.500000000 0.000000000 0.750000000
H 0.000000000 0.750000000 0.500000000
H 0.000000000 0.250000000 0.500000000
H 0.750000000 0.500000000 0.000000000
H 0.250000000 0.500000000 0.000000000
X 0.500000000 0.500000000 0.500000000
X 0.000000000 0.000000000 0.000000000
CELL_PARAMETERS (VASP中的坐标矢量直接引用)
3.2657001019 0.0000000000 0.0000000000
0.0000000000 3.2657001019 0.0000000000
0.0000000000 0.0000000000 3.2657001019
ecutwfc=60.0, ecutrho = 1000.0
occupations='smearing',
smearing='methfessel-paxton',
degauss=0.05
/
转载本文请联系原作者获取授权,同时请注明本文来自叶小球科学网博客。
链接地址:https://wap.sciencenet.cn/blog-567091-1088463.html?mobile=1
收藏