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/11/03 15:10] koudai [原子位置の構造緩和] |
quantumespresso:構造緩和 [2018/11/25 19:12] koudai [原子位置の構造緩和] |
||
---|---|---|---|
Line 31: | Line 31: | ||
/ | / | ||
& | & | ||
+ | | ||
/ | / | ||
&ions | &ions | ||
Line 46: | Line 47: | ||
^変数^初期値^説明^ | ^変数^初期値^説明^ | ||
- | |etot_conv_thr|1.d-4|構造緩和の各ステップで、エネルギーの変化がこの値よりも小さくなったときに収束したと判断する。実際の計算ではデフォルトの値から徐々に小さくしていって得られた原子位置の収束を見る。| | + | |etot_conv_thr|1.d-4|構造緩和の各ステップで、エネルギーの変化がこの値よりも小さくなったときに収束したと判断する。初期値のままだと粗いので小さくする。| |
- | |forc_conv_thr|1.d-3|構造緩和の各ステップで、原子に加わる力の変化がこの値よりも小さくなったときに収束したと判断する。こちらも収束を見るのを忘れないように。| | + | |forc_conv_thr|1.d-3|構造緩和の各ステップで、原子に加わる力の変化がこの値よりも小さくなったときに収束したと判断する。こちらも初期値だと粗い。| |
また、relax計算の際は& | また、relax計算の際は& | ||
Line 97: | Line 98: | ||
/ | / | ||
& | & | ||
+ | | ||
/ | / | ||
&ions | &ions | ||
Line 148: | Line 150: | ||
* vc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合は同時緩和はうまくいかないので、先にrelax計算を行った上でvc-relaxを実行しましょう。 | * vc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合は同時緩和はうまくいかないので、先にrelax計算を行った上でvc-relaxを実行しましょう。 | ||
* vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。収束するまでvc-relax計算を繰り返してください。 | * vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。収束するまでvc-relax計算を繰り返してください。 | ||
- | * 結果がエネルギーのカットオフやk点数に対してちゃんと収束しているか必ず確認しましょう(得られた構造に対してエネルギーのカットオフやk点数を増やしてもう一度構造緩和してみる)。 | + | * conv_thrが大きいと、出力ファイルに"SCF correction compared to forces is large: reduce conv_thr to get better values" |
* 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | ||
+ | < | ||
+ | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
+ | Error in routine checkallsym (2): | ||
+ | not orthogonal operation | ||
+ | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
+ | </ | ||
- | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
- | Error in routine checkallsym (2): | ||
- | not orthogonal operation | ||
- | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||