<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Mohnish<div><br><div><div>Il giorno 06/dic/2010, alle ore 07.51, mohnish pandey ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear Vicky!<div> Check the convergence of energy of isolated atoms also. Normally convergence threshold for energy for isolated species is higher than the crystal system. </div></blockquote><div><br></div><div>I'm sorry, but I hope you'll let me disagree on this point. What makes you say this?</div><br><blockquote type="cite"><div>Check for different values of ecut and use the converged value. </div></blockquote><div><br></div><div>One usually uses the same ecutwfc value for both bulk and atomic calculations, in order to take advantage of some error cancellation. Indeed, if you plot the difference between the bulk and atomic total energies as a function of ecutwfc, you should see that it converges much faster than any of the two energies taken alone.</div><br><blockquote type="cite"><div>And next thing is you are doing calculations for bulk system so all the degrees of freedom should be relaxed for geometry optimization so instead of flag "relax" use "vc-relax".. The former one will converge the forces but not the stress and latter one converges the stress also...<br></div></blockquote><div><br></div><div>I agree that the "vc-relax" option is useful, but maybe not so much in this case. I think that bulk W is bcc, hence it can be studied with one atom per cell and there are no internal parameters to relax. You can simply scan the total energy as a function of the lattice parameter to get the equilibrium value, as Vicky already did.</div><div><br></div><div>There are many other issues that come in, especially in the atomic calculation (see recent discussions on the forum). Even when you get a converged DFT value, in some cases this might be considerably far from the experimental value.</div><div><br></div><div>Best regards,</div><div><br></div><div>GS</div><div><br></div><blockquote type="cite"><div>
<div class="gmail_quote">On Sun, Dec 5, 2010 at 11:53 AM, vicky singh <span dir="ltr"><<a href="mailto:kirtinandan07@gmail.com">kirtinandan07@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi<br><br>I am trying to calculate the cohesive energy for tungsten. the code is attached. I am able to get the lattice parameter of 3.19 ang while experimental is 3.16 angs. Then i calculated the cohesive energy of w by calculating the energy of isolated w atom by changing lattice parameter up to 25. I got the cohesive energy of 11.1 ev/atom while reported value is 8.9 ev/atom. This is way too high. can anybody help me. I have varied the K point as well as The E_cut but no advantage.<br>
<font color="#888888">
<br>vicky singh<br>research student<br>Bangalore <br>
</font><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Regards,<br>MOHNISH,<br>-----------------------------------------------------------------<br>Mohnish Pandey<br>Y6927262,5th Year dual degree student,<br>Department of Chemical Engineering,<br>
IIT KANPUR, UP, INDIA<br>-----------------------------------------------------------------<br>
</div>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></body></html>