お世話になっております．
スピン軌道(以下SO)を考慮した上でhybrid汎関数バンド計算(cgmain, vbpef)を行うというのは現状可能でしょうか．
SOなしのhybrid汎関数計算と同じ手順でnumber_componentだけ=2とし，例えばシリコンSiについてPBE0バンド計算を行ってみますと，PBE計算のバンドからEFだけずらしたようなバンドが得られます(つまりバンドギャップおよびEfの評価が失敗しています)．
SOなしの場合と異なるインプットの設定が必要でしたらご教示いただけますと幸いです．
(仮に仕様でなく想定外の挙動だった場合，vbband.F90の129行目のif文でncomp=2の時はhybrid汎関数計算をスキップしているように見えるので，これが原因かと推察しております)