<div dir="ltr">Dear, Arshad<div><br></div><div>I am facing the same problem, did you find a way to fix it?</div><div><br>Best Regards,</div><div><br></div><div>Brandon<br><br>El viernes, 19 de julio de 2019, 4:57:17 (UTC-5), Arshad Mehmood  escribió:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fonlinelibrary.wiley.com%2Fdoi%2Fabs%2F10.1002%2Fanie.201702715\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG-xkYiaATgjU3auoen6moNP1LfFQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fonlinelibrary.wiley.com%2Fdoi%2Fabs%2F10.1002%2Fanie.201702715\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG-xkYiaATgjU3auoen6moNP1LfFQ';return true;">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" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!msg/cp2k/CYGzR4T2s5U/6yS0HXWxBAAJ';return true;" onclick="this.href='https://groups.google.com/forum/#!msg/cp2k/CYGzR4T2s5U/6yS0HXWxBAAJ';return true;">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></blockquote></div></div>