<div dir="ltr">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Dear Valerie, <br>
    <br>
    you may try a smaller EPS_DEFAULT value like 1.0E-12 or even
    1.0E-14. You may also switch off the smearing, if your systems shows
    a decent band gap. I would guess a cutoff of about 400 Ry should be
    sufficiently accurate.<br>
    <br>
    Matthias<br>
    <br>
    <div>On 24/03/15 17:40, Valérie wrote:<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">Hi all, <br>
        <br>
        I would like to run a geometry optimization of a single molecule
        with Quickstep but got stuck at the first step; i.e. get
        meaningful values of cutoff and relative cutoffs. <br>
        I followed this tutorial:
        <a href="http://www.cp2k.org/howto:converging_cutoff" target="_blank">http://www.cp2k.org/howto:<wbr>converging_cutoff</a> but the total energy
        does not converge to great accuracy with my system:<br>
        <br>
        # Grid cutoff vs total energy<br>
        # REL_CUTOFF = 100<br>
        # Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2
        | NG on grid 3 | NG on grid 4<br>
             50.00   -68.6591152269   51340     490       0       0<br>
            100.00   -71.4960226714   49235    2499      96       0<br>
            150.00   -71.3943596348   46461    4879     490       0<br>
            200.00   -71.4042102566   44831    6209     790       0<br>
            250.00   -71.4043241216   42771    7117    1846      96<br>
            300.00   -71.4043468150   41166    8069    2499      96<br>
            350.00   -71.4042673160   38450   10341    2847     192<br>
            400.00   -71.4043263577   36970   10901    3671     288<br>
            450.00   -71.4043454189   35016   11445    4879     490<br>
            500.00   -71.4043448309   34404   11749    5187     490<br>
            550.00   -71.4043354081   33007   11952    6081     790<br>
            600.00   -71.4043430530   32450   12381    6209     790<br>
        <br>
        I was wondering if there was something I could try to fix this
        or should I carry on with Cutoff = 400?<br>
        <br>
        Many thanks<br>
        Valerie<br>
        <br>
    </div></blockquote>
    <br>
    
  </div>

</div>