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/04/01 15:41]
koudai [atom 1 has a large sphere , consider setting HDLOs and/or larger LVNS]
wien2k:scf計算 [2020/07/07 15:33]
koudai [k点数の変更]
Line 137: Line 137:
 これで入力ファイルTiC.structができます。 これで入力ファイルTiC.structができます。
 ただし、このままだとマフィンチン半径$R_{\rm MT}$が正しく設定されていないので計算します。 ただし、このままだとマフィンチン半径$R_{\rm MT}$が正しく設定されていないので計算します。
 +拡張子の.structは不要です。
  
 <code> <code>
Line 163: Line 164:
 $ init_lapw -b -numk 1000 -rkmax 7.5 $ init_lapw -b -numk 1000 -rkmax 7.5
 </code> </code>
- +  * ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。 
-ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。 +  * 最後に init_lapw finished ok と出たら正孔です。このメッセージが出てこなかった場合、どこかで問題が起こっているので修正します。 
 +  
 オプションの意味は次のとおりです オプションの意味は次のとおりです
  
Line 212: Line 213:
  
 計算の収束を見るために、k点数を変更したいとします。 計算の収束を見るために、k点数を変更したいとします。
-再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。+再び init_lapw を走らせると全ての設定や計算結果が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。
  
 <code> <code>
Line 238: Line 239:
 最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。 最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。
 k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。 k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。
 +
 +  * run_lapw に新しいフラグをつけて再計算するときは、基本的に rm *.broyd* を使うと安全です
  
  
wien2k/scf計算.txt · Last modified: 2023/04/23 19:47 by koudai