使用mokit进行GVB计算，在进行轨道定域化时出现错误

linzz

在使用mokit进行GVB计算时，出现报错。以下是输入文件和出错信息。

Enter subroutine get_paired_LMO...
chem_core=0, ecp_core=0, Nskip_UNO=0
One set of MOs: invoke RHF virtual MO projection -> localization -> paring.
$python a_rhf_proj_loc_pair.py >a_rhf_proj_loc_pair.out 2>&1

ERROR in subroutine run_command: failed to run command line
python a_rhf_proj_loc_pair.py >a_rhf_proj_loc_pair.out 2>&1

感觉是在对轨道进行定域化的时候出错的？对于LiH分子确实是nbond=1，npair=2，请问是这种情况下无法进行GVB计算，还是哪里的设置出了问题？

linzz

是因为mokit把 Li 和 Be 的core_oeb定义为0，改成1就可以算了，不知道这算不算bug？还是定义问题。具体位置在read_fch.f90