Differences

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

--- quantumespresso:構造緩和 [2018/02/22 20:36]
koudai [金属の場合]
+++ quantumespresso:構造緩和 [2018/03/05 01:20]
koudai [格子定数の構造緩和]
@@ Line 44: / Line 44: @@
 構造緩和計算では、calculation='relax'を指定します。 関連するオプションの意味は次の通りです。
-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計算の際は&ionsという項目が必要です。 ここに必要に応じて構造緩和計算のオプションを書き加えますが、今はデフォルトのままとします。
-原子位置について、原点にあるSiは動いてほしくないので、座標の最後に0 0 0を付け加えます。 他にも動いてほしくない原子があれば、同様に書き加えます。
+原子位置について、原点にあるSiは動いてほしくないので、座標の最後に0 0 0を付け加えます。
+第3成分だけ動いて欲しい場合は0 0 1とします。
+何も指定していなければ1 1 1となります（すべての方向に動ける）。
 実行は次のようにします。
@@ Line 135: / Line 137: @@
 </file>
-もともと基本並進ベクトルが格子定数を単位としてa(1)=(-0.5, 0, 0.5)だったのが(-0.510132378, 0.000000000, 0.510132378)に拡大しています。 そのため、構造緩和によって得られた格子定数はボーア半径を単位として10.0*0.510132378/0.5=10.20・・・となります。
+もともと基本並進ベクトルが格子定数を単位としてa(1)=(-0.5, 0, 0.5)だったのが(-0.510132378, 0.000000000, 0.510132378)に拡大しています。そのため、構造緩和によって得られた格子定数はボーア半径を単位として10.0 * 0.510132378 / 0.5 = 10.20... となります。
 ===== 注意 =====
@@ Line 141: / Line 143: @@
   * 実はvc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合はうまくいかないので、先にrelax計算を行ってからvc-relaxを実行しましょう。
   * vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。計算が正しいか確認するため、入力ファイルを得られたパラメータでもう一度vc-relax計算をすると良いでしょう。
+  * 結果がエネルギーのカットオフやk点数に対してちゃんと収束しているか必ず確認しましょう（得られた構造に対してエネルギーのカットオフやk点数を増やしてもう一度構造緩和してみる）。
   * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は&systemに対称性を考慮しないというオプションnosym=.true.を加えてください。

雑多な記録

User Tools

Site Tools

Differences

Page Tools