<div dir="ltr">Dear Developers,<div>I need molecular orbitals in MOLDEN format for post-processing (not just visualization)  to calculate <a href="https://onlinelibrary.wiley.com/doi/abs/10.1002/anie.201702715">Orbital Overlap Distance Function</a>. </div><div>I have modified the lines in "molden_utlis" as recommended in the thread "<a href="https://groups.google.com/forum/#!msg/cp2k/CYGzR4T2s5U/6yS0HXWxBAAJ">MOS in Molden Format</a>". Namely;</div><div><span style="color: rgb(0, 0, 0);"><br></span></div><div><span style="color: rgb(0, 0, 0);">WRITE (iw, fmtstr1) irow_in+orbital - 1 , mo_coeff(orbmap(orbital))</span><br></div><div><span style="color: rgb(0, 0, 0);"><br></span></div><div><span style="color: rgb(0, 0, 0);">I am getting .MOLDEN output file containing orbitals in correct format.</span></div><div><font color="#000000">However, the molecular orbitals fail the normalization check. The sum of MOs occupancies are largely different from the analytically integrated number of electrons. For small molecules e.g. water at  PADE DZVP-GTH-PBE/GTH-PBE the difference is 0.0021, but for system like 6x6x1 graphene super-cell, the difference reaches to 10.348.</font></div><div><font color="#000000">I shall be very much thankful if you suggest me the reason and the possible fix. I need the correct number of electrons for my calculations.</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Best Regards,</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Arshad Mehmood</font></div><div><font color="#000000">Texas Christian University.   </font><span style="color: rgb(0, 0, 0);"> </span></div></div>