This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
quantumespresso:フォノンの状態密度と分散 [2019/01/07 21:01] koudai [実空間へのフーリエ変換] |
quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:41] koudai [実空間へのフーリエ変換] |
||
---|---|---|---|
Line 67: | Line 67: | ||
次のファイルを用意する。 | 次のファイルを用意する。 | ||
- | <file - q2r.in> | + | <file - Si.q2r.in> |
&input | &input | ||
fildyn = ' | fildyn = ' | ||
Line 77: | Line 77: | ||
出力先はSi.fcになる。 実行は次のようにする。 | 出力先はSi.fcになる。 実行は次のようにする。 | ||
- | $ mpirun | + | $ q2r.x < Si.q2r.in > Si.q2r.out |
以上で準備は終了である。 | 以上で準備は終了である。 | ||
Line 112: | Line 112: | ||
$ matdyn.x < matdyn.phdos.in > matdyn.phdos.out | $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out | ||
+ | |||
+ | 下図は Si.phdos に出力されたデータをプロットした。 | ||
{{: | {{: | ||
Line 118: | Line 120: | ||
===== 分散 ===== | ===== 分散 ===== | ||
- | 分散の場合は計算する波数点の経路を指定する。 | + | 分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/ |
<file - matdyn.freq.in> | <file - matdyn.freq.in> | ||
Line 145: | Line 147: | ||
$ plotband.x | $ plotband.x | ||
- | Input file > si.freq | + | Input file > Si.freq |
Reading 6 bands at 41 k-points | Reading 6 bands at 41 k-points | ||
Range: 0.0000 510.2482eV Emin, Emax > 0, 550 | Range: 0.0000 510.2482eV Emin, Emax > 0, 550 | ||
Line 162: | Line 164: | ||
{{: | {{: | ||
+ | * GNUPLOTの場合は Si.freq.gp に出力されているものを使うと便利です。 | ||
+ | |||
+ | |||
+ | ==== 固有振動モードの可視化 ==== | ||
+ | * [[https:// | ||
+ | * Si.scf.in と matdyn.modes をアップロードしてください。 | ||
===== 注意 ===== | ===== 注意 ===== | ||