This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
quantumespresso:収束判定 [2020/07/23 19:16] koudai [概要] |
quantumespresso:収束判定 [2021/06/27 22:04] (current) |
||
---|---|---|---|
Line 37: | Line 37: | ||
ecutwfcを変えながらSCF計算を行い、total energyの収束を見る(金属の場合はFermi energyも)ことでecutwfcの値を決定する。 | ecutwfcを変えながらSCF計算を行い、total energyの収束を見る(金属の場合はFermi energyも)ことでecutwfcの値を決定する。 | ||
- | ecutwfcに対してtotal energyのグラフを作成してみるとわかりやすい。 | + | ecutwfcに対するtotal energyのグラフを作成してみるとわかりやすい。 |
- | ほとんどの場合、ecutwfcが大きくなるにしたがってtotal energyはある値に向かって漸近的に収束していくはずである。 | + | ほとんどの場合、ecutwfcが大きくなるにしたがってtotal energyはある値に向かって収束していくはずである。 |
total energyに関して欲しい有効数値の桁数が得られる最小のecutwfcを探す。 | total energyに関して欲しい有効数値の桁数が得られる最小のecutwfcを探す。 | ||
- | ただし、ecutwfcを増やせば増やすほど精度が増すかというとそうでもなく、配布されている擬ポテンシャルの有効桁数から数値計算の精度の限界は1mRyである。 | + | ただし、ecutwfcを増やせば増やすほど精度が増すかというとそうでもなく、配布されている擬ポテンシャルの有効桁数から数値計算の精度の限界は1mRy程度である。 |
ちなみに、ecutwfcが小さいと必要な状態が足りずに収束にかかるiterationの回数が大きくなることを利用して、ecutwfcを決定する方法もある。 | ちなみに、ecutwfcが小さいと必要な状態が足りずに収束にかかるiterationの回数が大きくなることを利用して、ecutwfcを決定する方法もある。 |