Differences

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

--- quantumespresso:phonon:フォノンの計算 [2020/11/18 20:32]
koudai [Γ点のフォノン]
+++ quantumespresso:phonon:フォノンの計算 [2020/11/20 02:26]
koudai [X点のフォノン]
@@ Line 1: / Line 1: @@
-===== 電子状態計算 =====
+Siを例に、フォノンの簡単な計算を実行してみます。
-Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。
+===== SCF計算 =====
+事前にpw.xによるSCF計算を行います。
 <file - Si.scf.in>
@@ Line 26: / Line 28: @@
 </file>
-  * 金属の場合はsmearingを指定する必要がありますが、テトラヘドロン法はフォノン計算に使えません。
+  * 金属の場合はsmearingを指定する必要がありますが、テトラヘドロン法はフォノン計算に使えません（以下のph.x実行時にエラーが出ます）。
 <code>
@@ Line 68: / Line 70: @@
 $ ph.x < Si.phG.in > Si.phG.out
 </code>
-  * 金属の計算で |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。
 結果はfildynで指定したもの（今の場合はSi.dynG）に出力されます。
@@ Line 107: / Line 107: @@
 次のファイルを用意します。
-<file - Si.dynmat.in>
+<file - Si.dynmatG.in>
 &inputfil
   dyn ='Si.dynG',
+  filxsf = 'dynmatG.axsf'
   asr ='simple',
   lperm = .true.,
@@ Line 122: / Line 123: @@
 実行は次のようにします。
 <code>
-$ dynmat.x < Si.dynmat.in > Si.dynmat.out
+$ dynmat.x < Si.dynmatG.in > Si.dynmatG.out
 </code>
-すると dynmat.axsf というファイルができるので、以下の手順で変位ベクトルを表示します。
+すると dynmatG.axsf というファイルができるので、以下の手順で変位ベクトルを表示します。
-  - XCrysDenを起動して、[File] -> [Open Structure] -> [Open AXSF] よりdynmat.axsfを開く
+  - XCrysDenを起動して、[File] -> [Open Structure] -> [Open AXSF] よりdynmatG.axsfを開く
   - ウィンドウが開くので、矢印ボタンを押して変位ベクトルを表示したいフォノンを選ぶ（エネルギーの低い順です）
   - いったん [Hide] ボタンを押してウィンドウを最小化する
@@ Line 143: / Line 144: @@
   tr2_ph=1.0d-14,
   prefix='Si',
-  amass(1)=28.08,
   outdir='./tmp/',
   fildyn='Si.dynX',
@@ Line 183: / Line 183: @@
-<file - Si.dynmat.in>
+<file - Si.dynmatX.in>
 &inputfil
   dyn ='Si.dynX',
+  filxsf = 'dynmatX.axsf'
   asr ='no'
 /
 </file>
+<code>
+$ dynmat.x < Si.dynmatX.in > Si.dynmatX.out
+</code>

雑多な記録

User Tools

Site Tools

Differences

Page Tools