第一原理計算ソフトを使っていて、ほぼ確実に出会うのが計算途中のエラーである。 入力ファイルに誤りがないことをしっかりと確認してもエラーが出る場合は、次のように対処する。

もしそれでも解決しない場合は、がんばって解決しようとするより別の第一原理計算ソフトの使用をした方が早い。

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine  cdiaghg (218):
      problems computing cholesky
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

行列の対角化のさいに、何らかの特異点が生じているためにエラーが発生する。 バンドが多いときに起こりやすい。 以下のことを順番に行い、エラーが消えるか確認する。

• 並列計算をしている場合、プログラム実行時に -ndiag 1 という、対角化の際に並列化を抑制するオプションをつける

  $ mpirun -n 32 pw.x -ndiag 1 < case.scf.in > case.scf.out

• 入力ファイルの&ELECTRONSブロックで diagonalization = 'cg' に変更する。
• 入力ファイルの&ELECTRONSブロックで startingwfc = 'random' に変更する。
• ウルトラソフト型の擬ポテンシャルを使っている場合、ノルム保存型に変更する（ゴーストバンドができている可能性あり）。

以上のようにしても解決しない場合は、複数の対処法を組み合わせる。

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine cdiaghg (460):
     S matrix not positive definite
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

次のときに起こりやすいエラーです

• 原子間距離がrelaxやvc-relaxなどの構造最適化で求まるものよりも近すぎる
• 伝導電子帯のバンドが多すぎる

次のようにして解決します。

• relaxやvc-relaxで得られた構造を使う
• nbndを減らす
• ecutwfcとecutrhoを増やす

1. http://www.quantum-espresso.org/faq/frequent-errors-during-execution/
2. http://www.berkeleygw.org/releases/manual_v1.0.6.html#MeanField/ESPRESSO/README

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine c_bands (1):
     too many bands are not converged
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

NSCF計算でバンドが多いときに発生しやすい。

対策方法

• ecutwfcとecutrhoを増やす
• conv_thrをゆるくする

1. http://www.democritos.it/pipermail/pw_forum/2011-September/022066.html

  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     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