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/03/31 14:51]
koudai [k点数の変更]
wien2k:scf計算 [2020/06/24 10:56]
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 236: Line 237:
 でも可能です(電子密度の分布は前回の計算のものが引き継がれます)。 でも可能です(電子密度の分布は前回の計算のものが引き継がれます)。
  
-全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)が、k点数に対して収束しているか必ず確認しましょう。+最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。 
 +k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。
  
  
Line 269: Line 271:
 プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。 プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。
 例えば x_lapw は単に x だけで実行可能です。 例えば x_lapw は単に x だけで実行可能です。
 +
 +
 +====== エラーや警告への対処 ======
 +
 +===== atom 1 has a large sphere , consider setting HDLOs and/or larger LVNS =====
 +
 +init_lapw を実行したとき、lstartの段階でこのようなメッセージが出たときは、オプション -LVNS X (Xは4から10までの整数)をつけて、エラーが消えるまでXを増やしてください。
 +<code>
 +$ init_lapw -b -numk 1000 -rkmax 7.5 -lvns 5
 +</code>
 +
wien2k/scf計算.txt · Last modified: 2023/04/23 19:47 by koudai