射影バンド計算について

[fkomori]

射影バンド計算を行って、wfn_orb_proj.dataが出力されたので、以下のperlプログラムを実行しましたが、エラーメッセージが出ます。計算がうまくいっていないのでしょうか。

../bin/band-orbital-proj.pl nfenergy.data bandkpt_2D.in wfn_orb_proj.data

1.　オプションなしでは、バンド図を書くことができます。
2. 　オプション　-atom_range=1,4を入れると、以下のエラーが出ます。
Use of uninitialized value $atom in numeric lt (<) at ../bin/band-orbital-proj.pl line 218, line 147.
Use of uninitialized value $atom in numeric lt (<) at ../bin/band-orbital-proj.pl line 218, line 148.
Argument "=" isn't numeric in numeric lt (<) at ../bin/band-orbital-proj.pl line 218, line 149.
3. オプション -key=4を入れると、以下のエラーが出ます。
Use of uninitialized value $key in string ne at ../bin/band-orbital-proj.pl line 234, line 147.
（同じエラーが　.................................　　　　　　　　　　　　　　　　　　　　　 line 713.
まで出ます。）
4.　以下のオプションは正常に動作します。
-color -erange=-2,2 -einc=1 -cbrange=-1,1 -print_format=png
5.　wfn_orb_proj.dataのヘッダー部分を添付します。

小森文夫　生産技術研究所

Orbital Projection for bands.docx

[jkoga-asms]

株式会社アスムスの甲賀と申します．ご質問ありがとうございました．

ご提供いただいた情報のみでは問題を再現することができませんでした．
差し支えなければ以下のファイルをアップロードしていただけないでしょうか．

• SCF計算のnfinp.dataファイル
• 固定電荷計算(バンド計算)のnfinp.dataファイル
• bandkpt_2D.inファイル
• nfenergy.dataファイル
• （加工していない）wfn_orb_proj.dataファイル

最後の二つははぶいていただいてもよいです

ご承知の通りGitHubは特定の拡張子をもつファイルのみアップロードが可能なので，
たとえばzipにまとめていただいてからアップロードしていただければよいのではないかと思います．
よろしくお願いいたします

[fkomori]

甲賀様 ありがとうございます。後ほど入力します。 小森

…

On 2023/07/19 9:44, KOGA Junichiro wrote: 株式会社アスムスの甲賀と申します．ご質問ありがとうございました． ご提供いただいた情報のみでは問題を再現することができませんでした． 差し支えなければ以下のファイルをアップロードしていただけないでしょうか． - SCF計算のnfinp.dataファイル - 固定電荷計算(バンド計算)のnfinp.dataファイル - bandkpt_2D.inファイル - nfenergy.dataファイル - （加工していない）wfn_orb_proj.dataファイル 最後の二つははぶいていただいてもよいです ご承知の通りGitHubは特定の拡張子をもつファイルのみアップロードが可能なので， たとえばzipにまとめていただいてからアップロードしていただければよいのではないかと思います． よろしくお願いいたします

[fkomori]

甲賀様

zip ファイルを送ったつもりでしたが、見直してみるとここに入っていませんでした。以下にfileを入れてみます。
ProjBand.zip

小森

[jkoga-asms]

小森先生

ファイルをアップロードしていただきありがとうございました．

どうやら band-orbital-proj.pl はプロジェクターを全原子に割り当てた状態でバンド計算をしないと正しく動作しないようです．
先生の入力でプロジェクターが割り当てられていないのはPtだけなので，
入力ファイルの atoms ブロックにおいて，以下の要領でデフォルト設定を施しておけばよいはずです．

                atoms{
                #default weight = 1, element = Pt, mobile = 0, proj_group = 1

お手数ですがこの変更を施してやり直してみていただけないでしょうか．
バンド計算（固定電荷計算）だけでよいはずです．
よろしくお願いいたします

[fkomori]

甲賀様

早速のご返事ありがとうございます。また、ＧitＨubの使い方に慣れておらずファイル送付が遅れて申しわけありませんでした。

プロジェクターを残りのPtにも割り当てて計算すればよいとのこと了解しました。おなじPtでも割り当てていなかったPtは別分類にしたいので、それでやってみます。

小森

[jkoga-asms]

お手数ですがこの変更を施してやり直してみていただけないでしょうか． バンド計算（固定電荷計算）だけでよいはずです． よろしくお願いいたします

申し訳ありません，不正確なことを申し上げてしまいました

postprocessing{
  orbital_population{
   sw_diagonalize_population = on
  }
}

この設定を用いる場合固定電荷計算とSCF計算とでプロジェクターの設定が一致している必要があります．
sw_diagonalize_population = off とすれば固定電荷計算だけでうまくいくはずです

[fkomori]

了解しました。

sw_diagonalize_populationにも気をつけてやってみます。

小森

[junnara]

小森先生

ファイルをアップロードしていただきありがとうございました．

どうやら band-orbital-proj.pl はプロジェクターを全原子に割り当てた状態でバンド計算をしないと正しく動作しないようです． 先生の入力でプロジェクターが割り当てられていないのはPtだけなので， 入力ファイルの atoms ブロックにおいて，以下の要領でデフォルト設定を施しておけばよいはずです．

                atoms{
                #default weight = 1, element = Pt, mobile = 0, proj_group = 1

お手数ですがこの変更を施してやり直してみていただけないでしょうか． バンド計算（固定電荷計算）だけでよいはずです． よろしくお願いいたします

この点についてはマニュアルに記載しておく必要あり。