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/24 11:21]
koudai
Line 148: 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 163: Line 163:
 ====== 金属の場合 ====== ====== 金属の場合 ======
  
-金属の場合は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
  
 +
 +
 +====== 擬ポテンシャルについて ======
 +
 +一般に構造最適化で得られた格子定数はPZなどのLDA計算では過小評価され、PBEなどのGGA計算では過大評価されます。
 +最近ではPBEを改善したPBESOLという交換相関ポテンシャルが、実験の格子定数をよく再現するということで使われることが多いです。
quantumespresso/構造緩和.txt · Last modified: 2024/03/30 23:58 by koudai