QE动力学矩阵文件的主要内容及单位

matdyn0

不可约q点的位置

matdyn1

PHonon/PH/io_dyn_mat_old.f90 写文件部分代码
Mould/cell_base.f90
Mould/ions_base.f90 用到的变量及单位
基本形式

Dynamical matrix file
Electron-phonon coefficients for hydride
  2    4  0  1.8897269  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000
Basis vectors
     -1.551398696    1.551398696    1.551398696
      1.551398696   -1.551398696    1.551398696
      1.551398696    1.551398696   -1.551398696
           1  'H   '    918.58996499058958
           2  'S   '    29219.990958518523
    1    1      1.5513986956      0.0000000000     -0.0000000000
    2    1      0.0000000000      1.5513986956     -0.0000000000
    3    1     -0.0000000000      0.0000000000      1.5513986956
    4    2     -0.0000000000      0.0000000000      0.0000000000

注释

#固定的一行
Dynamical matrix file
#ph.in的头一行
Electron-phonon coefficients for Al
#原子类型 原子数 ibov celldm(6)单位bohr ibov是0的还有一行
#Basis vectors (晶格坐标，以celldm(1)为单位)

  1    1  2  7.5000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000


#（原子质量单位 a.m.u）
           1  'Al  '    24590.765652728711
#    原子指标 原子指标 原子位置（单位Bohr 直角坐标）
    1    1      0.0000000000      0.0000000000      0.0000000000
#共Number of q in the star =    x个，这个指标出现在ph.out.坐标也在这里
     Dynamical  Matrix in cartesian axes

     q = (    0.000000000   0.000000000   0.000000000 )
#这里是两个原子指标在循环只有一个原子就是1 1
#H3S四个原子就是11 到44
#里面的矩阵是phi(icar,jcar,na,nb)jcar=1,3
实部虚部分开了，所以是六个数
    1    1
 -0.00574348  0.00000000    0.00000000  0.00000000    0.00000000  0.00000000
  0.00000000  0.00000000   -0.00574348  0.00000000    0.00000000  0.00000000
  0.00000000  0.00000000    0.00000000  0.00000000   -0.00574348  0.00000000

     Diagonalizing the dynamical matrix

     q = (    0.000000000   0.000000000   0.000000000 )

 **************************************************************************
     freq (    1) =      -1.589924 [THz] =     -53.034171 [cm-1]
     ！这里面就是极化矢量，每个原子一行，cpol(1:3, iat, imode)复数
 ( -0.022744  0.037833 -0.039393  0.065529  0.105450 -0.175412 )
 ( -0.057950 -0.013241 -0.101851 -0.001847  0.143443  0.010062 )
 ( -0.059231  0.005021 -0.101111 -0.012390  0.143443  0.010062 )
 ( -0.044662 -0.000229 -0.048211  0.057854  0.124275 -0.046520 )