User Tools

Site Tools


quantumespresso:フォノンの状態密度と分散

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:23]
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 67: Line 67:
 次のファイルを用意する。 次のファイルを用意する。
  
-<file - q2r.in>+<file - Si.q2r.in>
 &input &input
   fildyn = 'Si.dyn'   fildyn = 'Si.dyn'
Line 77: Line 77:
 出力先はSi.fcになる。 実行は次のようにする。 出力先はSi.fcになる。 実行は次のようにする。
  
-  $ q2r.x < q2r.in > q2r.out+  $ q2r.x < Si.q2r.in > Si.q2r.out
  
 以上で準備は終了である。 以上で準備は終了である。
Line 88: Line 88:
 次の入力ファイルを用意する。 次の入力ファイルを用意する。
  
-<file - matdyn.phdos.in>+<file - Si.matdyn.phdos.in>
 &input &input
   flfrc = 'Si.fc'   flfrc = 'Si.fc'
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)を指定する。 分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/a)を指定する。
  
-<file - matdyn.freq.in>+<file - Si.matdyn.freq.in>
 &input &input
   asr = 'crystal'   asr = 'crystal'
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 163: Line 163:
  
 {{:quantumespresso:ph_dispersion.png?400|}} {{:quantumespresso:ph_dispersion.png?400|}}
 +
 +  * GNUPLOTでプロットしたい場合は Si.freq.gp に出力されているものを使うと便利です。
  
  
 ==== 固有振動モードの可視化 ==== ==== 固有振動モードの可視化 ====
  
-  * [[https://www.materialscloud.org/work/tools/interactivephonon|Interactive phonon visualizer|]]を使えば、フォノンのエネルギーに対応する固有振動モードを可視化することができます。+  * [[https://www.materialscloud.org/work/tools/interactivephonon|Interactive phonon visualizer]]を使えば、フォノンのエネルギーに対応する固有振動モードを可視化することができます。
   * Si.scf.in と matdyn.modes をアップロードしてください。   * Si.scf.in と matdyn.modes をアップロードしてください。
 ===== 注意 ===== ===== 注意 =====
quantumespresso/フォノンの状態密度と分散.txt · Last modified: 2021/06/27 22:04 (external edit)