Differences

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

--- quantumespresso:phonon:フォノンの計算 [2020/11/18 20:06]
koudai [電子状態計算]
+++ quantumespresso:phonon:フォノンの計算 [2020/11/18 20:12]
koudai [X点のフォノン]
@@ Line 39: / Line 39: @@
   prefix='Si',
   epsil=.true.,
-  outdir='./',
+  outdir='./tmp/',
   fildyn='Si.dynG',
  /
@@ Line 128: / Line 128: @@
 今度はX点(2π/a,0,0)のフォノンを調べる。 手順はさきほどと同じであるが、誘電率が計算できないことに注意する。
-<file - si.phX.in>
+<file - Si.phX.in>
 phonons of Si at X
  &inputph
   tr2_ph=1.0d-14,
-  prefix='si',
+  prefix='Si',
   amass(1)=28.08,
-  outdir='./',
+  outdir='./tmp/',
-  fildyn='si.dynX',
+  fildyn='Si.dynX',
  /
 .0 0.0 0.0
 </file>
-  $ ph.x < si.phX.in > si.phX.out
+  $ ph.x < Si.phX.in > Si.phX.out
-結果はsi.dynXに出力される。
+結果は Si.dynX に出力される。
-<file -si.dynX>
+<file -Si.dynX>
 （中略）
  **************************************************************************
@@ Line 169: / Line 169: @@
 重に縮退したものが3つあることがわかる。
+変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。
+<file - Si.dynmat.in>
+&inputfil
+  dyn ='Si.dynX',
+  asr ='no'
+/
+</file>

雑多な記録