User Tools

Site Tools


wien2k:scf計算

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
wien2k:scf計算 [2020/07/07 15:52]
koudai
wien2k:scf計算 [2021/06/27 22:04]
127.0.0.1 external edit
Line 165: Line 165:
 </code> </code>
   * ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。   * ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
-  * 最後に init_lapw finished ok と出たら正孔です。このメッセージが出てこなかった場合、どこかで問題が起こっているので修正します。+  * 最後に init_lapw finished ok と出たらうまく行っています。このメッセージが出てこなかった場合、どこかで問題が起こっているので修正します。
      
 オプションの意味は次のとおりです オプションの意味は次のとおりです
Line 202: Line 202:
  
 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。
- +<code> 
-  $ grep ":ENE" TiC.scf     # 全エネルギー +$ grep ":ENE" TiC.scf     # 全エネルギー 
-  $ grep ":FER" TiC.scf     # フェルミ準位。単位は Ry です。 +(略) 
-  $ grep ":GAP" TiC.scf     # エネルギーギャップ。金属の場合は0になります。+:ENE  : ********** TOTAL ENERGY IN Ry =        -1783.96056300 
 +$ grep ":FER" TiC.scf     # フェルミ準位。単位は Ry です。 
 +(略) 
 +:FER  : F E R M I - ENERGY(TETRAH.M.)=   0.7492912035 
 +$ grep ":GAP" TiC.scf     # エネルギーギャップ。金属の場合は0になります。 
 +(略) 
 +:GAP (global)   :    0.0    Ry =     0.0   eV (metal) 
 +</code>
  
 これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。
Line 218: Line 225:
 </code> </code>
  
-もしrun_lapwのフラグに変更がある場合は、ブロイデン法を使ってSCF計算を収束させているので、その履歴を削除してから再実行します。+もしrun_lapwのフラグに変更がある場合は、ブロイデン法の履歴を削除してから再実行します。
 (電子密度の分布は前回の計算のものが引き継がれます) (電子密度の分布は前回の計算のものが引き継がれます)
  
Line 255: Line 262:
  
  
-====== SCF計算結果の保存 ======+====== 計算結果の保存 ======
  
 計算の入出力ファイルを残しておきたい場合は save_lapw を使用します。 計算の入出力ファイルを残しておきたい場合は save_lapw を使用します。
wien2k/scf計算.txt · Last modified: 2023/04/23 19:47 by koudai