This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
quantumespresso:構造緩和 [2018/02/22 20:33] koudai |
quantumespresso:構造緩和 [2018/03/05 01:10] koudai [原子位置の構造緩和] |
||
---|---|---|---|
Line 44: | Line 44: | ||
構造緩和計算では、calculation=' | 構造緩和計算では、calculation=' | ||
- | etot_conv_thr 構造緩和の各ステップで、エネルギーの変化がこの値よりも小さくなったときに収束したと判断する。デフォルトは1.d-4だが、経験的にはもう少し厳しいほうが良い。 | ||
- | forc_conv_thr 構造緩和の各ステップで、原子に加わる力の変化がこの値よりも小さくなったときに収束したと判断する。デフォルトは1.d-3だが、経験的にはもう少し厳しいほうが良い。 | ||
- | エネルギーのカットオフecutwfcは普段よりも多めにとっておきます。 | + | |etot_conv_thr|構造緩和の各ステップで、エネルギーの変化がこの値よりも小さくなったときに収束したと判断する。デフォルトは1.d-4だが、経験的にはもう少し小さいほうが良い。実際の計算では徐々に小さくしていって収束を見る。| |
+ | |forc_conv_thr|構造緩和の各ステップで、原子に加わる力の変化がこの値よりも小さくなったときに収束したと判断する。デフォルトは1.d-3だが、こちらも経験的にはもう少し小さいほうが良い。| | ||
また、relax計算の際は& | また、relax計算の際は& | ||
- | 原子位置について、原点にあるSiは動いてほしくないので、座標の最後に0 0 0を付け加えます。 | + | 原子位置について、原点にあるSiは動いてほしくないので、座標の最後に0 0 0を付け加えます。 |
+ | 第3成分だけ動いて欲しい場合は0 0 1とします。 | ||
+ | 何も指定していなければ1 1 1となります(すべての方向に動ける)。 | ||
実行は次のようにします。 | 実行は次のようにします。 | ||
Line 139: | Line 141: | ||
===== 注意 ===== | ===== 注意 ===== | ||
- | | + | * 実はvc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合はうまくいかないので、先にrelax計算を行ってからvc-relaxを実行しましょう。 |
- | vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。計算が正しいか確認するため、入力ファイルを得られたパラメータでもう一度vc-relax計算をすると良いでしょう。 | + | |
- | 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | + | * 結果がエネルギーのカットオフやk点数に対してちゃんと収束しているか必ず確認しましょう(得られた構造に対してエネルギーのカットオフやk点数を増やしてもう一度構造緩和してみる)。 |
+ | * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | ||
- | | + | |
- | Error in routine checkallsym (2): | + | Error in routine checkallsym (2): |
- | not orthogonal operation | + | not orthogonal operation |
- | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | + | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
Line 159: | Line 162: | ||
このときのk点は、原点を含む | このときのk点は、原点を含む | ||
K_POINTS automatic | K_POINTS automatic | ||
- | 12 12 12 1 1 1 | + | 12 12 12 0 0 0 |
などとします。 | などとします。 | ||
* http:// | * http:// | ||