Differences

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

--- quantumespresso:フォノンの状態密度と分散 [2019/01/07 19:53]
koudai [dynamical matrixの計算]
+++ quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:40]
koudai [分散]
@@ Line 40: / Line 40: @@
 <file - Si.ph.in>
 phonons of Si
- &inputph
+&inputph
-  tr2_ph=1.0d-14,
+  tr2_ph = 1.0d-14,
-  prefix='Si',
+  prefix = 'Si',
-  outdir='./tmp/',
+  outdir = './tmp/',
-  fildyn='Si.dyn',
+  fildyn = 'Si.dyn',
-  ldisp=.true.,
+  ldisp  = .true.,
-  nq1=4,
+  nq1    = 4,
-  nq2=4,
+  nq2    = 4,
-  nq3=4,
+  nq3    = 4,
- /
+/
 </file>
@@ Line 58: / Line 58: @@
 実行は次のようにする。
-  $ ph.x < si.ph.in > si.ph.out
+  $ ph.x < Si.ph.in > Si.ph.out
 フォノンの計算には時間がかかるので注意。
@@ Line 90: / Line 90: @@
 <file - matdyn.phdos.in>
 &input
-  flfrc = 'si.fc'
+  flfrc = 'Si.fc'
-  fldos = 'si.phdos'
+  fldos = 'Si.phdos'
   asr = 'crystal'
   dos = .true.
@@ Line 112: / Line 112: @@
   $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out
+下図は Si.phdos に出力されたデータをプロットした。
 {{:quantumespresso:phdos.png?400|}}
@@ Line 118: / Line 120: @@
 ===== 分散 =====
-分散の場合は計算する波数点の経路を指定する。
+分散の場合は計算する波数点の経路（デカルト座標で、単位は格子定数aを使って2π/a）を指定する。
 <file - matdyn.freq.in>
 &input
   asr = 'crystal'
-  flfrc = 'si.fc'
+  flfrc = 'Si.fc'
-  flfrq = 'si.freq'
+  flfrq = 'Si.freq'
   q_in_band_form=.true.
 /
@@ Line 145: / Line 147: @@
   $ plotband.x
-  Input file > si.freq
+  Input file > Si.freq
   Reading 6 bands at 41 k-points
   Range: 0.0000 510.2482eV Emin, Emax > 0, 550
@@ Line 162: / Line 164: @@
 {{:quantumespresso:ph_dispersion.png?400|}}
+  * GNUPLOTの場合は Si.freq.gp に出力されているものを使うと便利です。
+==== 固有振動モードの可視化 ====
+  * [[https://www.materialscloud.org/work/tools/interactivephonon|Interactive phonon visualizer]]を使えば、各フォノンのエネルギーに対応する固有振動モードを可視化することができます。
+  * Si.scf.in と matdyn.modes をアップロードしてください。
 ===== 注意 =====

雑多な記録

User Tools

Site Tools

Differences

Page Tools