Differences

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

--- quantumespresso:phonon:フォノンの計算 [2020/11/18 19:48]
koudai [Γ点のフォノン]
+++ quantumespresso:phonon:フォノンの計算 [2020/11/18 20:08]
koudai [X点のフォノン]
@@ Line 7: / Line 7: @@
 &control
    calculation='scf',
-   prefix='si'
+   prefix='Si'
-   pseudo_dir = './',
+   pseudo_dir = './pseudo/',
-   outdir='./'
+   outdir='./tmp/'
 /
 &system
@@ Line 39: / Line 39: @@
   prefix='Si',
   epsil=.true.,
-  outdir='./',
+  outdir='./tmp/',
-  fildyn='si.dynG',
+  fildyn='Si.dynG',
  /
 .0 0.0 0.0
@@ Line 57: / Line 57: @@
 |fildyn|matdyn|結果であるdynamical matrixを出力するファイル名。|
-金属の場合は計算するk点を 0.01 0.0 0.0 のように少しずらしてやるとうまくいきます。
+金属の場合は、epsil=.false.として、計算するk点を 0.01 0.0 0.0 のようにΓ点から少しずらしてやるとうまくいきます。
 実行は次のようにします。
   $ ph.x < Si.phG.in > Si.phG.out
-結果はsi.dynGで見ることができる。
+結果はfildynで指定したもの（今の場合はSi.dynG）に出力されます。
-<file - si.dynG>
+<file - Si.dynG>
 （中略）
  **************************************************************************
@@ Line 94: / Line 94: @@
   * 振動モードの既約表現はSi.phG.out内に出力されている。
+フォノンの振動に対応する変位ベクトルはXcrysDenを使って見ることができます。
+次のファイルを用意します。
+<file - Si.dynmat.in>
+&inputfil
+  dyn ='Si.dynG',
+  asr ='simple',
+  lperm = .true.,
+  q(1)=1.0,
+  q(2)=0.0,
+  q(3)=0.0
+/
+</file>
+Γ点のフォノンは q->0 の極限で定義されます。lpermは q=(q(1), q(2), q(3)) に沿って q->0 の極限をとることを指定します。
+実行は次のようにします。
+<code>
+$ dynmat.x < Si.dynmat.in > Si.dynmat.out
+</code>
+すると dynmat.axsf というファイルができるので、以下の手順で変位ベクトルを表示します。
+  - XCrysDenを起動して、[File] -> [Open Structure] -> [Open AXSF] よりdynmat.axsfを開く
+  - ウィンドウが開くので、矢印ボタンを押して変位ベクトルを表示したいフォノンを選ぶ（エネルギーの低い順です）
+  - いったん [Hide] ボタンを押してウィンドウを最小化する
+  - [Display] から [Forces] を選ぶ
+  - このままだと矢印が長すぎるので、[Modify] -> [Force Settings] でLength Factorの大きさを調節する（デフォルトは200だが30くらいがちょうどよい）
 ===== X点のフォノン =====
@@ Line 140: / Line 169: @@
 重に縮退したものが3つあることがわかる。
+変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。
+<file - Si.dynmat.in>
+&inputfil
+  dyn ='Si.dynX',
+  asr ='no'
+/
+</file>

雑多な記録

User Tools

Site Tools

Differences

Page Tools