Differences

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

--- wien2k:scf計算 [2017/10/28 00:16]
koudai
+++ wien2k:scf計算 [2017/10/28 01:08]
koudai
@@ Line 3: / Line 3: @@
   * [[https://ja.wikipedia.org/wiki/%E7%82%AD%E5%8C%96%E3%83%81%E3%82%BF%E3%83%B3|炭化チタン - Wikipedia]]
-====== 準備 ======
+第一原理計算の基本的な流れは
+  - 結晶構造データの用意
+  - 計算方法の設定
+  - 実際の電子状態計算（self-consistent field計算、SCF計算）
+  - 各種物理量の計算
+となります。
+ここでは3番目までの方法を解説します。
+====== 結晶構造データの作成 ======
 WIEN2kでは結晶構造はcase.structというファイルに保存されます。
@@ Line 16: / Line 27: @@
 </code>
-structファイルの作り方は、プログラムmakestructを使って自分で結晶構造を入力する方法と、結晶構造の統一フォーマットであるCIFファイルを変換する方法の2通りがあります。
+structファイルの作り方は、プログラムmakestruct_lapwを使って自分で結晶構造を入力する方法と、結晶構造の統一フォーマットであるCIFファイルを変換する方法の2通りがあります。
 既知の物質を調べる場合、結晶構造は論文か物質科学系のデータベースサイトを参照すると良いでしょう。
@@ Line 24: / Line 35: @@
-===== makestructを使う方法 =====
+===== makestruct_lapwを使う方法 =====
 基本的には質問に答えていくだけで完了します。
@@ Line 30: / Line 41: @@
 <code>
-$ makestruct
+$ makestruct_lapw
 ************************************************
@@ Line 125: / Line 136: @@
 $ mv 5910091.struct TiC.struct
 </code>
+====== 計算の設定の初期化 ======
+SCF計算に入る前に、計算の設定の初期化をおこないます。
+ほとんどの場合はデフォルトのままでよいですが、ここではk点数だけを設定しましょう。
+k点数を増やせば計算の精度はあがりますが、計算時間もかかるようになります。
+最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。
+<code>
+$ init_lapw -b -numk 1000
+</code>
+ここでオプション-bはバッチモードで、これをつけることで指定のもの以外はすべてデフォルトの値に設定されます。
+なお、他の設定はオプション-hにより参照できます。
+<code>
+$ init -h
+</code>
+====== SCF計算の実行 ======
+これで準備が整いましたので、計算を実行します。
+プログラムrun_lapwを使用します。
+<code>
+$ run_lapw -ec 0.0001 -cc 0.0001
+</code>
+ここでオプション-ecはエネルギーの収束で、オプション-ccは電荷の収束です。
+デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。
+実行後、いろいろファイルができますが、重要なのはTiC.scfです。
+今回は金属なのでフェルミエネルギーが出力されています。
+半導体であればギャップサイズが出力されます。
+これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。

雑多な記録

User Tools

Site Tools

Differences

Page Tools