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計算 [2019/01/26 02:00]
koudai [k点数の変更]
wien2k:scf計算 [2020/03/28 18:45]
koudai [SCF計算結果の保存]
Line 114: Line 114:
 </code> </code>
  
-init.structというファイルができます(この他にもいろいろできますが、今は無視)ので、適当なに変更しておきます。+init.structというファイルを (プロジェクト).struct に変更しておきます。
  
 <code> <code>
Line 123: Line 123:
 ===== CIFファイルからstructファイルを生成する方法 ===== ===== CIFファイルからstructファイルを生成する方法 =====
  
-実際にWIEN2kを使った第一原理計算を行うときは、対話形式でstructファイルを作成するよりもCIFファイルから生成することが多くなると思います。+実際にWIEN2kを使った第一原理計算を行うときは、対話形式でstructファイルを作成するよりも、結晶構造データの統一フォーマットであるCIFファイルから生成することが多くなると思います。
 上記のサイトからCIFファイルをダウンロードし、ファイル名をTiC.cifに変更して作業ディレクトリに保存します。 上記のサイトからCIFファイルをダウンロードし、ファイル名をTiC.cifに変更して作業ディレクトリに保存します。
  
Line 161: Line 161:
 </code> </code>
  
-ここでオプション -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。 +ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。 
-設定はオプション -h により参照できます。+ 
 +オプション意味次のとおりです 
 + 
 +^オプション^初期値^説明^ 
 +|-rkmax|7.0|RKmax| 
 +|-numk|1000|k点数| 
 + 
 + 
 +他の設定はフラグ-hにより参照できます。
  
 <code> <code>
Line 175: Line 183:
  
 <code> <code>
-$ run_lapw -ec 0.0001 -cc 0.0001 -i 100+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
 </code> </code>
  
-ここでオプション -ec はエネルギーの収束で、オプション -cc 電荷の収束です +^オプション^初期値^説明^ 
-デフォルトエネルギーの収束か見ませんが、電荷の収束見ておいたほうが無難です。 +|-cc|なし|電荷の収束。単位e。指定なければ電荷の収束はチェックされないが、見ておいたほうが無難
-また、オプション -i で最大何回までSCF計算を行うかを指定きます+|-ec|0.0001|エネルギーの収束単位はRy。デフォルトだと少し荒い| 
 +|-i|40|SCF計算を行うときの、繰り返し回数の最大値。40回だと足りないこともあるの少し多めに|
  
-他のオプションはオプション-hで確認できます。+他のオプションはフラグ-hで確認できます。
  
 <code> <code>
Line 200: Line 209:
  
 計算の収束を見るために、k点数を変更したいとします。 計算の収束を見るために、k点数を変更したいとします。
-再び init_lapw を走らせると全ての設定が初期化されますので(今は全てデフォルトなのでさほど問題ないですが)、k点数だけ変更するために x_lapw kgen を使います。+再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。
  
 <code> <code>
Line 212: Line 221:
 </code> </code>
  
-再び計算を実行するには TiC.broydX 削除しから run_lapw を実行します。+再び計算を実行するには、フラグ-NIつけて run_lapw を実行します。
  
 <code> <code>
-$ rm *.broyd* +$ run_lapw -cc 0.0001 -ec 0.00001 -NI
-$ run_lapw -ec 0.0001 -cc 0.0001+
 </code> </code>
- 
-case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 
  
 全エネルギーやフェルミ準位(半導体の場合はギャップ)が、k点数に対して収束しているか必ず確認しましょう。 全エネルギーやフェルミ準位(半導体の場合はギャップ)が、k点数に対して収束しているか必ず確認しましょう。
 +
 +
 +
 +====== SCF計算結果の保存 ======
 +
 +構造最適化前のstructファイルや計算結果のファイルを残しておきたい場合は save_lapw を使用します。
 +例えばk点を増やして計算する際に、k点数を1000で計算したときの入力ファイルや結果を k1000 という名前のディレクトリに保存しておきたい場合は次のようにします。
 +
 +<code>
 +$ save_lapw -d k1000
 +</code>
 +
 +また、この保存したファイルを復活させる場合は次のようにします
 +<code>
 +$ restore_lapw -d k1000
 +</code>
 ====== ジョブの一時停止 ====== ====== ジョブの一時停止 ======
  
Line 230: Line 252:
 </code> </code>
  
 +この空のファイルは、ジョブが止まった後に自動的に削除されます。
 ====== プログラム名について ====== ====== プログラム名について ======
  
 プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。 プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。
-例えば x_lapw は単に x だけでも問題ありません+例えば x_lapw は単に x だけで実行可能です
wien2k/scf計算.txt · Last modified: 2023/04/23 19:47 by koudai