Differences

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

--- quantumespresso:収束判定 [2020/07/22 22:14] – [degauss] koudai
+++ quantumespresso:収束判定 [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 16: / Line 16: @@
 まず最初に、出発点となる値を決めなければならない。 経験的には、次の値から出発するとよい。
-  * k点数・・・単位胞内の原子数が1から5なら6*6*6、6から10なら4*4*4、10から20なら2*2*2、それ以上は1*1*1。物質の次元性に合わせて各方向にとるk点数を変えてもよい。
+  * k点数・・・単位胞内の原子数が1から3なら6*6*6、4から10なら4*4*4、10から20なら2*2*2、それ以上は1*1*1。物質の次元性に合わせて各方向にとるk点数を変えてもよい。
   * ecutwfc・・・まずは20, 30, 40, 60 Ryで計算して様子を見る。ウルトラソフト型の擬ポテンシャルを使う場合はその半分。
   * ecutrho・・・ecutwfcの4倍。ただし擬ポテンシャルがウルトラソフト型の場合はそれ以上（およそ8~12倍が必要と言われているが、とりあえず最初は12倍しておく）。
@@ Line 37: / Line 37: @@
 ecutwfcを変えながらSCF計算を行い、total energyの収束を見る（金属の場合はFermi energyも）ことでecutwfcの値を決定する。
-ecutwfcに対してtotal energyのグラフを作成してみるとわかりやすい。
+ecutwfcに対するtotal energyのグラフを作成してみるとわかりやすい。
-ほとんどの場合、ecutwfcが大きくなるにしたがってtotal energyはある値に向かって漸近的に収束していくはずである。
+ほとんどの場合、ecutwfcが大きくなるにしたがってtotal energyはある値に向かって収束していくはずである。
 total energyに関して欲しい有効数値の桁数が得られる最小のecutwfcを探す。
-ただし、ecutwfcを増やせば増やすほど精度が増すかというとそうでもなく、配布されている擬ポテンシャルの有効桁数から数値計算の精度の限界は1mRyである。
+ただし、ecutwfcを増やせば増やすほど精度が増すかというとそうでもなく、配布されている擬ポテンシャルの有効桁数から数値計算の精度の限界は1mRy程度である。
 ちなみに、ecutwfcが小さいと必要な状態が足りずに収束にかかるiterationの回数が大きくなることを利用して、ecutwfcを決定する方法もある。