Calling fortran routine from python: In the case of orbital conversion from PySCF to GAMESS
Unfortunately the definition of MO(Molecular Orbitals) depends on the particular program. Therefore , the MO format of PySCF is different from that of GAMESS. For example, in the case of Carbon atom with 6-31G* basis set, we need to rearrange the sequence of AO for each MO coefficients as # For the Second Row atoms with the 6-31G* basis set case # # We should transform the Second Row Elemenents as # # PySCF --> GAMESS # 1s --> 1s # 2s --> 2s # 3s --> 2px # 2px --> 2py # 2py --> 2pz # 2pz --> 3s # 3px --> 3px # 3py --> 3py # 3pz --> 3pz # 3dxx --> 3dxx # 3dxy --> 3dyy # 3dxz --> 3dzz # 3dyy --> 3dxy # 3dyz --> 3dxz # 3dzz --> 3dyz # This post will show you how to 1) perform PySCF calculation to obtain MO 2) transform the MO into GAMESS format 3) write the GAMESS formatted MO into $VEC format of GAMESS