<div>I ran some others these are the results:</div><div> </div><div><table width="590" style="width: 443pt; border-collapse: collapse;" border="0" cellSpacing="0" cellPadding="0">
<colgroup><col width="185" style="width: 139pt; mso-width-source: userset; mso-width-alt: 6765;">
<col width="135" style="width: 101pt; mso-width-source: userset; mso-width-alt: 4937;">
<col width="64" style="width: 48pt;" span="3">
<col width="78" style="width: 59pt; mso-width-source: userset; mso-width-alt: 2852;">
<tbody><tr height="20" style="height: 15pt;">
<td width="185" height="20" style="border: 0px black; width: 139pt; height: 15pt; background-color: transparent;"><font face="Calibri">Key Words</font></td>
<td width="135" style="border: 0px black; width: 101pt; background-color: transparent;"></td>
<td width="64" style="border: 0px black; width: 48pt; background-color: transparent;"><font face="Calibri">%HF</font></td>
<td width="64" style="border: 0px black; width: 48pt; background-color: transparent;"><font face="Calibri">CP2K</font></td>
<td width="64" style="border: 0px black; width: 48pt; background-color: transparent;"><font face="Calibri">G09</font></td>
<td width="78" style="border: 0px black; width: 59pt; background-color: transparent;"><font face="Calibri">diff (kJ/mol)</font></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_GGA_XC_HCTH_407</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.40895</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.40918</font></td>
<td align="right" class="xl63" style="border: 0px black; background-color: transparent;"><font face="Calibri">-0.6</font></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_HYB_GGA_XC_B97_1</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0.21</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.39399</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.39374</font></td>
<td align="right" class="xl63" style="border: 0px black; background-color: transparent;"><font face="Calibri">0.7</font></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_MGGA_C_M06_L</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-67.68914</font></td>
<td class="xl64" style="border: 0px black; background-color: transparent;"></td>
<td style="border: 0px black; background-color: transparent;"></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_MGGA_X_M06_L</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.04416</font></td>
<td class="xl64" style="border: 0px black; background-color: transparent;"></td>
<td style="border: 0px black; background-color: transparent;"></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_MGGA_C_M06_L</font></td>
<td style="border: 0px black; background-color: transparent;"><font face="Calibri">XC_MGGA_X_M06_L</font></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.43029</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.41007</font></td>
<td align="right" class="xl63" style="border: 0px black; background-color: transparent;"><font face="Calibri">53.1</font></td>
</tr>
<tr height="20" style="height: 15pt;">
<td height="20" style="border: 0px black; height: 15pt; background-color: transparent;"><font face="Calibri">XC_HYB_MGGA_XC_MPWB1K</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
<td align="right" style="border: 0px black; background-color: transparent;"><font face="Calibri">0.44</font></td>
<td class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">SCF crash</font></td>
<td align="right" class="xl64" style="border: 0px black; background-color: transparent;"><font face="Calibri">-76.38142</font></td>
<td style="border: 0px black; background-color: transparent;"></td>
</tr>
</tbody></table></div><div> </div><div> </div><div>The MGGA seems to be causing the problems.</div><div> </div><div> </div><div>Thanks, Geoff.</div><div> </div><div> </div><div><br>On Monday, July 28, 2014 11:59:38 AM UTC-4, Geoffrey Wood wrote:</div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;"><div dir="ltr"><div> </div><div>Hello</div><div> </div><div>I am using GAPW and the cutoff is 1000 (not 10000, sorry for the error). I used such a large cutoff because I was testing the convergence behavior of the energy wrt cutoff (I tried 200, 400, 800, 1000 and 4000) in the end I reported the 1000 numbers above as it seemed to give convergence to ~10^-5 Hartree, At this level the difference between the B3LYP and PBE0 energies between the two programs is very small. Knowing that GAPW is a different approximation to using the Gaussian basis sets this difference is to be expected (as reported in the original paper). Using the same settings the difference between the M06 energies seems to be rather large, which is the concern I have. ll try the M06L functional and provide feedback.</div><div> </div><div>Thanks, Geoff.</div><div> </div><div> </div><div> </div><div> </div><div> </div><div> </div><div><br>On Monday, July 28, 2014 10:44:49 AM UTC-4, jgh wrote:</div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;">Hi
<br>
<br>Are you running GPW or GAPW? With that cutoff one never knows.
<br>
<br>M06L might be another good test.
<br>
<br>Juerg
<br>------------------------------<wbr>------------------------------<wbr>--
<br>Juerg Hutter <wbr> Phone : ++41 44 635 4491
<br>Institut für Chemie C FAX : ++41 44 635 6838
<br>Universität Zürich E-<wbr>mail: <a>hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>-----<a>cp...@googlegroups.com</a> wrote: -----
<br>To: <a>cp...@googlegroups.com</a>
<br>From: Geoffrey Wood
<br>Sent by: <a>cp...@googlegroups.com</a>
<br>Date: 07/28/2014 03:27PM
<br>Subject: Re: [CP2K:5563] Re: Geometry optimization by using LIBXC caused no SCF iteration with error NaN
<br>
<br>Hello-
<br>
<br>Thanks for everyone's comments. I discovered that part of the problem I was having was an error in my "arch" file. With that fixed I have compiled and linked CP2K (SVN 14141) libint (1.1.4) and libxc (2.1.0) and tested the XC libraries. These are the results:
<br>
<br>system:
<br>CUTOFF=10000
<br>
<br> O 0.000000 0.000000 -0.065587
<br> H 0.000000 -0.757136 0.520545
<br> H 0.000000 0.757136 0.520545
<br> &KIND H
<br> BASIS_SET 6-31Gxx
<br> POTENTIAL ALL
<br> &END KIND
<br> &KIND O
<br> BASIS_SET 6-31Gxx
<br> POTENTIAL ALL
<br> &END KIND
<br>
<br>
<br>XC_HYB_GGA_XC_PBEH (+25% HF x), CP2k: -76.33577878183063 another popular QM code: -76.3357457
<br>XC_HYB_GGA_XC_B3LYP (+20% HF x), CP2k: -76.41801937267151 another popular QM code: -76.4180549
<br>XC_MGGA_C_M06 (+27% HF x), CP2K: -70.00722600780877
<br>XC_MGGA_X_M06 (+27% HF x) CP2K: -76.02248099520743
<br>XC_MGGA_X_M06 XC_MGGA_C_M06 (+27% HF x): CP2K: -76.42254621238931 another popular QM code: -76.3842635
<br>XC_HYB_MGGA_XC_M06 (+27% HF x): CP2K (crashes on SCF)
<br>XC_HYB_MGGA_XC_M06 (no HF x included): CP2K (crashes on SCF)
<br>
<br>So there still seems to be some problems with M06. I.e. the hybrid-meta XC doesn't seem to work and trying to combine the X and C parts together gives the incorrect energy. Thoughts?
<br>
<br>Geoff
<br>
<br>
<br>
<br>On Friday, July 25, 2014 8:22:08 AM UTC-4, Geoffrey Wood wrote:
<br>Hi Matthias and Michael..
<br>
<br>Thanks for your comments. I have previously run cp2k 2.5.1 with libint 1.1.5 and tested it without any problems i.e. correct energies, gradients 2nd derivatives, BOMD using GGAs, hyrbrid GGAs, with and without libxc. It is only with the current trunk that I'm running into problems. I will re-link with libint 1.1.4 and provide feedback.
<br>
<br>Thanks again.
<br>
<br>Geoff.
<br>
<br>
<br>
<br>On Thursday, July 24, 2014 10:58:59 AM UTC-4, Michael wrote:On Thu, Jul 24, 2014 at 07:43:03AM -0700, Matthias Krack wrote:
<br>> you write that you compiled libint v1.1.5, but CP2K is only validated for
<br>> libint v1.1.4.
<br>
<br>As a data point, Debian/Ubuntu ships libint-1.1.4 and we have not seen
<br>any problems with that, including building cp2k-2.5.1 against it.
<br>
<br>Michael
<br>
<br> --
<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>cp2k+...@googlegroups.com</a>.
<br> To post to this group, send email to <a>cp...@googlegroups.com</a>.
<br> Visit this group at <a onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;" href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/<wbr>group/cp2k</a>.
<br> For more options, visit <a onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;" href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/<wbr>optout</a>.
<br>
<br></blockquote></div></blockquote>