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:フォノンの状態密度と分散 [2019/01/07 21:02]
koudai [分散]
quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:49]
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になる。 実行は次のようにする。
  
-  $ mpirun 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:
 実行は次のようにする。 実行は次のようにする。
  
-  $ mpirun matdyn.x < matdyn.phdos.in > matdyn.phdos.out+  $ matdyn.x < Si.matdyn.phdos.in > Si.matdyn.phdos.out 
 + 
 +下図は Si.phdos に出力されたデータをプロットした。
  
 {{:quantumespresso:phdos.png?400|}} {{:quantumespresso:phdos.png?400|}}
Line 118: Line 120:
 ===== 分散 ===== ===== 分散 =====
  
-分散の場合は計算する波数点の経路を指定する。+分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/a)を指定する。
  
-<file - matdyn.freq.in>+<file - Si.matdyn.freq.in>
 &input &input
   asr = 'crystal'   asr = 'crystal'
Line 139: Line 141:
 実行はさきほどと同様。 実行はさきほどと同様。
  
-  $ mpirun matdyn.x < matdyn.freq.in > matdyn.freq.out+  $ matdyn.x < Si.matdyn.freq.in > Si.matdyn.freq.out
  
 結果はSi.freqに出力される。 結果はSi.freqに出力される。
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:
 {{:quantumespresso:ph_dispersion.png?400|}} {{:quantumespresso:ph_dispersion.png?400|}}
  
 +  * GNUPLOTでプロットしたい場合は Si.freq.gp に出力されているものを使うと便利です。
 +
 +
 +==== 固有振動モードの可視化 ====
  
 +  * [[https://www.materialscloud.org/work/tools/interactivephonon|Interactive phonon visualizer]]を使えば、フォノンの各エネルギーに対応する固有振動モードを可視化することができます。
 +  * Si.scf.in と matdyn.modes をアップロードしてください。matdyn.modesはmatdyn.xを実行すると自動的に生成されます。
 ===== 注意 ===== ===== 注意 =====
  
quantumespresso/フォノンの状態密度と分散.txt · Last modified: 2021/06/27 22:04 (external edit)