<div dir="ltr">Regarding the response of
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.800000190734863px;background-color:rgb(255,255,255)">Lorenzo Paulatto</span> <div><br></div><div>You are right on both accounts. The calculations involved 54 atoms of Li. For the cp.x I used Vanderbild uspp with 1 valence electron</div>
<div>and for the pw.x, a paw calculated with Holzwarth program with all 3 electrons as valence.</div><div><br></div><div>When I am looking in the literature Jarlborg (Physica Scripta 37, 795 (1988)) gives -14.84 Ry/atom and this is what I got in my pw calculation -800/54=-14.81, but
Dacorogna and Cohen (Phys. Rev. B 34, 4996 (1986) gave -0.617 Ry/atom, which is not close but resembles my result from cp.x -47*2/54=-1.74 Ry/atom.</div><div><br></div><div>Why I am bothered by these values ? I got involved in some EOS calculations and I am trying to solve what is called as first shock</div>
<div>Hugoniot equation (which takes you from an initial state P1, V1, U1 (pressure, volume, internal energy) to a final state P2, T2, U2 by a shock):</div><div><br></div><div> (U2 - U1) = 1/2*(P1+P2)*(V2-V1) if I remember correctly</div>
<div><br></div><div>in which you have to put the initial values. For U I use </div><div> </div><div> U=3/2kT + Etot/N (N=54) </div><div><br></div><div>following some papers by one Lenosky (I don't have references right here). So, if in this equation you use -1.74 Ry or -14.8 Ry it makes a hell of a difference.<br clear="all">
<div><br></div>-- <br>Sylvian (gmail: <a href="mailto:sylviankahane@gmail.com">sylviankahane@gmail.com</a>)<br>
</div></div>