<div dir="ltr">Thanks I will try it!<br>In the meanwhile I'm working on a script that reads the WFN and hopefully can combine the coefficient from different fragments.<br><br>Daniele<br><br>Il giorno martedì 2 ottobre 2018 10:18:32 UTC+2, Ole Schütt ha scritto:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Daniele,
<br>
<br>I once wrote a parser for the molecular orbital coefficients in Python:
<br>
<br>
<br><a href="https://github.com/ondrejkrejci/PPSTM/blob/master/pyPPSTM/ReadSTM.py#L520" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fondrejkrejci%2FPPSTM%2Fblob%2Fmaster%2FpyPPSTM%2FReadSTM.py%23L520\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHbfHvsho81jPvg43qMIDKczr7fDA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fondrejkrejci%2FPPSTM%2Fblob%2Fmaster%2FpyPPSTM%2FReadSTM.py%23L520\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHbfHvsho81jPvg43qMIDKczr7fDA';return true;">https://github.com/<wbr>ondrejkrejci/PPSTM/blob/<wbr>master/pyPPSTM/ReadSTM.py#L520</a>
<br>
<br>You can obtain the MO file like this:
<br>
<br> &DFT
<br> &PRINT
<br> &MO ON
<br> EIGVECS
<br> CARTESIAN
<br> FILENAME cartesian-mos
<br> &EACH
<br> QS_SCF 0
<br> &END EACH
<br> &END MO
<br> &END PRINT
<br> &END DFT
<br>
<br>
<br>-Ole
<br>
<br>
<br>On 2018-09-27 10:09, Daniele Ongari wrote:
<br>> Thank you for the hint, but I think something is going weir with the
<br>> conversion: the orbitals look fine if printed in a cube but not on
<br>> molden.
<br>> I will try to understand how the WFN is written and parse it directly
<br>> then. Thanks!
<br>>
<br>> Il giorno mercoledì 19 settembre 2018 12:21:20 UTC+2, Daniele Ongari
<br>> ha scritto:
<br>>
<br>>> Dear CP2K developers,
<br>>> I'm wondering if the binary WFN files generated by CP2K contain the
<br>>> coefficients for the gaussian basis set, so:
<br>>>
<br>>> 1) is there a tool to read the coefficients for the wave function
<br>>> from WFN files? I need it for educational purpose!
<br>>> 2) is it possible to combine the coefficients for weakly interacting
<br>>> fragments, to speed up calculations for interaction energies by
<br>>> using previously converged WFN of isolated fragments?
<br>>> So, I would need to make the WFN readable, understand the
<br>>> formatting, mix the two sets of coefficients, convert the file to a
<br>>> new binary WFN and restart the calculation from this.
<br>>> Is there already such a tool? I'm interested to work on it and
<br>>> possibly automatize the procedure using AiiDA.
<br>>>
<br>>> I found a similar old thread,
<br>>> <a href="https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A';return true;">https://groups.google.com/<wbr>forum/#!topic/cp2k/wab9r_ltk2A</a> [1], but it
<br>>> end up saying that RestartTools is not compatible with the newer
<br>>> versions of CP2K.
<br>>>
<br>>> Thanks
<br>>
<br>> --
<br>> You received this message because you are subscribed to the Google
<br>> Groups "cp2k" group.
<br>> To unsubscribe from this group and stop receiving emails from it, send
<br>> an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="GbJAsPJqAgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br>> To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="GbJAsPJqAgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>> Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br>> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>
<br>>
<br>> Links:
<br>> ------
<br>> [1] <a href="https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A';return true;">https://groups.google.com/<wbr>forum/#!topic/cp2k/wab9r_ltk2A</a>
<br>
<br></blockquote></div>