User Tools

Site Tools


quantumespresso:構造緩和

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
quantumespresso:構造緩和 [2018/11/26 18:06]
koudai [概要]
quantumespresso:構造緩和 [2020/03/18 16:28]
koudai [金属の場合]
Line 2: Line 2:
  
 構造緩和計算により、結晶の安定な構造を求めます。 構造緩和計算により、結晶の安定な構造を求めます。
- 
 構造緩和には次の2通りの方法があります。 構造緩和には次の2通りの方法があります。
  
Line 149: Line 148:
  
   * vc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合は同時緩和はうまくいかないので、先にrelax計算を行った上でvc-relaxを実行しましょう。   * vc-relax計算は原子位置の緩和も同時に行っています。ほとんどの場合は同時緩和はうまくいかないので、先にrelax計算を行った上でvc-relaxを実行しましょう。
-  * vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。収束するまでvc-relax計算を繰り返してください。+  * vc-relax計算で構造緩和を行うとき、計算に必要なパラメータは最初の格子定数と原子位置の計算から出したものを使い続けます。結果が変わらなくなるまでvc-relax計算を繰り返してください。
   * conv_thrが大きいと、出力ファイルに"SCF correction compared to forces is large: reduce conv_thr to get better values"と注意が出てくるので小さくしましょう。   * conv_thrが大きいと、出力ファイルに"SCF correction compared to forces is large: reduce conv_thr to get better values"と注意が出てくるので小さくしましょう。
   * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は&systemに対称性を考慮しないというオプションnosym=.true.を加えてください。   * 構造緩和の結果、結晶が別の対称性になってしまう場合、プログラムが以下のエラーを出して止まってしまいます。その場合は&systemに対称性を考慮しないというオプションnosym=.true.を加えてください。
Line 164: Line 163:
 ====== 金属の場合 ====== ====== 金属の場合 ======
  
-金属の場合はFermi面の情報が必要ですので、+金属の場合はFermi面の情報が必要ですので、&systemに
  
   occupations='tetrahedra_opt'   occupations='tetrahedra_opt'
  
-使用してください。 +追加してください。
-このときのk点は、原点を含む +
-  K_POINTS automatic +
-    12 12 12 0 0 0 +
-などとします+
  
 * http://qe-forge.org/pipermail/pw_forum/2017-October/114051.html * http://qe-forge.org/pipermail/pw_forum/2017-October/114051.html
  
quantumespresso/構造緩和.txt · Last modified: 2024/03/30 23:58 by koudai