<div dir="ltr">Hi,<div><br></div><div>your CUTOFF is still not converged. For PBE and a system with oxygen atoms to be fully converged in an individual energy you need to be at > 600 Ry, probably somewhat more.</div><div><br></div><div><table border="0" cellpadding="0" cellspacing="0" width="202" style="border-collapse:
 collapse;width:152pt">
 <colgroup><col width="64" style="width:48pt">
 <col width="138" style="mso-width-source:userset;mso-width-alt:5046;width:104pt">
 </colgroup><tbody><tr height="20" style="height:15.0pt">
  <td height="20" align="right" width="64" style="height:15.0pt;width:48pt">300</td>
  <td class="xl63" align="right" width="138" style="width:104pt">-1752.615532</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">350</td>
  <td class="xl63" align="right">-1752.601872</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">400</td>
  <td class="xl63" align="right">-1752.596534</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">450</td>
  <td class="xl63" align="right">-1752.596534</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">500</td>
  <td class="xl63" align="right">-1752.595309</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">550</td>
  <td class="xl63" align="right">-1752.594747</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">600</td>
  <td class="xl63" align="right">-1752.594059</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">700</td>
  <td class="xl63" align="right">-1752.594543</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">800</td>
  <td class="xl63" align="right">-1752.594127</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">900</td>
  <td class="xl63" align="right">-1752.594142</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">1000</td>
  <td class="xl63" align="right">-1752.594152</td>
 </tr>
 <tr height="20" style="height:15.0pt">
  <td height="20" align="right" style="height:15.0pt">1250</td>
  <td class="xl63" align="right">-1752.594154</td>
 </tr></tbody></table><br>Energies vs cutoff for your zeolite system. For a lot of uses you don't need to push that hard, but if you are comparing small relative energies, it might be needed. If the cost of the cutoff is too high, you can also try using the GAPW method or using pseudo potentials with non linear core correction. There are several threads mentioning these on the forum.</div><div><br></div><div>Matt<br><br>On Thursday, November 12, 2015 at 9:23:40 AM UTC, Guanna Li wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>I reported my test results here:<br></div><a href="https://groups.google.com/forum/#!topic/cp2k/UYaBsH6A7e0" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/UYaBsH6A7e0';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/UYaBsH6A7e0';return true;">https://groups.google.com/<wbr>forum/#!topic/cp2k/UYaBsH6A7e0</a><div><br></div><div>All single point calculations used the same input structure. But you can see the enerny didn't converge even the CUTOFF value is big. <br><div> <br><br>On Monday, November 9, 2015 at 4:33:32 PM UTC+1, katarína stančiaková wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2K users,<br>I am running my calculations and I have CUTOFF convergence problem - in fact I can not see any convergence in the total energy. I know that the relative differences in energies are important (not the absolute values of energies), therefore I am comparing the energy differences between two structures (I am working with zeolites, and my two structures differs in the position of hydrogen). Even with very high cutoff values I can not observe any convergence of the total energies as well as energy differences.<br>Do you have any idea where is the problem? I will be thankful for any help!<br><br>Here are my data:<br><br>CUTOFF [Ry]            Total Energy [a.u.]<br>                                structure 1                structure 2<br>600                           -3467.0403965424    -3467.0619253888<br>800                           -3467.0521866655    -3467.0622253178<br>1000                         -3467.0838714739    -3467.1096147317<br>1200                         -3467.0833499728    -3467.1103018678<br>1400                         -3467.0975445683    -3467.1102723936<br><br><br>ENERGY DIFFERENCES<br><br><br>CUTOFF [Ry]    E_str1 - E_str2 [kJ/mol]    <br><br>600                55.9750007029<br>800                26.1004960059<br>1000              66.932470253<br>1200              70.0749269477<br>1400              33.0923456248<br>    <br>Thank you<br>Katarina<br></div></blockquote></div></div></div></blockquote></div></div>