<div dir="ltr">Sorry, I have sent the email by accident without finishing it. <div><br></div><div>Attached (in the previous email) are the input and corresponding output files that I have used for the calculations with GAPW and GPW.</div><div><br></div><div>Yours faithfully, </div><div><br></div><div>Lucas Bandeira</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 28, 2021 at 2:45 PM Lucas Bandeira <<a href="mailto:bandeir...@gmail.com">bandeir...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Thank you for your reply, Fabian.<div><br></div><div>But I have performed the same kind of calculations for a Cu bulk, and the forces printed are much larger. I am attaching the input and output files. In addition, I am analyzing cutoff and relative cutoff convergence for a CO molecule. With the GPW, the printed forces were too low, but with the GAPW the forces are much higher. Could you help me with that? </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 28, 2021 at 2:09 PM Fabian Ducry <<a href="mailto:fabia...@gmail.com" target="_blank">fabia...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  
    
  
  <div>
    <p>Hi Lucas</p>
    <p>The forces are zero because of symmetry (the 1E-13 Ha is noise).
      You should check the pressure (stress tensor) instead.
      Alternatively, you can perform a volume optimization and check the
      convergence of the lattice parameter wrt. k-point grid.</p>
    <p>Cheers Fabian<br>
    </p>
    <div>On 28.06.2021 19:00, Lucas Bandeira
      wrote:<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">Dear cp2k community,
        <div><br>
        </div>
        <div>I am trying to analyze k-points convergence for a Au bulk,
          but I am facing a problem: the atomic forces are too small.
          They are of the order of 1E-13 Ha. Then, when I try to compute
          the atomic forces between two consecutive grids, they are of
          the order of 1E-14 Ha, much lower than the threshold I
          use, 4.5E-4 Ha (the default maximum atomic force used in cp2k
          geometry optimization). So basically it means that the results
          will be the same no matter the k-point mesh I choose, from a
          1x1x1 to a 15x15x15 mesh. Can somebody help me to solve this
          problem? Attached are one of the inputs that I have used and
          the corresponding output.</div>
        <div><br>
        </div>
        <div>Yours sincerely,</div>
        <div><br>
        </div>
        <div>Lucas Bandeira</div>
      </div>
      -- <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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
      To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAM1%3DHZm0WxJwKD4LrhCCER75TMOYryC6Ka685NydKzP6S8eyQQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAM1%3DHZm0WxJwKD4LrhCCER75TMOYryC6Ka685NydKzP6S8eyQQ%40mail.gmail.com</a>.<br>
    </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/80e6f8fc-f5a9-e2b8-80c2-406eeac7b1b2%40gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/80e6f8fc-f5a9-e2b8-80c2-406eeac7b1b2%40gmail.com</a>.<br>
</blockquote></div>
</blockquote></div>