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/25 19:03] koudai [注意] |
quantumespresso:構造緩和 [2018/11/26 18:06] koudai [概要] |
||
---|---|---|---|
Line 5: | Line 5: | ||
構造緩和には次の2通りの方法があります。 | 構造緩和には次の2通りの方法があります。 | ||
- | - 格子定数はそのままで、単位胞内の原子の位置についての構造緩和 | + | - 格子定数はそのままで、単位胞内の原子の内部座標の最適化 |
- | - 単位胞内の原子位置はそのままで、格子定数についての構造緩和 | + | - 内部座標とともに格子定数も変化させる構造緩和 |
通常は2つを繰り返すことで、安定な結晶構造を探します。 | 通常は2つを繰り返すことで、安定な結晶構造を探します。 | ||
Line 47: | 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 152: | Line 152: | ||
* conv_thrが大きいと、出力ファイルに" | * conv_thrが大きいと、出力ファイルに" | ||
* 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は& | ||
- | | + | < |
- | Error in routine checkallsym (2): | + | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
- | not orthogonal operation | + | Error in routine checkallsym (2): |
- | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | + | not orthogonal operation |
+ | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
+ | </ | ||