Hi Daniele,<br /><br />I'm encountering a similar issue, and I'm interested in reading a binary .wfn file, making some changes, and then restarting the calculation with a newly created .wfn file. I was wondering if you have found a way or have a script for this task. Any insights on creating a new .wfn for restarting calculations would be greatly appreciated.<br /><br />Best regards,<br />Lili<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, October 4, 2018 at 9:11:39 AM UTC-4 Daniele Ongari wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><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</div><div dir="ltr"><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" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://github.com/ondrejkrejci/PPSTM/blob/master/pyPPSTM/ReadSTM.py%23L520&source=gmail&ust=1695220295326000&usg=AOvVaw1_g4aiz2O9wbAr-44wS_Fw">https://github.com/ondrejkrejci/PPSTM/blob/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" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/forum/%23!topic/cp2k/wab9r_ltk2A&source=gmail&ust=1695220295326000&usg=AOvVaw1kbUuFpityCCcxHDQuyT7U">https://groups.google.com/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 rel="nofollow">cp2k+uns...@googlegroups.com</a>.
<br>> To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.
<br>> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/group/cp2k&source=gmail&ust=1695220295326000&usg=AOvVaw3NvrzjsMEA-JLSKADXKCTh">https://groups.google.com/group/cp2k</a>.
<br>> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/optout&source=gmail&ust=1695220295326000&usg=AOvVaw2O8nMYmVfSNYRNS9zm8vpl">https://groups.google.com/d/optout</a>.
<br>> 
<br>> 
<br>> Links:
<br>> ------
<br>> [1] <a href="https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/forum/%23!topic/cp2k/wab9r_ltk2A&source=gmail&ust=1695220295326000&usg=AOvVaw1kbUuFpityCCcxHDQuyT7U">https://groups.google.com/forum/#!topic/cp2k/wab9r_ltk2A</a>
<br>
<br></blockquote></div></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/2e7297aa-54ab-4205-b649-2a6670af8a43n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/2e7297aa-54ab-4205-b649-2a6670af8a43n%40googlegroups.com</a>.<br />