Differences

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

--- quantumespresso:よくあるエラー [2019/02/13 19:24] – [NSCF計算で収束しない] koudai
+++ quantumespresso:よくあるエラー [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 4: / Line 4: @@
 入力ファイルに誤りがないことをしっかりと確認してもエラーが出る場合は、次のように対処する。
+もしそれでも解決しない場合は、がんばって解決しようとするより別の第一原理計算ソフトの使用をした方が早い。
-====== エラーの種類 ======
+====== pw.x ======
 ===== Error in routine cdiaghg (あるいはrdiaghg) =====
@@ Line 19: / Line 21: @@
   * 並列計算をしている場合、プログラム実行時に -ndiag 1 という、対角化の際に並列化を抑制するオプションをつける<code>$ mpirun -n 32 pw.x -ndiag 1 < case.scf.in > case.scf.out</code>
-  * 入力ファイルで diagonalization = 'cg' に変更する。
+  * 入力ファイルの&ELECTRONSブロックで diagonalization = 'cg' に変更する。
-  * 入力ファイルで startingwfc = 'random' に変更する。
+  * 入力ファイルの&ELECTRONSブロックで startingwfc = 'random' に変更する。
   * ウルトラソフト型の擬ポテンシャルを使っている場合、ノルム保存型に変更する（ゴーストバンドができている可能性あり）。
 以上のようにしても解決しない場合は、複数の対処法を組み合わせる。
-もしそれでも解決しない場合は、このエラーはかなり致命的であるため、別の第一原理計算ソフトの使用を推奨する。
+   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+       Error in routine cdiaghg (460):
+       S matrix not positive definite
+   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+次のときに起こりやすいエラーです
+  * 原子間距離がrelaxやvc-relaxなどの構造最適化で求まるものよりも近すぎる
+  * 伝導電子帯のバンドが多すぎる
+次のようにして解決します。
+  * relaxやvc-relaxで得られた構造を使う
+  * nbndを減らす
+  * ecutwfcとecutrhoを増やす
 ==== 参考 ====
@@ Line 49: / Line 70: @@
 対策方法
-  * ecutwfcを増やす
+  * ecutwfcとecutrhoを増やす
-  * conv_thrを大きくする（1d-8程度で十分）
+  * conv_thrをゆるくする
 ==== 参考 ====
   - http://www.democritos.it/pipermail/pw_forum/2011-September/022066.html
+====== ph.x ======
+===== DFPT with the Blochl correction is not implemented =====
+    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+       task #        18
+       from phq_readin : error #         1
+       DFPT with the Blochl correction is not implemented
+   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+テトラヘドロン法でSCF計算をしたあと、フォノン計算をしたときに発生する。
+フォノン計算で使うテトラヘドロン法のプログラムがXML形式にしか対応していないためである。
+対策方法
+プログラムをコンパイルしなおしてXMLを使えるようにする。
+  $ ./configure --enable-xml
+  $ make pw ph
+参考
+https://lists.quantum-espresso.org/pipermail/developers/2017-October/001760.html