<div dir="ltr">Try using the MO section rather than MO_CUBES...it may not work.<div><br></div><div>The semi-empirical methods do not have the full functionality of the DFT methods. I'd suggest trying DFT (for very small systems, but then you know the input works), then seeing if you can do it with PM6 or similar.</div><div><br></div><div>Matt<br><br>On Monday, July 29, 2019 at 6:03:31 PM UTC+1, Richard Marchese Robinson wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">





<div lang="EN-GB" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Hi Patrick,</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Thanks for the suggestion, but this does not solve the problem.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Best regards,</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Richard</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> <a href="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a> [mailto:<a href="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a>]
<b>On Behalf Of </b>Patrick Gono<br>
<b>Sent:</b> 29 July 2019 17:47<br>
<b>To:</b> <a href="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a><br>
<b>Subject:</b> Re: [CP2K:12004] Problems computing LUMO energy for a non-periodic metallic cluster using PM6</span></p>
<p class="MsoNormal"> </p>
<div>
<div>
<p class="MsoNormal">Dear Richard,</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Try to switch off the printing of the cube files, while maintaining NHOMO and NLUMO equal to 1 in the PRINT section of DFT:</p>
</div>
<div>
<p class="MsoNormal">&DFT</p>
</div>
<div>
<p class="MsoNormal">...</p>
</div>
<div>
<p class="MsoNormal">    &PRINT<br>
        &MO_CUBES<br>
            NHOMO  1<br>
            NLUMO  1</p>
</div>
<div>
<p class="MsoNormal">            WRITE_CUBE .FALSE.</p>
</div>
<div>
<p class="MsoNormal">        &END MO_CUBES<br>
    &END PRINT<br>
&END DFT</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Does it help?</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Yours sincerely,</p>
</div>
<div>
<p class="MsoNormal">Patrick Gono</p>
</div>
</div>
<p class="MsoNormal"> </p>
<div>
<div>
<p class="MsoNormal">On Mon, 29 Jul 2019 at 18:25, Richard Marchese Robinson <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">r.l.march...@leeds.<wbr>ac.uk</a>> wrote:</p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Dear all,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Firstly, thanks to the developers for making this tool available to the community.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I have only recently started using CP2K and I am hoping to be able to compute HOMO – LUMO gaps and individual energies for non-periodic, metallic clusters using different levels
 of theory. (It is my understanding from the literature that this gap should tend to zero as the clusters become larger.)</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am trying to perfom this calculation at the PM6 level of theory.
</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I provide my input and output files.
</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">It is my understanding that the 22 eigenvalues in the output files refer to the energy levels of the occupied valence orbitals, with the Fermi energy corresponding to the energy
 of the HOMO.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">However, as far as I can tell, I cannot obtain the LUMO energy. Am I overlooking something?</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I have tried to address this using both the MO and MO_CUBES settings in the input files. I get a warning about the printing of MO cube files not being implemented for semi-empirical
 methods, e.g. PM6. Since this section is the only section which explicitly refers to the LUMO in the input file, does this mean PM6 LUMO values cannot be computed using CP2K?</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am running cp2k.7.0.190717.for.winxp.<wbr>fftw3 on Windows 7.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Thanks in advance for any assistance you have time to provide.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Best regards,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Richard</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">-----
</p>
<p class="MsoNormal">Dr Richard Marchese Robinson</p>
<p class="MsoNormal">Post-Doctoral Research Associate</p>
<p class="MsoNormal">BioRima Project [<a href="https://www.biorima.eu/" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.biorima.eu%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFgkM5aSsqf2Btt3l-fF_T_Yj6-xA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.biorima.eu%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFgkM5aSsqf2Btt3l-fF_T_Yj6-xA';return true;">https://www.biorima.eu/</a>]</p>
<p class="MsoNormal">University of Leeds</p>
<p class="MsoNormal">United Kingdom</p>
</div>
<p class="MsoNormal">-- <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="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/2c10e923-9c2c-42dc-bc1c-82662b696665%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/2c10e923-9c2c-42dc-bc1c-82662b696665%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/2c10e923-9c2c-42dc-bc1c-82662b696665%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">
https://groups.google.com/d/<wbr>msgid/cp2k/2c10e923-9c2c-42dc-<wbr>bc1c-82662b696665%<wbr>40googlegroups.com</a>.</p>
</blockquote>
</div>
<p class="MsoNormal">-- <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="javascript:" target="_blank" gdf-obfuscated-mailto="_1FIN3j1BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPL1_v-ZbDPGTiWNewdYNVkVkLGmd2DTjFboMC4vL-vMTmQq8A%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/CAPL1_v-ZbDPGTiWNewdYNVkVkLGmd2DTjFboMC4vL-vMTmQq8A%40mail.gmail.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/CAPL1_v-ZbDPGTiWNewdYNVkVkLGmd2DTjFboMC4vL-vMTmQq8A%40mail.gmail.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">
https://groups.google.com/d/<wbr>msgid/cp2k/CAPL1_v-<wbr>ZbDPGTiWNewdYNVkVkLGmd2DTjFboM<wbr>C4vL-vMTmQq8A%40mail.gmail.com</a><wbr>.</p>
</div>
</div>

</blockquote></div></div>