关注：

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
/