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:フォノンの状態密度と分散 [2020/11/19 02:41] koudai [実空間へのフーリエ変換] |
quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:48] koudai [固有振動モードの可視化] |
||
---|---|---|---|
Line 60: | Line 60: | ||
$ ph.x < Si.ph.in > Si.ph.out | $ ph.x < Si.ph.in > Si.ph.out | ||
- | フォノンの計算には時間がかかるので注意。 | + | フォノンの計算には時間がかかるので、メッシュのとりすぎに注意。 |
===== 実空間へのフーリエ変換 ===== | ===== 実空間へのフーリエ変換 ===== | ||
Line 88: | Line 88: | ||
次の入力ファイルを用意する。 | 次の入力ファイルを用意する。 | ||
- | <file - matdyn.phdos.in> | + | <file - Si.matdyn.phdos.in> |
&input | &input | ||
flfrc = ' | flfrc = ' | ||
Line 111: | Line 111: | ||
実行は次のようにする。 | 実行は次のようにする。 | ||
- | $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out | + | $ matdyn.x < Si.matdyn.phdos.in > Si.matdyn.phdos.out |
下図は Si.phdos に出力されたデータをプロットした。 | 下図は Si.phdos に出力されたデータをプロットした。 | ||
Line 122: | Line 122: | ||
分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/ | 分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/ | ||
- | <file - matdyn.freq.in> | + | <file - Si.matdyn.freq.in> |
&input | &input | ||
asr = ' | asr = ' | ||
Line 141: | Line 141: | ||
実行はさきほどと同様。 | 実行はさきほどと同様。 | ||
- | $ matdyn.x < matdyn.freq.in > matdyn.freq.out | + | $ matdyn.x < Si.matdyn.freq.in > Si.matdyn.freq.out |
結果はSi.freqに出力される。 | 結果はSi.freqに出力される。 | ||
Line 164: | Line 164: | ||
{{: | {{: | ||
- | * GNUPLOTの場合は Si.freq.gp に出力されているものを使うと便利です。 | + | * GNUPLOTでプロットしたい場合は Si.freq.gp に出力されているものを使うと便利です。 |
==== 固有振動モードの可視化 ==== | ==== 固有振動モードの可視化 ==== | ||
- | * [[https:// | + | * [[https:// |
* Si.scf.in と matdyn.modes をアップロードしてください。 | * Si.scf.in と matdyn.modes をアップロードしてください。 | ||
===== 注意 ===== | ===== 注意 ===== |