<div>Dear all</div>  <div>Hello</div>  <div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I am studying about relaxation of&nbsp; Na(100) .I use "GGA(PBE)" &nbsp;for XC .I have&nbsp;worked&nbsp;on that issue(relaxation) &nbsp;with k-point mesh&nbsp;equal 32×32×1 ,and&nbsp; Ecut range of 25 to 35Ry&nbsp; for wavefunction and Ecut=140Ry up to 350 Ry&nbsp; for charge density .I&nbsp;run&nbsp;program based on presented ranges&nbsp;,and also used&nbsp;two&nbsp;smearing model(gauss and methfessel-paxton). However,forces&nbsp;had not been&nbsp;reach below 0.0001 Ry/a.u&nbsp; ,but approximately&nbsp;0.0005 Ry/a.u&nbsp;. It should be stated that I reduced dt&nbsp;from 20.DO&nbsp;to 2.D0&nbsp;step by step&nbsp;in relaxation ,but I&nbsp;have not&nbsp;reached forces below 0.0001 Ry/a.u yet.</div>  <div>I should state that I have this difficult for slab with 7 layers and more, that number of atoms in supercell is more than 7 atoms.</div>  <div>I should state that my main question is about reaching
 0.0001 Ry/a.u and what can I do for that?, and I&nbsp;was wondering if you could&nbsp;help me. Thank you for your attention.</div>  <div>Best regards,</div>  <div>Nasser Zare Dehnavi&nbsp;</div><p>&#32;__________________________________________________<br>Do You Yahoo!?<br>Tired of spam?  Yahoo! Mail has the best spam protection around <br>http://mail.yahoo.com